Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,512,106 G→T 45.5% E68* (GAA→TAA)  ydcY → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,512,1060GT45.5% 166.4 / 229.7 201E68* (GAA→TAA) ydcYuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (53/56);  new base T (50/42);  total (103/98)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

TCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACTGCGGAAGATATCCAGCATC  >  CP009273/1511960‑1512253
                                                                                                                                                  |                                                                                                                                                   
 cATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTaa                                                                                                                                                   >  1:1761872/1‑148 (MQ=255)
 cATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGaa                                                                                                                                                   >  1:3047068/1‑148 (MQ=255)
    tGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAg                                                                                                                                                >  1:1027487/1‑148 (MQ=255)
       ccATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTa                                                                                                                                             <  2:1400901/148‑1 (MQ=255)
       ccATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTa                                                                                                                                             <  2:1453643/148‑1 (MQ=255)
           atGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCa                                                                                                                                         <  2:1034907/148‑1 (MQ=255)
           atGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCa                                                                                                                                         >  2:3457990/1‑148 (MQ=255)
            tGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCaa                                                                                                                                        <  1:3457990/148‑1 (MQ=255)
             gAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCaaa                                                                                                                                       <  1:3169333/148‑1 (MQ=255)
             gAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCaaa                                                                                                                                       >  1:2687045/1‑148 (MQ=255)
              aaCGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAg                                                                                                                                      >  1:1376529/1‑148 (MQ=255)
               aCGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAgg                                                                                                                                     >  1:948143/1‑148 (MQ=255)
               aCGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAgg                                                                                                                                     >  1:934127/1‑148 (MQ=255)
                cGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAggc                                                                                                                                    <  1:1986645/148‑1 (MQ=255)
                    ttattaATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCa                                                                                                                                >  2:3482250/1‑148 (MQ=255)
                    ttattaATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCa                                                                                                                                <  2:1496613/148‑1 (MQ=255)
                     tattaATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCAc                                                                                                                               <  1:333245/148‑1 (MQ=255)
                     tattaATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCAc                                                                                                                               >  1:2044476/1‑148 (MQ=255)
                          aTCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACTAAAGGCGGCACCATcc                                                                                                                          >  2:3015248/1‑148 (MQ=255)
                            cGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCctc                                                                                                                        >  1:1222592/1‑148 (MQ=255)
                            cGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCctc                                                                                                                        >  1:3002447/1‑148 (MQ=255)
                             gCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCctct                                                                                                                       <  1:2942428/148‑1 (MQ=255)
                               ccTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGa                                                                                                                     <  2:286461/148‑1 (MQ=255)
                                cTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGAt                                                                                                                    <  2:280725/148‑1 (MQ=255)
                                  gAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTa                                                                                                                  >  1:2737673/1‑148 (MQ=255)
                                   aGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCa                                                                                                                             <  1:824523/136‑1 (MQ=255)
                                   aGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCa                                                                                                                             >  2:824523/1‑136 (MQ=255)
                                    gCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTaaa                                                                                                                <  2:1286396/148‑1 (MQ=255)
                                       aTGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATg                                                                                                             <  2:1902848/148‑1 (MQ=255)
                                       aTGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATg                                                                                                             >  1:1524999/1‑148 (MQ=255)
                                             cgGAGTTAAGTCGGGAAGATCGCTACACCCCGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTc                                                                                                       <  2:888920/148‑1 (MQ=255)
                                             cgGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTc                                                                                                       >  1:2638721/1‑148 (MQ=255)
                                               gAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTg                                                                                                     >  2:1735325/1‑148 (MQ=255)
                                                     aGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTacc                                                                                               >  2:762770/1‑148 (MQ=255)
                                                       tCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTaccaa                                                                                             <  1:1229258/148‑1 (MQ=255)
                                                          ggAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTaccaacca                                                                                          >  1:756201/1‑148 (MQ=255)
                                                          ggAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTaccaacca                                                                                          >  1:3488605/1‑148 (MQ=255)
                                                          ggAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTaccaacca                                                                                          <  1:2322335/148‑1 (MQ=255)
                                                           gAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAg                                                                                         >  2:37388/1‑148 (MQ=255)
                                                               aTCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGgcgc                                                                                     >  2:2898976/1‑148 (MQ=255)
                                                               aTCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGgcgc                                                                                     <  2:1484543/148‑1 (MQ=255)
                                                                  gCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCaa                                                                                  <  1:1196703/148‑1 (MQ=255)
                                                                  gCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCaa                                                                                  <  1:384558/148‑1 (MQ=255)
                                                                    tACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAAcc                                                                                >  2:3465282/1‑148 (MQ=255)
                                                                     acacCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACca                                                                               >  1:623184/1‑148 (MQ=255)
                                                                     acacCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACca                                                                               <  2:703455/148‑1 (MQ=255)
                                                                       acCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACcaca                                                                             <  2:605414/148‑1 (MQ=255)
                                                                        cccAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACaa                                                                            <  1:1570709/148‑1 (MQ=255)
                                                                          cagcagCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGt                                                                          <  2:596425/148‑1 (MQ=255)
                                                                          cagcagCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGt                                                                          <  2:2252709/148‑1 (MQ=255)
                                                                             cagcaACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTaaa                                                                       <  1:3280305/148‑1 (MQ=255)
                                                                              agcaACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTaaaa                                                                      <  2:900502/148‑1 (MQ=255)
                                                                              agcaACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTaaaa                                                                      <  2:2733073/148‑1 (MQ=255)
                                                                                caACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAAtt                                                                    >  1:3065505/1‑148 (MQ=255)
                                                                                caACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAAtt                                                                    >  2:2241534/1‑148 (MQ=255)
                                                                                   cGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTgcg                                                                 <  1:3327523/148‑1 (MQ=255)
                                                                                   cGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTgcg                                                                 <  1:876938/148‑1 (MQ=255)
                                                                                      cTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCca                                                              <  2:1769902/148‑1 (MQ=255)
                                                                                      cTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCACCCACAAGTAAAATTGCGCca                                                              >  1:2839387/1‑148 (MQ=255)
                                                                                      cTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCca                                                              >  1:566974/1‑148 (MQ=255)
                                                                                      cTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCca                                                              <  2:2044476/148‑1 (MQ=255)
                                                                                           cacaGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAACGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAAt                                                         <  2:3145257/148‑1 (MQ=255)
                                                                                           cacaGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAATTAAAATTGCGCCACACAt                                                         <  2:1372319/148‑1 (MQ=255)
                                                                                           cacaGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAAt                                                         <  1:2867275/148‑1 (MQ=255)
                                                                                            acaGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAAtt                                                        <  1:2966852/148‑1 (MQ=255)
                                                                                               gCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTggg                                                     <  1:762770/148‑1 (MQ=255)
                                                                                                cGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGc                                                    >  1:915713/1‑148 (MQ=255)
                                                                                                 gATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGcc                                                   >  2:2554667/1‑148 (MQ=255)
                                                                                                      cccATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGaaca                                              >  1:1142974/1‑148 (MQ=255)
                                                                                                         aTCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACg                                           <  2:2391320/148‑1 (MQ=255)
                                                                                                         aTCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACg                                           <  2:413143/148‑1 (MQ=255)
                                                                                                          tCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACNACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGc                                          <  1:2090573/148‑1 (MQ=255)
                                                                                                          tCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGc                                          >  2:1565780/1‑148 (MQ=255)
                                                                                                           cACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGcc                                         >  2:2591721/1‑148 (MQ=255)
                                                                                                            aCGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGccc                                        >  2:1670135/1‑148 (MQ=255)
                                                                                                              ggTCGAAAGAATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGt                                      <  1:3482424/148‑1 (MQ=255)
                                                                                                               gTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTa                                     >  2:626591/1‑148 (MQ=255)
                                                                                                                tCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAg                                    <  2:2687045/148‑1 (MQ=255)
                                                                                                                 cGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGt                                   >  1:1664052/1‑148 (MQ=255)
                                                                                                                  gaaaGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTg                                  <  2:1219343/146‑1 (MQ=255)
                                                                                                                   cAAGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTgg                                 <  1:522436/148‑1 (MQ=255)
                                                                                                                   cAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTgg                                 <  1:2589833/148‑1 (MQ=255)
                                                                                                                    aaGCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCGGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGt                                <  1:2207021/148‑1 (MQ=255)
                                                                                                                      gCATAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTAc                              >  2:3111978/1‑148 (MQ=255)
                                                                                                                      gCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTAc                              >  1:3130793/1‑148 (MQ=255)
                                                                                                                       cATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACg                             >  2:2649657/1‑148 (MQ=255)
                                                                                                                         tAAAGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGcc                           <  2:3326567/148‑1 (MQ=255)
                                                                                                                          aaaGAAGATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCa                          <  2:1578525/148‑1 (MQ=255)
                                                                                                                           aagaagATATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCAt                         >  2:324819/1‑148 (MQ=255)
                                                                                                                            agaagaTATGGAAGCGCGCCACTAACAGTTAACCCAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATg                        <  2:3030570/148‑1 (MQ=255)
                                                                                                                             gaagaTATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGa                       <  2:1223434/148‑1 (MQ=255)
                                                                                                                               agaTATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACt                     >  1:1420744/1‑148 (MQ=255)
                                                                                                                                gaTATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACTg                    >  2:3406581/1‑148 (MQ=255)
                                                                                                                                gaTATGGAAGCGCGCCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACTg                    >  2:893698/1‑148 (MQ=255)
                                                                                                                                 atatGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACTGc                   >  1:3523976/1‑148 (MQ=255)
                                                                                                                                        aGCGCGCAACTAACAGTTAACCAAAGGCGTCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCACCCACAAGTAAAATGGCCCCACAAATTGGGCCGAAAAACGCCCGTAGTGGTACGCCATGACTGCGGAAGat            <  2:3002447/148‑1 (MQ=255)
                                                                                                                                         gcgcgcCACGAACAGTTAACAAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACTGCGGAAGata           <  1:1604584/148‑1 (MQ=255)
                                                                                                                                           gcgcCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACTGCGGAAGATATc         <  2:1362608/148‑1 (MQ=255)
                                                                                                                                             gcCACTAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACTGCGGAAGATATCCa       >  2:1720793/1‑148 (MQ=255)
                                                                                                                                                  gAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACTGCGGAAGATATCCAGCATc  >  1:2397859/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAGAAGATATGGAAGCGCGCCACGAACAGTTAACCAAAGGCGGCACCATCCTCTGATTAAAATGAGCGTCTGGCTACCAACCAGGCGCCAACCACAAGTAAAATTGCGCCACAAATTGGGCCGAACAACGCCCGTAGTGGTACGCCATGACTGCGGAAGATATCCAGCATC  >  CP009273/1511960‑1512253

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