Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi526,7960CG31.1% 112.0 / 62.9 74A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base G (10/13);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTAT  >  NC_000913_3_pae_pgi/526649‑526944
                                                                                                                                                   |                                                                                                                                                    
tttttcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCAGCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                     <  2:783362/144‑1 (MQ=255)
ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                     <  2:653127/149‑1 (MQ=255)
    gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGagt                                                                                                                                                 <  1:716208/149‑3 (MQ=255)
      aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCCCTGCGACCACCGGGGACTGCCCCTGGCCCTTATCCGCGAACACGGCCATACGGCGCGGCGGCCGGCatat                                                                                                                                               >  2:149505/1‑149 (MQ=14)
           gggCTGACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATAtgatga                                                                                                                                          <  2:233458/149‑1 (MQ=255)
            ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGGCCACCGGGGACTGCCGCGGGCGCTTATCAGCGACGACGGCAATACGGCGTGGCGCGGGGAGCATGGTGaa                                                                                                                                         >  2:479669/1‑149 (MQ=11)
             gCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATAACGCGTGGCGCGGGGCGTATGATGAAt                                                                                                                                        >  1:731032/1‑149 (MQ=11)
              cTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTTGCGCGGGGAGTATGATGAATg                                                                                                                                       <  2:694159/149‑1 (MQ=11)
               ggACGGTGGCGCAACGGGCCAGACGGGTGGAGGCGGAATACAGACCGGCGCGAAAAGCTCATCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                                      <  2:488590/148‑1 (MQ=11)
                gatcgTGGGGTAACTGGCCAGACATGTGGAGCCGGAATACACACCGGCGCGAAAATCTCATCTGTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCATCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                     <  1:197880/145‑1 (MQ=11)
                gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAGAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                     <  2:669922/149‑1 (MQ=11)
                gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                     <  1:187414/149‑1 (MQ=255)
                     tgtatcAACTGGCCAGACAGGTGGAGCCGGAATTCACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                <  1:493903/144‑1 (MQ=255)
                     tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                                <  1:52438/149‑1 (MQ=255)
                     tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                <  2:707343/149‑1 (MQ=255)
                        aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGTGCGAAGAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAg                                                                                                                             <  2:561283/149‑1 (MQ=255)
                         gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGc                                                                                                                            <  1:176682/149‑1 (MQ=255)
                            aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTa                                                                                                                         >  1:244394/1‑149 (MQ=255)
                                 ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAGCAGCTTAATgag                                                                                                                    >  1:431809/1‑149 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCaa                                                                                                                                 >  1:445430/1‑135 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCaa                                                                                                                                 <  2:445430/135‑1 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTACTgagg                                                                                                                   >  1:47398/1‑149 (MQ=11)
                                   aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                  >  1:194208/1‑149 (MQ=255)
                                   aGACAGGGGGGGGCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCATCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGGGCGCTGGATATGATGCATGGGGCAACCAGCTTAATgagga                                                                                                                  <  2:723617/149‑1 (MQ=17)
                                       aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCCGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCCGCTTAATGAGGAGAAc                                                                                                              >  1:317611/1‑149 (MQ=11)
                                         gtgaggCCGGAATACACAACGGCGCGAAAAGCTCATCTTGATCACTGCGACCACCCGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATCCGGCGTGGCGCGGGGCGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                            <  1:381352/144‑1 (MQ=11)
                                         gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCACCGAAGACGGCAATACGCCCTGGCCCGGGGGGGCTGATCACTAGGGCACACAGCCTCATGAGGAGaaaca                                                                                                            >  2:117804/1‑146 (MQ=11)
                                         gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGAACACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCGTAATGAGGAGAAccc                                                                                                            >  1:407159/1‑149 (MQ=255)
                                          tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                           >  2:764907/1‑149 (MQ=255)
                                          tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAAACGGAGACGGGCAGGGCGGAGTATGAATGGGGCAACAAGCTTAATGAGGAGaaacag                                                                                                           >  2:478928/1‑145 (MQ=14)
                                           ggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                          <  1:764907/149‑1 (MQ=255)
                                            gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATg                                                                                                                                       >  2:211313/1‑119 (MQ=255)
                                            gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATg                                                                                                                                       <  1:211313/119‑1 (MQ=255)
                                            gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                         <  1:558308/149‑1 (MQ=255)
                                                  gAATACCCACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                                   <  2:203540/149‑1 (MQ=255)
                                                  gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacct                                                                                                   <  2:207231/149‑4 (MQ=11)
                                                  gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                                   <  2:705895/149‑1 (MQ=255)
                                                  gAAGACACACCGGCGCGAAAAGCTCATCTTGATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                                   <  2:34062/149‑1 (MQ=255)
                                                      acacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctgcac                                                                                               >  1:507395/1‑142 (MQ=1)
                                                        acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                             >  1:711457/1‑149 (MQ=255)
                                                        acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                             >  1:447732/1‑149 (MQ=255)
                                                                  gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACGCGCTACATGTGTATTAGCCGGACCgc                                                                                   >  1:656430/1‑148 (MQ=18)
                                                                  gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                                 >  2:459438/1‑135 (MQ=255)
                                                                  gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                                 <  1:459438/135‑1 (MQ=255)
                                                                  gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                                   <  1:81352/149‑1 (MQ=255)
                                                                  gAAAAGCTCATCTTTATCACTGCGAAAACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                                   >  2:486836/1‑149 (MQ=255)
                                                                        cTCATCTTGATCACTGCCAGCACCGGGGACGGCCGCAGGCGGTTATGGGCTAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                             <  2:83039/149‑1 (MQ=18)
                                                                         tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGCGAAATATGACGAATGGGGCAACCAGCTTGTTGGGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                            >  2:716813/1‑149 (MQ=17)
                                                                         tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                            >  1:187827/1‑149 (MQ=255)
                                                                               ttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagc                                                                      >  1:20283/1‑149 (MQ=255)
                                                                                     cTGCGACCACCGGGGACTGCCGCTGGCGCTGATGATCGAAGGCGGCAATACGGCGTGGCGCCGGGGTTATGATGGACGGTGCAACCAGCTTAATTCCGATAACCGGCATCACCTGCGCTATCCGTACCGCCCGCCCGGGCAGCCGCatg                                                                >  2:41052/1‑149 (MQ=255)
                                                                                      tGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCctgc                                                               >  1:86287/1‑148 (MQ=255)
                                                                                      tGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatga                                                               >  1:213517/1‑149 (MQ=255)
                                                                                         gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                            <  2:187827/149‑1 (MQ=255)
                                                                                          accaccGGGGACTGTCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                           <  1:197423/149‑1 (MQ=255)
                                                                                          accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                           <  1:588432/149‑1 (MQ=255)
                                                                                                 gggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg                                                    >  1:405925/1‑149 (MQ=255)
                                                                                                  ggACTGCCGCTGGCGCTTAGCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAGCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAgg                                                   <  2:381694/149‑1 (MQ=255)
                                                                                                   gACCGCCGGTGGCGCTTATCAGAGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                                  <  2:641355/149‑1 (MQ=255)
                                                                                                        cgcgtGGCGCTTATCAGCGAAGACGGCAATGCGGCGTGGGGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                             <  2:257721/145‑1 (MQ=255)
                                                                                                        ccGCTGGGGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCAGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                             <  2:460654/149‑1 (MQ=255)
                                                                                                         cGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTAc                                            <  1:236896/149‑1 (MQ=255)
                                                                                                            tGGCGCTTATCAGCGAAGACGGCAATACGGAGAGGCGCGCGTAATGTGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACtat                                         >  1:189268/1‑149 (MQ=17)
                                                                                                               cgcTTATCAGCGAAGACGGCAATACGGCGGGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAc                                      <  2:586053/149‑1 (MQ=255)
                                                                                                                 cTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCg                                    >  2:165873/1‑149 (MQ=255)
                                                                                                                  ttcttcgCGAAGACGTCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGt                                   >  2:405749/7‑149 (MQ=255)
                                                                                                                     tCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg                                                          >  2:701570/1‑123 (MQ=255)
                                                                                                                     tCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg                                                          <  1:701570/123‑1 (MQ=255)
                                                                                                                       aGCGAAGACGGCAATACGGCGTGGAGAGCGGAATATGATGAATGGGGCAACCAGCGTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCTTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCg                              <  2:149442/149‑1 (MQ=255)
                                                                                                                             gACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTAc                                            >  2:541323/1‑129 (MQ=255)
                                                                                                                             gACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTAc                                            <  1:541323/129‑1 (MQ=255)
                                                                                                                              aCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagca                                                                     <  1:620887/103‑1 (MQ=255)
                                                                                                                              aCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagca                                                                     >  2:620887/1‑103 (MQ=255)
                                                                                                                                ggCAATACGGCGTGGAGCGCGGAAAATGATTAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGa                     >  2:289734/1‑149 (MQ=255)
                                                                                                                                     tACGGCGTGGAGCGGGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGt                <  2:711457/149‑1 (MQ=255)
                                                                                                                                         gCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCa            <  2:405925/149‑1 (MQ=255)
                                                                                                                                                  gcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGta   >  2:376873/1‑149 (MQ=255)
                                                                                                                                                   cgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtat  >  1:358832/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTAT  >  NC_000913_3_pae_pgi/526649‑526944

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑