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,7520AG98.4% 70843.6 / inf 17855V331A (GTG→GCG) lacIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (130/159);  new base G (8410/9147);  total (8546/9309)
Fisher's exact test for biased strand distribution p-value = 3.42e-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 17855 total aligned reads displayed.

TGATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATA  >  CP009273/360605‑360898
                                                                                                                                                   |                                                                                                                                                  
   ataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAg                                                                                                                                                >  2:5016250‑M1/45‑149 (MQ=255)
    taagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAgg                                                                                                                                               >  1:3816482‑M1/44‑149 (MQ=255)
     aagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAggg                                                                                                                                              <  1:2286031‑M1/107‑1 (MQ=255)
      agctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggt                                                                                                                                             <  1:1656375‑M1/108‑1 (MQ=255)
        ctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtgg                                                                                                                                           >  1:2770285‑M1/40‑149 (MQ=255)
         tgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtggt                                                                                                                                          <  2:3634207‑M1/111‑1 (MQ=255)
          gtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtt                                                                                                                                         <  1:793796‑M1/112‑1 (MQ=255)
               accagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                                                                                                    >  1:713610‑M1/33‑149 (MQ=255)
                ccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctt                                                                                                                                   >  1:1174005‑M1/32‑149 (MQ=255)
                 cagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttt                                                                                                                                  <  2:4467313‑M1/119‑1 (MQ=255)
                  agatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctttt                                                                                                                                 >  2:164864‑M1/30‑149 (MQ=255)
                   gatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                                                                                                <  2:2668743‑M1/121‑1 (MQ=255)
                    atcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCa                                                                                                                               >  2:4185226‑M1/28‑149 (MQ=255)
                     tcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                                                                                              <  1:4463855‑M1/123‑1 (MQ=255)
                       aattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                            <  2:2301725‑M1/125‑1 (MQ=255)
                       aattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                            <  2:3457668‑M1/125‑1 (MQ=255)
                         ttcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                          <  2:3810438‑M1/127‑1 (MQ=255)
                           cgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                        <  2:2855404‑M1/129‑1 (MQ=255)
                           cgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                        >  1:3644616‑M1/21‑149 (MQ=255)
                           cgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTCACCAGTga                                                                                                                        <  2:3993991‑M1/129‑1 (MQ=255)
                            gcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                       >  2:1191455‑M1/20‑149 (MQ=255)
                              gctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                     >  2:2117116‑M1/18‑149 (MQ=255)
                              gctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTCACCAGTGAGAc                                                                                                                     <  1:1850470‑M1/132‑1 (MQ=255)
                               ctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                    <  1:3350109‑M1/133‑1 (MQ=255)
                                taacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                   >  1:1570278‑M1/16‑149 (MQ=255)
                                   gTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCa                                                                                                                <  1:3263992/149‑1 (MQ=255)
                                   cttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCa                                                                                                                >  1:5016144‑M1/13‑149 (MQ=255)
                                     tacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                                                                                              <  2:1006246‑M1/139‑1 (MQ=255)
                                      acattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTAAACAGGGAGACGGGCcaca                                                                                                             <  1:2665985‑M1/140‑1 (MQ=255)
                                        attaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTAGTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc                                                                                                           <  2:4923340‑M1/142‑1 (MQ=255)
                                         ttaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATGGGGCGCCAGGGTGGTGTTTCTTTTAACCAGTGAGACGGGCCACAgca                                                                                                          <  1:4813456‑M1/143‑2 (MQ=255)
                                           aattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGa                                                                                                        <  2:726464‑M1/145‑1 (MQ=255)
                                             ttGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTGTTTTCACCAGTGAGACGGGCCACAGCTGAtt                                                                                                      <  1:126886‑M1/147‑1 (MQ=255)
                                             ttGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                                                                                      >  2:2485001‑M1/3‑149 (MQ=255)
                                             ttGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGCGCCAGGGTGGTGTTTATTTTAAACAGTGAGAAGGGCAACAGCTGAtt                                                                                                      <  1:528647‑M1/147‑1 (MQ=255)
                                               gcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGc                                                                                                    <  2:1478964/149‑1 (MQ=255)
                                                 gTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGccc                                                                                                  <  2:192532/149‑1 (MQ=255)
                                                 gTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTAACCAGTGAGAAGGGCAACAGCAGATTTccc                                                                                                  >  1:622718/1‑149 (MQ=255)
                                                      gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAc                                                                                             <  1:2563890/149‑1 (MQ=255)
                                                        tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCg                                                                                           >  2:4094981/1‑149 (MQ=255)
                                                            tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                                                                                       <  1:4117160/149‑1 (MQ=255)
                                                             gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTgg                                                                                      >  1:3728174/1‑149 (MQ=255)
                                                               ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGcc                                                                                    >  1:4719069/1‑149 (MQ=255)
                                                                 gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCt                                                                                  <  1:4020240/149‑1 (MQ=255)
                                                                     tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgaga                                                                              <  1:1000192/149‑1 (MQ=255)
                                                                      ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagag                                                                             <  2:1739058/149‑1 (MQ=255)
                                                                       cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTGTTTTAACCAGTGAGACGGGCAAAAGCTGATTTCCCATCACCGCCAGGCCCGgagaga                                                                            >  1:1089769/1‑149 (MQ=255)
                                                                       cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagaga                                                                            >  2:3722468/1‑149 (MQ=255)
                                                                       cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTGTTTTTTTTAACCAGTGAGAAGGGGAACACCGGATTTCCCATAACCCCCAGGCCCCgagaga                                                                            >  1:3918500/1‑149 (MQ=255)
                                                                        aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGagagag                                                                           <  1:4197980/149‑1 (MQ=255)
                                                                         gTCGGGAAACCTTTCGTGCCAGCTGCATTAATTAATCGGCCAACGAGCGGGGAGAGGCGGTTTGCGTATTGGTCTCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGAAAAAGCTGAATGACCTTCAACGCCTGGCCCTGAGAGAGt                                                                          >  2:2575798/1‑149 (MQ=255)
                                                                             ggACACCTGTCGTGCCAGCTGCATTAATGAATCGGCCCACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                      <  2:3262640/149‑1 (MQ=255)
                                                                              gAAACCTGTCGTGCCAGCTGCATTAATGAATCGTCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGCTCCAGGGTGGTTTTTATTTTAACCACTGATACGGGCAACACCAGATTTACCTTCACCGCCAGGCCCCGACAGAGTTtcac                                                                     >  1:4558726/1‑148 (MQ=255)
                                                                              gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGCGCCAGGGTGGTTTTTATTTTAACCAGTGAGACGGGCAACACCTGTTTTCCCAACACCCCCAGGCCCGCAGAGAGTTggag                                                                     >  1:899339/1‑149 (MQ=255)
                                                                               aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagc                                                                    >  1:2475289/1‑149 (MQ=255)
                                                                                aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTAACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCAGGCCCGGAGAGTGTTgcagcc                                                                   >  2:900106/1‑148 (MQ=255)
                                                                                 aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                  <  2:634160/149‑1 (MQ=255)
                                                                                   cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCATAACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                >  1:4644828/1‑149 (MQ=255)
                                                                                      tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                             <  2:2399906/149‑1 (MQ=255)
                                                                                         tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                          <  2:4566104/149‑1 (MQ=255)
                                                                                               cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTgg                                                    >  2:1100969/1‑149 (MQ=255)
                                                                                               cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTgg                                                    <  2:335416/149‑1 (MQ=255)
                                                                                               cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTAACCAGTGAGAAGGGCAAAACCAGATTTCCCGTCACCCCCAGGCCCGGAGAGAGTTTCAGCAAGCGGTCCAcccagg                                                    >  1:1379307/1‑144 (MQ=255)
                                                                                                tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTCTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACCGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCGGAGAGAGTTGCAGCAAGCGGTCCACGCTGGt                                                   >  2:4837596/1‑149 (MQ=255)
                                                                                                tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGt                                                   >  2:1561488/1‑149 (MQ=255)
                                                                                                 gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTTTTTTAACCACTGATAAGGGCAACAGCGGATTTCCCATAACCGCCAGGCCCAGAGGGAGTTGCAGCAACCCGGCCAcccatgtt                                                  >  1:3540019/1‑142 (MQ=255)
                                                                                                  cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTAACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGttt                                                 >  2:2764342/1‑149 (MQ=255)
                                                                                                   aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTCACCAGTGAGACGGGCAAAAGCTTATTGCCCATCACCGCCTGGCCCAGAGAGAGTTGCAGCAAGCGGTCCACGCAGGTTTg                                                >  2:1834187/1‑149 (MQ=255)
                                                                                                      aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGccc                                             >  2:1385225/1‑149 (MQ=255)
                                                                                                       aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcccc                                            <  2:3173198/149‑1 (MQ=255)
                                                                                                       aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCGGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcccc                                            >  1:3445168/1‑149 (MQ=255)
                                                                                                            tCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCGTCACCGCCTGGCCCGGGGAGAGGTGCAGCAAGCGGTCCACGCTGGGTTGCCCccgccg                                       >  1:2387149/1‑147 (MQ=255)
                                                                                                              ggCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                     <  1:1762352/149‑1 (MQ=255)
                                                                                                                ccAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                   <  2:535892/149‑1 (MQ=255)
                                                                                                                   aCGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaa                                <  2:2024207/149‑1 (MQ=255)
                                                                                                                    cgcgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaaa                               >  2:4381376/1‑148 (MQ=255)
                                                                                                                     gcgcgGGGAGAGGCGGTGTGCGGATTGGGCGCCAGGGTGGTGTTTCTTTTCACCAGTGAGACGGGCAACAGCGGATTGCCCGTCACCGCCTGGCCCGGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATc                              >  2:818931/1‑149 (MQ=255)
                                                                                                                      cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATcc                             <  2:3080007/149‑1 (MQ=255)
                                                                                                                      cgcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATcc                             <  1:2113544/149‑1 (MQ=255)
                                                                                                                       gcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                            <  1:3078809/149‑1 (MQ=255)
                                                                                                                       gcgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                            <  1:4364080/149‑1 (MQ=255)
                                                                                                                          gggAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTCACCACTGATACGGGCAACAGCTGATTTCCCTTCACCGCCTGGCCCTGAGAGAGTTTCAGCAAGCGGTCCACCCAGGTTTTCCCCAACAGGCGAAAATCCAGtt                         >  2:3899825/1‑149 (MQ=255)
                                                                                                                             agagGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa                      >  2:4657771/1‑149 (MQ=255)
                                                                                                                             agagGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTATTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGACCACGCTTGTTTGCCCCAGCAGGCGAAAATACTGtttta                      >  2:1292083/1‑147 (MQ=255)
                                                                                                                              gagGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGATACGGGCAACAGCTGATTGCCCTTCACCGCCAGGCCCGGAGAGAGTTGCAGCAAGCGGTCCACCCTGGTTTTCCCCACCAGGCGAAAATCCTGTTTGAt                     >  1:4284350/1‑149 (MQ=255)
                                                                                                                                ggCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                   >  2:4744926/1‑149 (MQ=255)
                                                                                                                                 gCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                  <  1:3167776/149‑1 (MQ=255)
                                                                                                                                 gCGGTGTGCGTATTGGGCGCCAGGGTGGTGTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                  <  2:5112301/149‑1 (MQ=255)
                                                                                                                                  cGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                 <  1:148252/149‑1 (MQ=255)
                                                                                                                                   ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                <  1:5120493/149‑1 (MQ=255)
                                                                                                                                    gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt               <  1:4908368/149‑1 (MQ=255)
                                                                                                                                    gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACACCAGATTGCCCTTCACCCCCTGGCCCAGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTTCCCCACCAGGCGAAAATCCAGTTTGATggtggt               >  1:1937618/1‑149 (MQ=255)
                                                                                                                                    gTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt               <  1:2194202/149‑1 (MQ=255)
                                                                                                                                     tttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt              <  2:4280150/149‑1 (MQ=255)
                                                                                                                                     tttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt              <  2:3545981/149‑1 (MQ=255)
                                                                                                                                           tATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc        >  2:2971665/1‑149 (MQ=255)
                                                                                                                                            aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCAGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg       >  2:2208766/1‑149 (MQ=255)
                                                                                                                                              tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg     <  1:2881989/149‑1 (MQ=255)
                                                                                                                                               gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa    <  1:2977554/149‑1 (MQ=255)
                                                                                                                                               gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa    <  1:2020285/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATA  >  CP009273/360605‑360898

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: ‑