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.4% 52891.2 / inf 13621V331A (GTG→GCG) lacIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (182/177);  new base G (6557/6700);  total (6741/6880)
Fisher's exact test for biased strand distribution p-value = 6.69e-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 13625 total aligned reads displayed.

CTGATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATA  >  CP009273/360604‑360900
                                                                                                                                                    |                                                                                                                                                    
   gataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCa                                                                                                                                                   >  1:1311024‑M1/46‑149 (MQ=255)
    ataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAg                                                                                                                                                  >  2:3607228‑M1/45‑149 (MQ=255)
    ataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAg                                                                                                                                                  <  1:494278‑M1/105‑1 (MQ=255)
      aagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAggg                                                                                                                                                >  1:1000016‑M1/43‑149 (MQ=255)
      aagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAggg                                                                                                                                                >  2:2494122‑M1/43‑149 (MQ=255)
       agctgtcaaaccagatcacttcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggt                                                                                                                                               >  1:3610151‑M1/42‑149 (MQ=255)
        gctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtg                                                                                                                                              <  1:2263162‑M1/109‑1 (MQ=255)
          tgtcaaaccagatcaattcgcgctaacttacattaagtGCGGTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtggt                                                                                                                                            >  2:1146162‑M1/39‑149 (MQ=255)
           gtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtt                                                                                                                                           >  1:3515698‑M1/38‑149 (MQ=255)
           gtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtt                                                                                                                                           >  1:3799931‑M1/38‑149 (MQ=255)
            tcaaaccagatcaaatcgcgctaacttacataaaatGCGTAGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttt                                                                                                                                          <  2:1812864‑M1/113‑1 (MQ=255)
             caaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtttt                                                                                                                                         <  2:1915946‑M1/114‑1 (MQ=255)
              aaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCACCGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttttt                                                                                                                                        >  2:4364082‑M1/35‑149 (MQ=255)
              aaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttttt                                                                                                                                        <  1:1426037‑M1/115‑1 (MQ=255)
               aaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                                                                                       >  1:4733221‑M1/34‑149 (MQ=255)
               aaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                                                                                       <  2:734890‑M1/116‑1 (MQ=255)
               aaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACACGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                                                                                       <  2:3707074‑M1/116‑1 (MQ=255)
                accagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                                                                                                      >  1:2177077‑M1/33‑149 (MQ=255)
                accagatcaagtcgcgctaacttaccttaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                                                                                                      >  2:630859‑M1/33‑149 (MQ=255)
                 ccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctt                                                                                                                                     >  2:3890765‑M1/32‑149 (MQ=255)
                 ccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctt                                                                                                                                     <  2:4556733‑M1/118‑1 (MQ=255)
                  cagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttt                                                                                                                                    <  1:1907172‑M1/119‑1 (MQ=255)
                       caattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAcc                                                                                                                               <  2:2584711‑M1/124‑1 (MQ=255)
                       caattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAcc                                                                                                                               >  1:3702117‑M1/26‑149 (MQ=255)
                        aattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                              <  2:3161066‑M1/125‑1 (MQ=255)
                        aattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                              <  1:2096230‑M1/125‑1 (MQ=255)
                           tcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTg                                                                                                                           >  1:690784‑M1/22‑149 (MQ=255)
                               gctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                       >  1:3979732‑M1/18‑149 (MQ=255)
                               gctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                       <  2:2278496‑M1/132‑1 (MQ=255)
                               gctaacttacattaattGCGTTGCGCTCACTGCCCGATTTCCAGTCGTGAAACCTGTCGTGCAAGCTGCATTATTGAATCGGACAACGCGCGGGGAGAGTCGGTTTGCGTATTGGGCGCAAGGGTGGATTTTCTTTTCACCAGTGAGAc                                                                                                                       >  2:5015870‑M1/18‑149 (MQ=255)
                                 taacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                     <  1:4341073‑M1/134‑1 (MQ=255)
                                     ttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCaa                                                                                                                 <  2:1639528‑M1/138‑1 (MQ=255)
                                      tacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                                                                                                >  2:2970119‑M1/11‑149 (MQ=255)
                                         cAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                             >  2:825124/1‑149 (MQ=255)
                                         attaattGCGGTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                             >  2:4915167‑M1/8‑149 (MQ=255)
                                          ttaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTAGTGGGCGCCAGGGTGGTTTTTATTTTAACCAGTGAGACGGGCCACAGCt                                                                                                            <  1:960541‑M1/143‑1 (MQ=255)
                                            aattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGa                                                                                                          >  2:1230626‑M1/5‑149 (MQ=255)
                                               tGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTg                                                                                                       >  2:4089485‑M1/2‑149 (MQ=255)
                                                  gTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGccc                                                                                                    >  1:2447821/1‑149 (MQ=255)
                                                  gTTGCGCTCACTGCCCGCTGTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGccc                                                                                                    <  2:3246328/149‑1 (MQ=255)
                                                   ttGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCt                                                                                                   >  2:2885885/1‑149 (MQ=255)
                                                     gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTc                                                                                                 >  2:343486/1‑149 (MQ=255)
                                                        cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCCGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAcc                                                                                              <  1:2723648/149‑1 (MQ=255)
                                                             tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCATCACCCCCTg                                                                                         >  1:4157288/1‑149 (MQ=255)
                                                              gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTgg                                                                                        >  1:3430575/1‑149 (MQ=255)
                                                                 cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGccc                                                                                     >  2:353034/1‑149 (MQ=255)
                                                                   atttCCAGTCGGGAAACCTGTCGTGACAGCTGCATTAATGAATCGGACAAAGCGCGGGGAGAGGCTGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTAAAAAGTGAGACGGGCAACAGCTGAATGACAATCAACGCCTGTCCCtt                                                                                   >  2:5105630/2‑148 (MQ=255)
                                                                        cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagaga                                                                              <  1:3135311/149‑1 (MQ=255)
                                                                          gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGt                                                                            <  2:5184312/149‑1 (MQ=255)
                                                                           tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGtt                                                                           <  1:599334/149‑1 (MQ=255)
                                                                              ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                        >  2:4640398/1‑149 (MQ=255)
                                                                              ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                        <  2:902810/149‑1 (MQ=255)
                                                                               gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTAACCCGTGAGACGGGCAACAGCTGATTTCCCATCACCGCCGGGCCCTGAGAGAGTTgcag                                                                       >  1:2347225/1‑149 (MQ=255)
                                                                                 aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                     <  2:2008095/149‑1 (MQ=255)
                                                                                  aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                    >  2:4456867/1‑149 (MQ=255)
                                                                                   ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg                                                                   >  2:1720165/1‑149 (MQ=255)
                                                                                       tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                               <  1:2917768/149‑1 (MQ=255)
                                                                                        cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTc                                                              <  2:2169774/149‑1 (MQ=255)
                                                                                             cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGc                                                         <  2:2765469/149‑1 (MQ=255)
                                                                                             cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGc                                                         >  1:1523247/1‑149 (MQ=255)
                                                                                             cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGATTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGc                                                         <  2:1442115/149‑1 (MQ=255)
                                                                                              aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCGCCGGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                        >  1:179659/1‑149 (MQ=255)
                                                                                              aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                        <  1:4256635/149‑1 (MQ=255)
                                                                                              aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTAAACAGTGAGACGGGCAACATCTGATTTCCCGTCACCCCCAGGCCCAGAGATAGTTGCAGCAAGCGGTCCACGca                                                        >  1:4849473/1‑148 (MQ=255)
                                                                                                 tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTATTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGt                                                     <  1:4520737/149‑1 (MQ=255)
                                                                                                   cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGttt                                                   <  1:4447193/149‑1 (MQ=255)
                                                                                                   cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCGGAGAGAGTTGCAGCAAGCGGTCCACGCTGGttt                                                   >  2:1050904/1‑149 (MQ=255)
                                                                                                   cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAACAGTGAGACGGGCAACAGCTGATTTCCCATCACCGCCAGGCCCGGAGAGAGTTGCAGCAAGCGGTCCACGCTGGttt                                                   >  1:1108280/1‑149 (MQ=255)
                                                                                                     ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGc                                                 >  1:874849/1‑149 (MQ=255)
                                                                                                      tAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcc                                                >  1:2619729/1‑149 (MQ=255)
                                                                                                        aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcccc                                              <  2:4256739/149‑1 (MQ=255)
                                                                                                        aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcccc                                              <  1:1705635/149‑1 (MQ=255)
                                                                                                         tGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTAACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTTCAGCAAGCGGTCCACGCTGGTTTGCCCca                                             >  1:4626276/1‑149 (MQ=255)
                                                                                                          gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCAGGCCCGGAGAGAGTTGCAGCAAGCGGTCCACGCGGGGTTGcccccg                                            >  1:5032272/1‑147 (MQ=255)
                                                                                                            aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                          <  1:2528439/149‑1 (MQ=255)
                                                                                                             tCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTAACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                         >  1:281798/1‑149 (MQ=255)
                                                                                                                   aaCGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                   >  2:544702/1‑149 (MQ=255)
                                                                                                                   aaCGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                   >  2:4720806/1‑149 (MQ=255)
                                                                                                                   aaCGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCAGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                   >  2:4824843/1‑149 (MQ=255)
                                                                                                                       cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATcc                               <  1:1613403/149‑1 (MQ=255)
                                                                                                                       cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATcc                               <  2:1538212/149‑1 (MQ=255)
                                                                                                                       cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATcc                               <  1:1212626/149‑1 (MQ=255)
                                                                                                                        gcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                              <  1:5225988/149‑1 (MQ=255)
                                                                                                                         cgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTg                             <  2:2381493/149‑1 (MQ=255)
                                                                                                                         cgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTg                             <  2:2098895/149‑1 (MQ=255)
                                                                                                                          ggggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCAGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGt                            >  1:799693/1‑149 (MQ=255)
                                                                                                                           gggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                           >  2:3794490/1‑149 (MQ=255)
                                                                                                                           gggAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                           <  1:3015547/149‑1 (MQ=255)
                                                                                                                            ggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGttt                          <  2:267635/149‑1 (MQ=255)
                                                                                                                             gagaGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTATTTTAACCACTGAGACGGGCAACACCTGATTGCCCATCACCGCCAGGCCCTGAGAGAGTTGCAGCAAGCCGTCCACCCAGGTTTGCCCCACCACGCGAAAATCCAGTTTg                         >  1:5122738/1‑149 (MQ=255)
                                                                                                                              agagGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa                        <  2:1343571/149‑1 (MQ=255)
                                                                                                                                 ggCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                     >  2:399018/1‑149 (MQ=255)
                                                                                                                                 ggCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                     >  2:3336293/1‑149 (MQ=255)
                                                                                                                                   cGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                   <  1:2007651/149‑1 (MQ=255)
                                                                                                                                    ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                  <  1:3245776/149‑1 (MQ=255)
                                                                                                                                        tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa              <  1:3876439/149‑1 (MQ=255)
                                                                                                                                         gCGTATTGGGCGCCAGGGTGGTTTTTATTTTCACCAGTGAGAAGGGCAACAGCATATTTCCCATCACCCCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACCCTGGTTTTCCCCAGCAGGCGAAAATCCTGTTTGATTGTGGTTaaa             >  1:4071700/1‑148 (MQ=255)
                                                                                                                                             aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg         >  2:3066415/1‑149 (MQ=255)
                                                                                                                                             aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCAGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAActgcg         >  2:4167106/1‑149 (MQ=255)
                                                                                                                                               tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg       <  1:3340789/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATA  >  CP009273/360604‑360900

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