Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,204,284 C→A 89.9% D98Y (GAT→TAT)  tfaE ← e14 prophage; putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,204,2840CA89.9% 321.8 / 33.1 119D98Y (GAT→TAT) tfaEe14 prophage; putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base A (51/56);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 119 total aligned reads displayed.

TTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAC  >  CP009273/1204137‑1204430
                                                                                                                                                   |                                                                                                                                                  
ttGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAAc                                                                                                                                                  <  2:2920600/150‑1 (MQ=255)
 tGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAAcc                                                                                                                                                 >  1:4039909/1‑150 (MQ=255)
   cagcaTCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCg                                                                                                                                               >  2:2977081/1‑150 (MQ=255)
       atcaatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGc                                                                                                                                           <  1:286244/150‑1 (MQ=255)
       atcaatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGGATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGc                                                                                                                                           >  2:3965498/1‑150 (MQ=255)
         caatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCtt                                                                                                                                         >  1:2509508/1‑150 (MQ=255)
          aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                        <  1:2362467/150‑1 (MQ=255)
          aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                        <  2:1079767/150‑1 (MQ=255)
          aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                        <  2:2880409/150‑1 (MQ=255)
          aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                        <  2:3493959/150‑1 (MQ=255)
           atcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGGTGTATTTTCAGGGTAATAACCCGGAGCTTTg                                                                                                                                       >  2:2015961/1‑150 (MQ=255)
              aGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAtt                                                                                                                                    >  1:734425/1‑150 (MQ=255)
              aGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAtt                                                                                                                                    <  2:3456196/150‑1 (MQ=255)
              aGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAtt                                                                                                                                    <  2:3456000/150‑1 (MQ=255)
              aGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAtt                                                                                                                                    >  1:3489602/1‑150 (MQ=255)
                tGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTc                                                                                                                                  >  2:1370190/1‑150 (MQ=255)
                 gACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCt                                                                                                                                 >  1:3114156/1‑150 (MQ=255)
                  aCTGGCGCTGTGCTTGTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAATGGGGCTATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCtt                                                                                                                                <  1:3230106/150‑1 (MQ=255)
                  aCTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATACAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCTGAGCTTTGATTTCtt                                                                                                                                >  2:2803618/1‑150 (MQ=255)
                  aCTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCtt                                                                                                                                <  1:2856708/150‑1 (MQ=255)
                   cTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCttt                                                                                                                               <  2:2813398/150‑1 (MQ=255)
                       cgcTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGAt                                                                                                                           <  1:1934547/150‑1 (MQ=255)
                              cTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAg                                                                                                                    <  2:2963001/150‑1 (MQ=255)
                              cTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAg                                                                                                                    >  1:456584/1‑150 (MQ=255)
                              cTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAg                                                                                                                    >  1:1079288/1‑150 (MQ=255)
                              cTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAg                                                                                                                    >  1:1450434/1‑150 (MQ=255)
                                 cTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGttt                                                                                                                 >  1:173664/1‑150 (MQ=255)
                                   gCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTg                                                                                                               <  2:2155422/150‑1 (MQ=255)
                                   gCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTg                                                                                                               <  2:2155485/150‑1 (MQ=255)
                                    ccGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTgg                                                                                                              >  1:2833025/1‑150 (MQ=255)
                                    acgagTCTACTGCGGCGAAAAGCTGTGCTTCAGAATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAAGGGCGATAGTGGTATTATTTTCAGGGTAATAACCCGTAGCATTGATTACATTTGATTCTCAAGTTTTgg                                                                                                              >  2:1490090/5‑150 (MQ=25)
                                      gcgTCTACTGCGGCGCTATGCTGTGATTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGTGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTg                                                                                                            <  1:430977/150‑1 (MQ=255)
                                           tACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCtata                                                                                                       <  2:2577898/150‑1 (MQ=255)
                                           tACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCtata                                                                                                       >  1:490182/1‑150 (MQ=255)
                                           tACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCtata                                                                                                       >  1:3939708/1‑150 (MQ=255)
                                              tGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGAc                                                                                                    <  2:887544/150‑1 (MQ=255)
                                              tGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGAc                                                                                                    >  1:2626349/1‑150 (MQ=255)
                                                    gcTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTc                                                                                              >  1:1492800/1‑150 (MQ=255)
                                                               cTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGt                                                                                   <  2:1461668/150‑1 (MQ=255)
                                                               cTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGt                                                                                   <  2:1461955/150‑1 (MQ=255)
                                                                ttCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTc                                                                                  <  2:1740868/150‑1 (MQ=255)
                                                                 tCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCt                                                                                 <  1:4156145/150‑1 (MQ=255)
                                                                 tCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCt                                                                                 <  1:4156215/150‑1 (MQ=255)
                                                                   aGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTgg                                                                               <  1:2682643/150‑1 (MQ=255)
                                                                     tATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAGTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGca                                                                             >  2:1223310/1‑150 (MQ=255)
                                                                      aTCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGcac                                                                            <  2:1725623/150‑1 (MQ=255)
                                                                       tCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGcaca                                                                           >  1:1587942/1‑150 (MQ=255)
                                                                        cGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACat                                                                          >  1:1805297/1‑150 (MQ=255)
                                                                        cGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACat                                                                          >  1:1804925/1‑150 (MQ=255)
                                                                        cGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTGGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGGTTCACCGCGATGGTCTGGcacct                                                                          >  2:2840837/1‑148 (MQ=255)
                                                                         ggTCACCCATTTCTCACAATCCCAATTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATTGACCGTTTCACAGCGATGGTCTGGCACata                                                                         >  2:3551468/1‑150 (MQ=255)
                                                                          gTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGTTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACatat                                                                        <  2:2009130/150‑1 (MQ=255)
                                                                           tCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCGCATAtt                                                                       >  1:2429618/1‑150 (MQ=255)
                                                                                     ctcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGtt                                                             <  2:2245202/150‑1 (MQ=255)
                                                                                     ctcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGtt                                                             <  2:2245194/150‑1 (MQ=255)
                                                                                     ctcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGtt                                                             <  2:2245196/150‑1 (MQ=255)
                                                                                      tcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTa                                                            >  2:1114967/1‑150 (MQ=255)
                                                                                          cATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAAt                                                        <  1:2473495/150‑1 (MQ=255)
                                                                                           aTCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATc                                                       <  1:3959671/150‑1 (MQ=255)
                                                                                           aTCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATc                                                       <  1:3959816/150‑1 (MQ=255)
                                                                                            tCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCt                                                      <  2:3467903/150‑1 (MQ=255)
                                                                                            tCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCt                                                      <  2:3467892/150‑1 (MQ=255)
                                                                                              ccATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGc                                                    <  1:1302748/150‑1 (MQ=255)
                                                                                               cATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCa                                                   >  2:2127534/1‑150 (MQ=255)
                                                                                               cATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCa                                                   >  2:2219134/1‑150 (MQ=255)
                                                                                                 tttATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGa                                                 <  1:3038726/150‑1 (MQ=255)
                                                                                                      cGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGc                                            >  1:163708/1‑150 (MQ=255)
                                                                                                           ggAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGAtt                                       >  1:2054212/1‑150 (MQ=255)
                                                                                                               aTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCAt                                   <  1:213425/150‑1 (MQ=255)
                                                                                                               aTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCAt                                   <  1:213452/150‑1 (MQ=255)
                                                                                                                 aaaGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATaa                                 >  2:809421/1‑150 (MQ=255)
                                                                                                                 aaaGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATaa                                 >  2:1849163/1‑150 (MQ=255)
                                                                                                                   aGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAAcc                               <  2:826915/150‑1 (MQ=255)
                                                                                                                     gggCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAg                             >  2:3150326/1‑150 (MQ=255)
                                                                                                                     gggCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAg                             >  2:3150314/1‑150 (MQ=255)
                                                                                                                      ggCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGc                            >  1:710680/1‑150 (MQ=255)
                                                                                                                       gCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGcc                           >  1:3270379/1‑150 (MQ=255)
                                                                                                                         taTAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                         <  2:1667503/149‑1 (MQ=255)
                                                                                                                         gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                         >  1:2155424/1‑150 (MQ=255)
                                                                                                                         gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                         >  2:1790463/1‑150 (MQ=255)
                                                                                                                         gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                         >  2:1790394/1‑150 (MQ=255)
                                                                                                                           tAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAt                       <  1:2268497/150‑1 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     <  1:3726442/150‑1 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     >  2:4003375/1‑150 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     <  1:3726463/150‑1 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     >  2:3551823/1‑150 (MQ=255)
                                                                                                                             gTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgt                     <  1:3726459/150‑1 (MQ=255)
                                                                                                                                gTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAAACAGCCTTATGTGTa                  >  1:2148279/1‑150 (MQ=255)
                                                                                                                                  tGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTAcc                <  2:382696/150‑1 (MQ=255)
                                                                                                                                   gTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCAGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCa               >  2:2958049/1‑150 (MQ=255)
                                                                                                                                      ttttCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAggg            <  1:3381132/150‑1 (MQ=255)
                                                                                                                                         tCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCa         <  1:3086838/150‑1 (MQ=255)
                                                                                                                                          cAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCAt        >  2:2989417/1‑150 (MQ=255)
                                                                                                                                           aGGGTAATCACCCGGAGCTGTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTACCATGAGTTAAAATCTGCAGAATGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATc       >  2:1476022/1‑150 (MQ=255)
                                                                                                                                           aGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATc       >  2:1468094/1‑150 (MQ=255)
                                                                                                                                            gggTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCt      <  2:4063094/150‑1 (MQ=255)
                                                                                                                                            gggTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCt      <  2:557746/150‑1 (MQ=255)
                                                                                                                                             ggTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTa     <  1:2990435/150‑1 (MQ=255)
                                                                                                                                             ggTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTa     <  1:2990431/150‑1 (MQ=255)
                                                                                                                                                aaTAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAc  <  1:918957/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAC  >  CP009273/1204137‑1204430

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