Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,011,932 G→A 6.6% intergenic (+61/‑69) uup → / → pqiA replication regulatory ABC‑F family DNA‑binding ATPase/paraquat‑inducible, SoxRS‑regulated inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,011,9320GA 6.6% 281.1 / 7.2 90intergenic (+61/‑69)uup/pqiAreplication regulatory ABC‑F family DNA‑binding ATPase/paraquat‑inducible, SoxRS‑regulated inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (42/42);  new base A (4/2);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

TTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCA  >  NC_000913/1011786‑1012080
                                                                                                                                                  |                                                                                                                                                    
ttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGt                                                                                                                                                    <  2:1461678/149‑1 (MQ=255)
  gctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCTTTAAAAAATGGTGGCTGATCGCAATATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGttt                                                                                                                                                  <  2:1034990/149‑1 (MQ=255)
    tGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTAt                                                                                                                                                <  2:925068/149‑1 (MQ=255)
      atatGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCa                                                                                                                                              >  2:1548663/1‑149 (MQ=255)
       tatGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAg                                                                                                                                             <  2:255451/149‑1 (MQ=255)
         tGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGc                                                                                                                                           <  2:479384/149‑1 (MQ=255)
          ggCTGCTGCAGTGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCa                                                                                                                                          >  1:268604/1‑149 (MQ=255)
          ggCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCa                                                                                                                                          <  2:1504868/149‑1 (MQ=255)
          ggCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCa                                                                                                                                          <  1:485058/149‑1 (MQ=255)
          ggCTGCTGCAGAGCAGGAGCGGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGCa                                                                                                                                          <  1:308273/149‑1 (MQ=255)
           gctgctGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCAt                                                                                                                                         <  2:650330/149‑1 (MQ=255)
           gctgctGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGCaa                                                                                                                                         <  2:26475/149‑2 (MQ=255)
             tgctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGa                                                                                                                                       >  1:1141843/1‑149 (MQ=255)
                tgcAGAGCAGGAGCGGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGaaaa                                                                                                                                    <  2:249325/149‑1 (MQ=255)
                  cAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAcc                                                                                                                                  >  2:249827/1‑149 (MQ=255)
                       cAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAacac                                                                                                                             >  2:1117078/1‑149 (MQ=255)
                           aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTg                                                                                                                         <  1:456645/149‑1 (MQ=255)
                           aGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGATAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTg                                                                                                                         <  2:1140145/149‑1 (MQ=255)
                            gCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGt                                                                                                                        <  2:1410751/149‑1 (MQ=255)
                             cTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGtt                                                                                                                       <  2:1519196/149‑1 (MQ=255)
                              tGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTa                                                                                                                      <  1:1219791/149‑1 (MQ=255)
                                gAGCAAGCCTTTGCACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAAGCATGAAAACCAACACATTCTTAGt                                                                                                                    >  1:276828/1‑149 (MQ=255)
                                gAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGt                                                                                                                    >  1:1515034/1‑149 (MQ=255)
                                gAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGt                                                                                                                    >  1:1098237/1‑149 (MQ=255)
                                gAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGt                                                                                                                    >  2:174993/1‑149 (MQ=255)
                                gAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGAAAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGt                                                                                                                    >  2:798669/1‑149 (MQ=255)
                                  gCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCa                                                                                                                  >  2:322009/1‑149 (MQ=255)
                                  gCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCa                                                                                                                  >  1:457719/1‑149 (MQ=255)
                                  gCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCa                                                                                                                  >  2:164499/1‑149 (MQ=255)
                                  gCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCa                                                                                                                  <  1:1115008/149‑1 (MQ=255)
                                   cAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCAt                                                                                                                 >  2:120149/1‑149 (MQ=255)
                                         tttGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATTAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTaa                                                                                                           <  1:914042/149‑1 (MQ=255)
                                         tttGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTaa                                                                                                           >  1:166727/1‑149 (MQ=255)
                                               cGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTTCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAgc                                                                                                     <  1:1155533/149‑1 (MQ=255)
                                                  tGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTATAATAGCGca                                                                                                  >  2:1050909/1‑149 (MQ=255)
                                                    ggAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCGTCTGTAAAATAGCGcatc                                                                                                >  1:632278/1‑149 (MQ=255)
                                                       gTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCAtt                                                                                             >  2:197605/1‑149 (MQ=255)
                                                        tATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACATAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGATTATAAGGCATGAAAACCAACACTTTGTTAGACATCTGTAAAATACCGCATCATTa                                                                                            >  2:783548/1‑149 (MQ=255)
                                                              gAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCTTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGt                                                                                      >  2:1238315/1‑149 (MQ=255)
                                                                aGCGTTAAAAAATGGTGGCTGATCACAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGGTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTAc                                                                                    <  1:976734/149‑1 (MQ=255)
                                                                    ttAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAt                                                                                >  2:584442/1‑149 (MQ=255)
                                                                    ttAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAAAAGCTCATCATTAAGGAGTACCAAt                                                                                >  1:779796/1‑149 (MQ=255)
                                                                           aTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACTGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGaa                                                                         <  1:174993/149‑1 (MQ=255)
                                                                           aTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGaa                                                                         >  2:575680/1‑149 (MQ=255)
                                                                                 gCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcatcat                                                                   <  2:1515034/149‑1 (MQ=255)
                                                                                 gCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcatcat                                                                   <  2:1098237/149‑1 (MQ=255)
                                                                                   tGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAt                                                                                <  2:703973/134‑1 (MQ=255)
                                                                                   tGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAt                                                                                >  1:703973/1‑134 (MQ=255)
                                                                                       ctaaaaaTATTATATTTAAACAGACACAGCTGACAACGGTAAAATTGCCTGATGAACTTGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGcc                                                             <  1:1419173/146‑1 (MQ=255)
                                                                                                gTCGATTTAAACAGATACAGCTTAAAACTGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAATcac                                                    <  1:680596/149‑1 (MQ=255)
                                                                                                 tCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGcaca                                                   >  1:934832/1‑149 (MQ=255)
                                                                                                 tCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGcaca                                                   >  2:1571464/1‑149 (MQ=255)
                                                                                                     tttAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGcacaggat                                               <  1:164499/149‑5 (MQ=255)
                                                                                                     tttAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCt                                               >  1:819315/1‑149 (MQ=255)
                                                                                                     tttAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACAAAATCATGCCGCGAATCACATcat                                               >  1:187827/1‑147 (MQ=255)
                                                                                                     tttAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCAGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCt                                               >  1:835477/1‑149 (MQ=255)
                                                                                                      ttAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCtg                                              <  1:120149/149‑1 (MQ=255)
                                                                                                         aaCAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGc                                           >  2:13627/1‑149 (MQ=255)
                                                                                                         aaCAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGc                                           <  1:575680/149‑1 (MQ=255)
                                                                                                          aCAGAGACAGCTGACACCGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCt                                          <  2:300214/149‑1 (MQ=255)
                                                                                                          aCAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCt                                          >  1:929386/1‑149 (MQ=255)
                                                                                                          aCAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACAAACACTTTGTTAGTCATCTGTATAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCt                                          >  2:981920/1‑149 (MQ=255)
                                                                                                              agaCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACAAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCa                                      >  2:1180689/1‑149 (MQ=255)
                                                                                                               gaCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCGTCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAGCACAGGATGTGCttcgcgg                                    <  2:632278/149‑7 (MQ=255)
                                                                                                               gaCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAg                                     <  2:268604/149‑1 (MQ=255)
                                                                                                                aCAGCTGACAACGGTAAAATTGCCTGATGAACTGGTTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAgt                                    <  2:1415159/149‑1 (MQ=255)
                                                                                                                aCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAgt                                    <  2:457719/149‑1 (MQ=255)
                                                                                                                  aGCTGACAACGGTAAATTTGCCTTATGAACTGGGTTTATCAGGCATGAAAACCTACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAgtgt                                  <  1:624804/149‑1 (MQ=255)
                                                                                                                    cTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGa                                >  1:653360/1‑149 (MQ=255)
                                                                                                                     tGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGAc                               <  2:653360/149‑1 (MQ=255)
                                                                                                                     tGACAACGGTAAAATTGCCTGATGAACTGGGTTTATAAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCAACATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGAc                               >  1:1090923/1‑149 (MQ=255)
                                                                                                                        cAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATg                            <  1:1117078/149‑1 (MQ=255)
                                                                                                                         aaCGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGc                           <  2:160149/149‑1 (MQ=255)
                                                                                                                          aCGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCt                          <  1:1548663/149‑1 (MQ=255)
                                                                                                                             gTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGACGCACATCCTGTGCTCGCAGTGTGACATGCtggt                       >  1:760096/1‑149 (MQ=255)
                                                                                                                             gTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTGTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAATCACATCGTGTGCTCGCAGTGTGACATGCtggt                       >  1:295321/1‑149 (MQ=255)
                                                                                                                              ttttcatGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtggtg                      >  2:803410/7‑149 (MQ=255)
                                                                                                                                aaaTTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGc                    <  2:166727/149‑1 (MQ=255)
                                                                                                                                  aTTGCTTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGt                  <  1:1238315/149‑1 (MQ=255)
                                                                                                                                    tGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTa                <  2:1411621/149‑1 (MQ=255)
                                                                                                                                    tGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTa                >  2:582727/1‑149 (MQ=255)
                                                                                                                                    tGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTa                >  1:1080471/1‑149 (MQ=255)
                                                                                                                                     gCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTAc               >  1:787935/1‑149 (MQ=255)
                                                                                                                                      ccTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTAcc              >  1:418909/1‑149 (MQ=255)
                                                                                                                                       cTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcg             >  2:1067014/1‑149 (MQ=255)
                                                                                                                                        tgatgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgc            >  1:1541180/1‑149 (MQ=255)
                                                                                                                                         gatgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgcg           >  1:1199854/1‑149 (MQ=255)
                                                                                                                                         gatgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgcg           >  1:1024141/1‑149 (MQ=255)
                                                                                                                                               cTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGa     <  1:322009/149‑1 (MQ=255)
                                                                                                                                               cTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGa     <  2:929386/149‑1 (MQ=255)
                                                                                                                                                  ggTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCa  >  1:373137/1‑149 (MQ=255)
                                                                                                                                                  ggTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCa  >  2:7469/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAGCA  >  NC_000913/1011786‑1012080

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