Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,752 A→G 100% V331A (GTG→GCG)  lacI ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273360,7520AG97.7% 32075.0 / inf 8115V331A (GTG→GCG) lacIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (88/98);  new base G (3787/4137);  total (3877/4238)
Fisher's exact test for biased strand distribution p-value = 9.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 8116 total aligned reads displayed.

GATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATAT  >  CP009273/360606‑360899
                                                                                                                                                  |                                                                                                                                                   
  ataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCa                                                                                                                                                  >  1:323484‑M1/45‑148 (MQ=255)
  ataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCa                                                                                                                                                  <  2:3294543‑M1/104‑1 (MQ=255)
  ataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCa                                                                                                                                                  <  1:2675220‑M1/104‑1 (MQ=255)
      gctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggt                                                                                                                                              <  1:1503726‑M1/108‑1 (MQ=255)
      gctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggt                                                                                                                                              >  2:3361128‑M1/41‑148 (MQ=255)
      gctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggt                                                                                                                                              >  1:2119025‑M1/41‑148 (MQ=255)
       ctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtg                                                                                                                                             >  2:1130541‑M1/40‑148 (MQ=255)
       ctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAACGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtg                                                                                                                                             >  2:156337‑M1/40‑148 (MQ=255)
         gtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtggt                                                                                                                                           >  2:1442256‑M1/38‑148 (MQ=255)
           caaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttt                                                                                                                                         >  2:3419320‑M1/36‑148 (MQ=255)
            aaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtttt                                                                                                                                        >  2:1505311‑M1/35‑148 (MQ=255)
             aaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttttt                                                                                                                                       <  1:1074025‑M1/115‑1 (MQ=255)
             aaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttttt                                                                                                                                       >  2:1691739‑M1/34‑148 (MQ=255)
              accagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                                                                                      >  1:2463514‑M1/33‑148 (MQ=255)
              accagataaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                                                                                      >  2:1988976‑M1/33‑148 (MQ=255)
                cagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctt                                                                                                                                    <  1:1768611‑M1/118‑1 (MQ=255)
                cagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAAAGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTttttctt                                                                                                                                    <  2:1803313‑M1/118‑1 (MQ=255)
                cagataaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctt                                                                                                                                    <  2:2317302‑M1/118‑1 (MQ=255)
                 agatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttt                                                                                                                                   <  1:57538‑M1/119‑1 (MQ=255)
                 agatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttt                                                                                                                                   >  1:1199182‑M1/30‑148 (MQ=255)
                  gatcaattcgcgctaacttacatgacttGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctttt                                                                                                                                  >  2:2246348‑M1/29‑148 (MQ=255)
                     caattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                                                                                               >  2:1746761‑M1/26‑148 (MQ=255)
                     aaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                                                                                               >  2:2769419‑M1/26‑148 (MQ=255)
                        ttcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAg                                                                                                                            >  2:2915612‑M1/23‑148 (MQ=255)
                         tcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCTGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGCGCCAGGGTGGTTTTTCTTTtaacccct                                                                                                                           <  1:16443‑M1/127‑4 (MQ=255)
                            cgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                        >  1:2528691‑M1/19‑148 (MQ=255)
                             gctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                       <  2:271146‑M1/131‑1 (MQ=255)
                              ctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                      <  2:836984‑M1/132‑1 (MQ=255)
                                 acttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACggg                                                                                                                   <  2:1007802‑M1/135‑1 (MQ=255)
                                  cttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGc                                                                                                                  <  1:1383026‑M1/136‑1 (MQ=255)
                                     acattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                                                                                               <  1:710953‑M1/139‑1 (MQ=255)
                                         taattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCTCCAGTGAAACGGGCAACAGCt                                                                                                           <  1:2193695‑M1/143‑1 (MQ=255)
                                         taattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCt                                                                                                           <  1:2916739‑M1/143‑1 (MQ=255)
                                           attGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGa                                                                                                         <  2:2455482‑M1/145‑1 (MQ=255)
                                           attGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGa                                                                                                         >  1:522188‑M1/4‑148 (MQ=255)
                                           attGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTGTTTGTTTTAACCAGTGAGACGGGCCACACCTGa                                                                                                         <  1:1258812‑M1/145‑1 (MQ=255)
                                            ttGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAt                                                                                                        <  2:2178563‑M1/146‑1 (MQ=255)
                                            ttGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGCGTGGGCGCCAGGGTGGTTTTTATTTTCACCAGTTAGACCGGCAACAGCTGAt                                                                                                        <  2:1632183‑M1/146‑1 (MQ=255)
                                            ttGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCAGTAATGAATCGGCCAACGCGCCGGGAGAGGCGGTTTGCGGAGTGGGCGCCAGGGTGGGTTTTATTTTCACCAGTGTGACGGGGAACAGCTgaa                                                                                                        <  1:2866695‑M1/146‑2 (MQ=255)
                                              gcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTg                                                                                                      <  2:2973256/148‑1 (MQ=255)
                                               cgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGc                                                                                                     >  2:887908/1‑148 (MQ=255)
                                               cgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGc                                                                                                     <  1:1973229/148‑1 (MQ=255)
                                                   gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCtt                                                                                                 <  2:1865152/148‑1 (MQ=255)
                                                       tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAcc                                                                                             >  2:1077218/1‑148 (MQ=255)
                                                         aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGc                                                                                           <  2:2054508/148‑1 (MQ=255)
                                                         aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGc                                                                                           <  1:3518890/148‑1 (MQ=255)
                                                         aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGc                                                                                           <  2:1316510/148‑1 (MQ=255)
                                                          cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGGGAGACGGGCAACAGCTGATTGGCCGTCACCGcc                                                                                          >  1:1849202/1‑148 (MQ=255)
                                                             cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTgg                                                                                       >  2:1258614/1‑148 (MQ=255)
                                                              ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTGTTTTAACCAGTGAGACGGGCAACAGCTGATTTCCCGTCACCCCCGGGc                                                                                      >  1:23964/1‑148 (MQ=255)
                                                                  tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTg                                                                                  <  1:3301894/148‑1 (MQ=255)
                                                                      cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagag                                                                              <  2:78271/148‑1 (MQ=255)
                                                                       aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagaga                                                                             <  2:3475304/148‑1 (MQ=255)
                                                                        gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGagagag                                                                            <  2:1372579/148‑1 (MQ=255)
                                                                         tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGt                                                                           <  1:1906570/148‑1 (MQ=255)
                                                                            ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgc                                                                        >  1:3165101/1‑148 (MQ=255)
                                                                               aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagc                                                                     >  2:1927435/1‑148 (MQ=255)
                                                                                    gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                <  2:2652975/148‑1 (MQ=255)
                                                                                     tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCgg                                                               <  1:1001059/148‑1 (MQ=255)
                                                                                      cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                              <  2:1201133/148‑1 (MQ=255)
                                                                                      cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                              <  2:2388041/148‑1 (MQ=255)
                                                                                        tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACATCAGATTGCCCTACACCCCCTTGCCCTGAGAGAGTTGCAGCAAGCGGTcc                                                            >  1:1435349/1‑148 (MQ=255)
                                                                                        tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCGCCGGGCCCTGAGAGAGGTGCAGCAAGCGGTcc                                                            >  1:642194/1‑148 (MQ=255)
                                                                                        tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTcc                                                            <  1:776826/148‑1 (MQ=255)
                                                                                            aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTCACCAGTGAGACGGGCAACACCTGATTGGCCTTAACCGCCTGGCCCGGAGAGAGTTGGAGCAAGCGGTCCACGc                                                        >  1:1312856/1‑148 (MQ=255)
                                                                                             gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGTGAGAGTTGCAGCAAGCGGTCCACGCt                                                       >  2:938014/1‑148 (MQ=255)
                                                                                                gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGt                                                    <  1:837683/148‑1 (MQ=255)
                                                                                                gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACCCTGGt                                                    >  1:1709698/1‑148 (MQ=255)
                                                                                                   ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTg                                                 >  2:516623/1‑148 (MQ=255)
                                                                                                      aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGccc                                              <  2:406786/148‑1 (MQ=255)
                                                                                                       tGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCCCCTGGCCCTGAGAGAGTTGCAGCAACCGGTCCACGCTGGTTTGcccc                                             >  1:3107118/1‑148 (MQ=255)
                                                                                                       tGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcccc                                             <  2:45762/148‑1 (MQ=255)
                                                                                                        gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCGGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGccccc                                            >  1:2807672/1‑147 (MQ=255)
                                                                                                         aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcag                                           >  2:2524600/1‑148 (MQ=255)
                                                                                                              gCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                      >  1:2050882/1‑148 (MQ=255)
                                                                                                               ccAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTGTTTTGACCACTGAGACCGGGCACCCCGGTTTTCCCATCACCGCCGGGCCCCGTGAGAGTTGCAGCAAGCGGTCCACCCTGGGTTTGCCCACCAGGCg                                     >  1:945474/1‑148 (MQ=255)
                                                                                                                  aCGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                  >  1:890081/1‑148 (MQ=255)
                                                                                                                  aCGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                  <  2:578349/148‑1 (MQ=255)
                                                                                                                   cgcgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaa                                 >  1:3442855/1‑148 (MQ=255)
                                                                                                                   cgcgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTAACCAGTGAGACGGGCCACAGCAGATTTCCCTTAACCCCCAGGCCCGGAGATAGTTTCAGCAAGCCGTCCACGCTGGTTTTCCCCACCAGGCGaaaa                                 >  1:387805/1‑148 (MQ=255)
                                                                                                                   cgcgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTGTTTCTTTTCACCACTGAGACCGGCAACACCTGATTTCCCATAACCCCCGGGCCCGGAGTGAGTTTCAGCAAGCGGTCCACCCTGGTTTGCCCCAGCAGGCGcaaa                                 >  1:1580574/1‑148 (MQ=255)
                                                                                                                     cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACTCTGGTTTGCCCCAGCAGGCGAAAATc                               <  2:3163889/148‑1 (MQ=255)
                                                                                                                     cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATc                               <  2:654918/148‑1 (MQ=255)
                                                                                                                     cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATc                               <  2:3543420/148‑1 (MQ=255)
                                                                                                                      gcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTGTTTCTTTTCACCAGTGAGAAGGGCAACAGCTGATTGCCCATCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCCGTCCACCCTGGTTTGCCCCAGCAGGCGAAAATcc                              >  1:2259800/1‑148 (MQ=255)
                                                                                                                        ggggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTg                            <  1:232722/148‑1 (MQ=255)
                                                                                                                         gggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGt                           <  2:727966/148‑1 (MQ=255)
                                                                                                                         gggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGt                           <  1:455214/148‑1 (MQ=255)
                                                                                                                          ggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                          >  2:2846723/1‑148 (MQ=255)
                                                                                                                            agagGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg                        <  2:3044617/148‑1 (MQ=255)
                                                                                                                             gagGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa                       <  2:2596006/148‑1 (MQ=255)
                                                                                                                              agGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                      >  2:3227607/1‑148 (MQ=255)
                                                                                                                                  ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                  <  1:2603058/148‑1 (MQ=255)
                                                                                                                                   gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                 <  1:3042238/148‑1 (MQ=255)
                                                                                                                                   gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                 <  2:211721/148‑1 (MQ=255)
                                                                                                                                   gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCATCACCGCCAGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCACCAGGCGAAAATCCTGTTTGAtggtgg                 >  1:1132382/1‑148 (MQ=255)
                                                                                                                                   gTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTAAACACTGAGACGGGCCACACCGGATTGCCCTTCACCCCCGGGCCCTGAGCGAGTTGCAGCAACCGGTCCACCCAGGTTTGCCCCAGCAGGCGAAAATCCAGTTTGAtggtgg                 >  1:337496/1‑148 (MQ=255)
                                                                                                                                             tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg       <  1:2752705/148‑1 (MQ=255)
                                                                                                                                               ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa     >  2:3100795/1‑148 (MQ=255)
                                                                                                                                               ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGATGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa     <  1:2981124/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATAT  >  CP009273/360606‑360899

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: