Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,706,777 C→A 8.2% intergenic (+362/‑249) nth → / → dtpA DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,706,7770CA 8.2% 937.2 / 35.2 250intergenic (+362/‑249)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base C (104/125);  new base A (3/18);  total (107/143)
Fisher's exact test for biased strand distribution p-value = 5.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

CAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATAGTGAACATTACC  >  CP009273/1706628‑1706926
                                                                                                                                                     |                                                                                                                                                     
cAATAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCctc                                                                                                                                                       >  2:1633249/1‑150 (MQ=255)
cAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCctc                                                                                                                                                       >  1:1180231/1‑150 (MQ=255)
cAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCctc                                                                                                                                                       <  2:3343749/150‑1 (MQ=255)
cAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCctc                                                                                                                                                       <  2:2894028/150‑1 (MQ=255)
cAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCctc                                                                                                                                                       >  1:1620857/1‑150 (MQ=255)
 tagtaTCAGGTTATATTGATTTAAATTTACACACGCAAAAAATCTTGACCTACAAAACATGACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATTGCAGAACATATCGCCGTTAAGCCATCACCCCCctct                                                                                                                                                      <  1:1291216/146‑1 (MQ=255)
    aaTCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTcag                                                                                                                                                   >  1:843485/1‑150 (MQ=255)
     aTCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTcagt                                                                                                                                                  <  1:4330237/150‑1 (MQ=255)
     aTCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTcagt                                                                                                                                                  >  1:2185213/1‑150 (MQ=255)
         ggTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTcagtgcag                                                                                                                                              <  1:2751729/150‑1 (MQ=255)
          gTTAAAATCATTTAAATTTAAACACGCAACAAATATTGACCTACAAAAAATTACAATTGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGAAGAACATATCGCCGTTAAGCCATAACCCCACTCAcagtgcagt                                                                                                                                             >  1:781588/1‑150 (MQ=255)
            tAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGa                                                                                                                                           <  1:1740179/150‑1 (MQ=255)
             aaaaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGaa                                                                                                                                          >  2:3566283/1‑150 (MQ=255)
             aaaaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGaa                                                                                                                                          >  2:1840997/1‑150 (MQ=255)
             aaaaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGaa                                                                                                                                          >  1:4251022/1‑150 (MQ=255)
               aaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGaaaa                                                                                                                                        >  1:3187405/1‑150 (MQ=255)
                aTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGaaaaa                                                                                                                                       <  2:1816058/150‑1 (MQ=255)
                  cATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAAt                                                                                                                                     <  1:2573683/150‑1 (MQ=255)
                    tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCt                                                                                                                                   <  1:3153350/150‑1 (MQ=255)
                      tAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGc                                                                                                                                 <  1:2499777/150‑1 (MQ=255)
                      aaaaTTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGc                                                                                                                                 <  2:1373465/149‑1 (MQ=255)
                       aaaTTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGcc                                                                                                                                >  1:1832947/1‑150 (MQ=255)
                           ttACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTa                                                                                                                            <  1:2421511/150‑1 (MQ=255)
                           ttACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTa                                                                                                                            <  1:2129884/150‑1 (MQ=255)
                           ttACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTa                                                                                                                            <  2:2008213/150‑1 (MQ=255)
                             acacacGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACg                                                                                                                          >  2:4432717/1‑150 (MQ=255)
                                 acGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttt                                                                                                                      >  2:2418980/1‑150 (MQ=255)
                                  cGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                                                                                                     <  2:3671182/150‑1 (MQ=255)
                                     aacaaATATTGAACTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGa                                                                                                                  <  2:1117438/150‑1 (MQ=255)
                                       caaATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaa                                                                                                                <  2:338859/150‑1 (MQ=255)
                                       caaATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGATTTTTGaaa                                                                                                                <  2:392669/150‑1 (MQ=255)
                                          atatTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAt                                                                                                             <  1:3339546/150‑1 (MQ=255)
                                            atTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAttt                                                                                                           >  1:4415971/1‑150 (MQ=255)
                                            atTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAttt                                                                                                           >  1:2528601/1‑150 (MQ=255)
                                                  cTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                                                                                                     <  2:1004444/150‑1 (MQ=255)
                                                  cTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                                                                                                     >  2:223378/1‑150 (MQ=255)
                                                    aCAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTgg                                                                                                   <  2:1965796/150‑1 (MQ=255)
                                                     cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                                                                                  >  1:3480167/1‑150 (MQ=255)
                                                      aaaaCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGAt                                                                                                 >  2:2087583/1‑150 (MQ=255)
                                                       aaaCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATa                                                                                                >  1:581059/1‑150 (MQ=255)
                                                           aTTACACTGGATATTTTTCAGAAACTGGACTATCTCACGAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAt                                                                                            <  2:4033831/150‑1 (MQ=255)
                                                             tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAttt                                                                                          <  1:4076781/150‑1 (MQ=255)
                                                             tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTATCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGTAAAATTTAACGCTGGATAACAttt                                                                                          <  2:2611517/150‑1 (MQ=255)
                                                             tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTATCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAttt                                                                                          <  2:3431894/150‑1 (MQ=255)
                                                                 cTGGCTATTATTCAGAAACTGGACTATCTCACAAGTAAAAACGCTAAATAGCAGAACATATCGACGTTAAGCCATCACCACCCTCTCAGTTCAGTGAAAAAATCTGAAGATACGTTTTTTGAAAAAAATAACGCTTGATAACATTACCCg                                                                                      >  1:4142019/1‑150 (MQ=255)
                                                                    gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                                                   >  2:2329573/1‑150 (MQ=255)
                                                                     cTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAAt                                                                                  <  2:4199918/150‑1 (MQ=255)
                                                                      tATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                                                                                 >  2:424115/1‑150 (MQ=255)
                                                                      tATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                                                                                 >  1:3847007/1‑150 (MQ=255)
                                                                        tttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGt                                                                               >  2:55426/1‑150 (MQ=255)
                                                                         ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGtt                                                                              <  1:1122248/150‑1 (MQ=255)
                                                                         ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGtt                                                                              >  2:3812319/1‑150 (MQ=255)
                                                                         ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGtt                                                                              <  1:1750997/150‑1 (MQ=255)
                                                                            tCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                                                                           <  2:4399115/150‑1 (MQ=255)
                                                                             cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                                                                          <  1:3390017/150‑1 (MQ=255)
                                                                             cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTATCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                                                                          <  2:83956/150‑1 (MQ=255)
                                                                              aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                                                         >  1:4404878/1‑150 (MQ=255)
                                                                               gAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                                                        <  2:1459952/150‑1 (MQ=255)
                                                                               gAAACTGGACTATCTAACTAGTAAATACGCTAAATAGCAGAAGATAACGCCGTTAAGCCATCACCCGCCTCTCAGTGCAGTGAAAAAATCTTCGGTTACGTATTTTGAAAATTTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                                                        <  1:63631/150‑1 (MQ=255)
                                                                                    tGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGTTTGAATTCCCCGc                                                                   <  1:1981818/150‑1 (MQ=255)
                                                                                     ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                                                  >  2:2980203/1‑150 (MQ=255)
                                                                                      gACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                                                 <  1:2312508/150‑1 (MQ=255)
                                                                                      gACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                                                 <  1:301314/150‑1 (MQ=255)
                                                                                       aCTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTc                                                                >  1:1300391/1‑150 (MQ=255)
                                                                                          aTCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGt                                                             >  1:1559208/1‑150 (MQ=255)
                                                                                          aTCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGt                                                             >  1:3796054/1‑150 (MQ=255)
                                                                                             tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtat                                                          <  1:4026707/150‑1 (MQ=255)
                                                                                               aCTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTatat                                                        <  1:2016865/150‑1 (MQ=255)
                                                                                                cTAGTAAAAACTCTAAATAGCAGAACATATCGCCTTTAAGCCATCACCCCCCTCTGATTTCATTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTATCGCTGGATAACATTTCCAGGAATGGTTGAATTCCCCGCCTCAGTTATATg                                                       <  2:2751559/150‑1 (MQ=255)
                                                                                                 tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGt                                                      <  2:3294508/150‑1 (MQ=255)
                                                                                                    tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAAc                                                   >  2:2684401/1‑150 (MQ=255)
                                                                                                      aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCACTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAg                                                 >  1:3053010/1‑150 (MQ=255)
                                                                                                           gCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGattat                                            >  2:1049896/1‑150 (MQ=255)
                                                                                                            cTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGattatt                                           >  1:1841995/1‑150 (MQ=255)
                                                                                                            cTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGattatt                                           >  1:3247875/1‑150 (MQ=255)
                                                                                                                 tAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACaaa                                      <  2:3663670/150‑1 (MQ=255)
                                                                                                                 tAGCAGAACATATCGCCGTTAAGCCATCACCCACCTCTCAGTGCAGTGAAAAAATCTGACGTTACGTTTTTTGAAAATTTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACaaa                                      <  1:3548364/150‑1 (MQ=255)
                                                                                                                      gAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGATTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACt                                 <  1:997279/150‑1 (MQ=255)
                                                                                                                          atatCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTc                             <  1:2951829/150‑1 (MQ=255)
                                                                                                                          atatCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTc                             <  1:2951812/150‑1 (MQ=255)
                                                                                                                           tatCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCt                            <  1:1356184/150‑1 (MQ=255)
                                                                                                                               gCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGaaa                        >  2:3060370/1‑150 (MQ=255)
                                                                                                                               gCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGaaa                        >  2:1285671/1‑150 (MQ=255)
                                                                                                                               gCCGTTAAGCCATCACCCCCCTATCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGaaa                        <  2:4076279/150‑1 (MQ=255)
                                                                                                                                ccGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAg                       <  1:4502069/150‑1 (MQ=255)
                                                                                                                                 cGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGt                      >  1:2968522/1‑150 (MQ=255)
                                                                                                                                    tAAGCCATCACCCCCCGATCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCa                   <  2:714626/150‑1 (MQ=255)
                                                                                                                                     aaGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCaa                  >  1:2084418/1‑150 (MQ=255)
                                                                                                                                      aGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAg                 <  2:497765/150‑1 (MQ=255)
                                                                                                                                      aGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAGTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAg                 <  2:3169011/150‑1 (MQ=255)
                                                                                                                                        ccATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGAt               <  2:3745380/150‑1 (MQ=255)
                                                                                                                                         cATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATa              <  1:1648060/150‑1 (MQ=255)
                                                                                                                                         cATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATa              <  1:4200132/150‑1 (MQ=255)
                                                                                                                                           tCACCCCCCTCTCAGTGCAGTGAACAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGATTGGTTGAATTCCCCGCCTCCTTTATCTGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATatt            <  2:619244/150‑3 (MQ=255)
                                                                                                                                             aCCCCCCTATCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATAGTGa          >  2:1201143/1‑150 (MQ=255)
                                                                                                                                              ccccccTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATAGTGaa         <  2:447201/150‑1 (MQ=255)
                                                                                                                                                ccccTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATAGTGAACa       >  1:4044287/1‑150 (MQ=255)
                                                                                                                                                 cccTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATAGTGAACAt      >  2:241601/1‑150 (MQ=255)
                                                                                                                                                   ctctCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATAGTGAACATTa    >  2:374624/1‑150 (MQ=255)
                                                                                                                                                    tctcAGTGCAGTGAAAAAATCTGCCGTTACATTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATAGTGAACATTAc   >  1:2378948/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGAAAGTGCAAGATAGTGAACATTACC  >  CP009273/1706628‑1706926

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