Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,340,038 C→A 6.5% I34I (ATC→ATA ptsN → sugar‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,340,0380CA 6.5% 707.2 / 16.5 190I34I (ATC→ATAptsNsugar‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (85/92);  new base A (7/6);  total (92/98)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GAGCCTGTTCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCC  >  CP009273/3339893‑3340185
                                                                                                                                                 |                                                                                                                                                   
gAGCCTGTTCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAatcat                                                                                                                                                   >  1:1242069/1‑148 (MQ=255)
     tGTTCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCg                                                                                                                                              <  1:2089230/148‑1 (MQ=255)
        tCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGc                                                                                                                                           >  1:2053355/1‑148 (MQ=255)
           aCTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTgg                                                                                                                                        <  1:2345659/148‑1 (MQ=255)
           aCTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTgg                                                                                                                                        <  1:11984/148‑1 (MQ=255)
              gTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTggcgg                                                                                                                                     >  2:1170444/1‑148 (MQ=255)
                         cAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACtt                                                                                                                          >  2:1551185/1‑148 (MQ=255)
                           ggTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAg                                                                                                                        >  1:608865/1‑148 (MQ=255)
                                    gttaaaTTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTTCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTgccgcc                                                                                                               <  1:1782873/145‑1 (MQ=255)
                                    gTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTgccgcc                                                                                                               >  2:2041561/1‑148 (MQ=255)
                                     tGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCt                                                                                                              <  2:2436850/148‑1 (MQ=255)
                                      gAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTc                                                                                                             <  2:411837/148‑1 (MQ=255)
                                      gAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTc                                                                                                             <  2:1724346/148‑1 (MQ=255)
                                       aaaTTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCa                                                                                                            <  2:2644430/148‑1 (MQ=255)
                                         aTTATGTCTAATAATGATACAACTCTAAAGTTTATCAGTTTTCTTAACAGGGAATGTTCGCGAAGCAGCGTCTAGTGTCTTATGAATAAACGATGCCTGGAAATCATCAGCGAGCTGGCGGGGATTCATCTTTGCCTGACGCCTCAgg                                                                                                          <  1:1391520/148‑1 (MQ=255)
                                         aTTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAgg                                                                                                          >  1:688785/1‑148 (MQ=255)
                                          ttATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAggt                                                                                                         <  1:1243960/148‑1 (MQ=255)
                                          ttATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAggt                                                                                                         >  1:254182/1‑148 (MQ=255)
                                          ttATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAggt                                                                                                         >  1:253889/1‑148 (MQ=255)
                                            aTGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAggtgg                                                                                                       >  1:1498336/1‑148 (MQ=255)
                                             tGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAggtggt                                                                                                      <  1:1144261/148‑1 (MQ=255)
                                              gACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGtt                                                                                                     >  1:2073053/1‑148 (MQ=255)
                                              gACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGtt                                                                                                     >  1:2129887/1‑148 (MQ=255)
                                              gACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGtt                                                                                                     >  2:35223/1‑148 (MQ=255)
                                               aCAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGttt                                                                                                    >  2:1747451/1‑148 (MQ=255)
                                                cAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGtttt                                                                                                   <  1:2950588/148‑1 (MQ=255)
                                                     aatGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAg                                                                                              <  2:2504773/148‑1 (MQ=255)
                                                        gATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTa                                                                                           <  1:449173/148‑1 (MQ=255)
                                                         aTACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTAt                                                                                          <  1:547873/148‑1 (MQ=255)
                                                          tACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATc                                                                                         <  1:238332/148‑1 (MQ=255)
                                                          tACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATc                                                                                         <  2:2731660/148‑1 (MQ=255)
                                                          tACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATc                                                                                         <  2:910523/148‑1 (MQ=255)
                                                            cAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCt                                                                                       <  2:2863370/148‑1 (MQ=255)
                                                             aaCTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTg                                                                                      <  2:271667/148‑1 (MQ=255)
                                                             aaCTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTg                                                                                      <  2:271705/148‑1 (MQ=255)
                                                             aaCTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTg                                                                                      <  2:1044866/148‑1 (MQ=255)
                                                               ctctACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTTTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGAc                                                                                    >  1:100656/1‑148 (MQ=255)
                                                                 ctACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGAcgc                                                                                  >  1:1036199/1‑148 (MQ=255)
                                                                    cAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCggg                                                                               >  2:1118151/1‑148 (MQ=255)
                                                                     aGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGa                                                                              >  2:768900/1‑148 (MQ=255)
                                                                      gCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaa                                                                             >  1:2780130/1‑148 (MQ=255)
                                                                       cTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaa                                                                            >  1:2999025/1‑148 (MQ=255)
                                                                       cTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaa                                                                            >  1:1352518/1‑148 (MQ=255)
                                                                         tAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaaaa                                                                          >  1:510997/1‑148 (MQ=255)
                                                                         tAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaaaa                                                                          >  1:511247/1‑148 (MQ=255)
                                                                         tAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaaaa                                                                          >  2:1041046/1‑148 (MQ=255)
                                                                          aGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaaaaa                                                                         <  2:1880548/148‑1 (MQ=255)
                                                                               tgtTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAACTGGGc                                                                    >  1:128204/1‑148 (MQ=255)
                                                                               tgtTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGc                                                                    >  2:1406886/1‑148 (MQ=255)
                                                                                gtTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCa                                                                   <  2:3058714/148‑1 (MQ=255)
                                                                                  tCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGt                                                                 >  2:2110424/1‑148 (MQ=255)
                                                                                   cTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTa                                                                <  1:789282/148‑1 (MQ=255)
                                                                                     tAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTAcc                                                              >  1:3241426/1‑148 (MQ=255)
                                                                                      aaCAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCg                                                             <  2:3308024/148‑1 (MQ=255)
                                                                                       aCAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCgg                                                            >  2:338728/1‑148 (MQ=255)
                                                                                       aCAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCgg                                                            <  1:1084459/148‑1 (MQ=255)
                                                                                        cAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGt                                                           <  2:670003/148‑1 (MQ=255)
                                                                                         aGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTa                                                          <  2:2020291/148‑1 (MQ=255)
                                                                                          gggAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTAt                                                         <  1:1598382/148‑1 (MQ=255)
                                                                                          gggAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTAt                                                         <  2:2205658/148‑1 (MQ=255)
                                                                                            gAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCg                                                       >  2:1240473/1‑148 (MQ=255)
                                                                                                gTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCaa                                                   <  1:1787766/148‑1 (MQ=255)
                                                                                                gTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCaa                                                   <  1:116027/148‑1 (MQ=255)
                                                                                                 gaCGCGCAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAAt                                                  <  2:1770100/147‑1 (MQ=255)
                                                                                                  aCGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATg                                                 <  1:1529335/148‑1 (MQ=255)
                                                                                                  aCGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATg                                                 <  2:2450511/148‑1 (MQ=255)
                                                                                                  aCGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATg                                                 <  2:2331697/148‑1 (MQ=255)
                                                                                                   cgcgAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATgg                                                <  2:2995016/148‑1 (MQ=255)
                                                                                                       aaGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTAtt                                            >  1:3072413/1‑148 (MQ=255)
                                                                                                        aGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTg                                           <  1:2924949/148‑1 (MQ=255)
                                                                                                        aGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTg                                           <  1:955999/148‑1 (MQ=255)
                                                                                                             cgTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCAtt                                      >  2:772064/1‑148 (MQ=255)
                                                                                                                 cACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCACTGGTATTGCCATTCCGc                                  >  1:2697723/1‑148 (MQ=255)
                                                                                                                      tCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGc                             >  2:2011932/1‑148 (MQ=255)
                                                                                                                        agagCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCaa                           <  2:2053355/148‑1 (MQ=255)
                                                                                                                        agagCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCaa                           <  2:445909/148‑1 (MQ=255)
                                                                                                                         gagCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCaaa                          <  1:768900/148‑1 (MQ=255)
                                                                                                                         gagCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCaaa                          >  2:2499996/1‑148 (MQ=255)
                                                                                                                         gagCAAAAAACGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCaaa                          <  2:2256593/148‑1 (MQ=255)
                                                                                                                          agCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAAc                         <  1:264089/148‑1 (MQ=255)
                                                                                                                              aaaaaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGga                     <  2:2687764/148‑1 (MQ=255)
                                                                                                                              aaaaaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGga                     <  2:1298497/148‑1 (MQ=255)
                                                                                                                               aaaaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaa                    <  1:1406886/148‑1 (MQ=255)
                                                                                                                               aaaaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaa                    >  1:251480/1‑148 (MQ=255)
                                                                                                                                aaaCGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaag                   >  2:2836767/1‑148 (MQ=255)
                                                                                                                                 aaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaaga                  >  2:3202110/1‑148 (MQ=255)
                                                                                                                                 aaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaaga                  >  1:2187609/1‑148 (MQ=255)
                                                                                                                                    gcgcCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAagaaga               <  2:2413196/148‑1 (MQ=255)
                                                                                                                                       cccTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATAc            >  2:3049423/1‑148 (MQ=255)
                                                                                                                                       cccTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATAc            <  1:1118151/148‑1 (MQ=255)
                                                                                                                                         cTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATActc          >  2:326147/1‑148 (MQ=255)
                                                                                                                                          tGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATActct         <  2:100656/148‑1 (MQ=255)
                                                                                                                                          tGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATActct         >  1:487732/1‑148 (MQ=255)
                                                                                                                                           ggAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCGTTCCGCATGGCAAACTGGAAGAAGATACTCTg        <  2:1242088/148‑1 (MQ=255)
                                                                                                                                            gAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTgc       >  2:1649156/1‑148 (MQ=255)
                                                                                                                                            gAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTgc       >  2:1649090/1‑148 (MQ=255)
                                                                                                                                            gAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTgc       >  1:1975446/1‑148 (MQ=255)
                                                                                                                                             aaaTCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTgcg      >  1:207787/1‑148 (MQ=255)
                                                                                                                                                 catccGCGAGCTGGCGGCGACACAACTGAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCAGTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGcc  <  2:2606456/148‑1 (MQ=255)
                                                                                                                                                 catcaGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGcc  >  2:761686/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GAGCCTGTTCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCC  >  CP009273/3339893‑3340185

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: