Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,714,372 G→T 6.7% intergenic (+142/+67) srmB → / ← yfiE ATP‑dependent RNA helicase SrmB/putative LysR‑type DNA‑binding transcriptional regulator YfiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,714,3720GT 6.7% 524.3 / 15.5 156intergenic (+142/+67)srmB/yfiEATP‑dependent RNA helicase SrmB/putative LysR‑type DNA‑binding transcriptional regulator YfiE
Reads supporting (aligned to +/- strand):  ref base G (83/62);  new base T (2/9);  total (85/71)
Fisher's exact test for biased strand distribution p-value = 2.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 162 total aligned reads displayed.

AAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGCATGATGTG  >  NC_000913/2714223‑2714518
                                                                                                                                                     |                                                                                                                                                  
aagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                                    <  1:1454447/150‑1 (MQ=255)
     tAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                               <  2:651182/150‑1 (MQ=25)
      aaTCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                              <  2:726990/150‑1 (MQ=255)
       aTCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                             <  1:1982141/150‑1 (MQ=25)
       aTCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                             <  1:191306/150‑1 (MQ=25)
       aTCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgt                                                                                                                                             <  2:1959249/150‑2 (MQ=255)
         ctcAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                           >  2:2208087/1‑150 (MQ=255)
          tcAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCAGCATCCGGCAt                                                                                                                                          >  1:1636544/1‑150 (MQ=255)
           cAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGccta                                                                                                                                         <  1:767546/150‑4 (MQ=21)
            aaTGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                                                                                                        >  1:933669/1‑150 (MQ=255)
               gCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATa                                                                                                                                     >  1:1359953/1‑150 (MQ=255)
               gCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATa                                                                                                                                     >  1:2041741/1‑150 (MQ=255)
                ccAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATaa                                                                                                                                    <  2:862638/150‑1 (MQ=255)
                 cAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATaac                                                                                                                                   >  1:1850000/1‑149 (MQ=255)
                  aGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTTTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAg                                                                                                                                  <  2:2008670/150‑1 (MQ=255)
                    gTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGATTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAgcg                                                                                                                                >  2:385623/1‑150 (MQ=255)
                     tttAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGAATAAAgcgc                                                                                                                               <  2:1162163/150‑1 (MQ=255)
                     tttAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAgcgc                                                                                                                               <  1:443265/150‑1 (MQ=255)
                         aGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACtt                                                                                                                           <  2:924075/150‑1 (MQ=255)
                            cTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGt                                                                                                                        >  2:810636/1‑150 (MQ=255)
                              ggTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTTTAGGCCGTATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCg                                                                                                                      <  2:414302/150‑1 (MQ=255)
                              ggTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCg                                                                                                                      <  2:1765614/150‑1 (MQ=255)
                                 aTTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGTCa                                                                                                                   <  2:2138480/150‑1 (MQ=255)
                                 aTTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCACGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACa                                                                                                                   >  2:1001115/1‑150 (MQ=255)
                                   tAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAg                                                                                                                 <  2:1520691/150‑1 (MQ=255)
                                          gCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATg                                                                                                          >  2:2186033/1‑150 (MQ=255)
                                          gCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATg                                                                                                          <  1:2282560/150‑1 (MQ=255)
                                               cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCggg                                                                                                     >  2:1442037/1‑150 (MQ=255)
                                               cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCggg                                                                                                     >  2:348218/1‑150 (MQ=255)
                                                 aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGtt                                                                                                   <  1:1997806/150‑1 (MQ=255)
                                                   cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTa                                                                                                 >  2:1656641/1‑150 (MQ=255)
                                                     cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACa                                                                                                                   >  1:1668517/1‑130 (MQ=255)
                                                     cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACa                                                                                                                   <  2:1668517/130‑1 (MQ=255)
                                                     cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAAt                                                                                               >  1:1155883/1‑150 (MQ=255)
                                                     cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAAt                                                                                               >  2:235358/1‑150 (MQ=255)
                                                         aCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCg                                                                                           >  1:846184/1‑150 (MQ=255)
                                                               tGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACatta                                                                                     >  2:1217926/1‑150 (MQ=255)
                                                                 aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACattatt                                                                                   >  1:1161493/1‑150 (MQ=255)
                                                                 aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACattatt                                                                                   >  1:2251188/1‑150 (MQ=255)
                                                                  tGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTAttt                                                                                  >  2:202801/1‑150 (MQ=255)
                                                                   gcgcTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTAtttt                                                                                 <  1:1524587/150‑1 (MQ=255)
                                                                   gcgcTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTAtttt                                                                                 <  2:1047001/150‑1 (MQ=255)
                                                                    cgcTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTAttttt                                                                                >  2:275799/1‑150 (MQ=255)
                                                                    cgcTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTAttttt                                                                                >  1:2288093/1‑150 (MQ=255)
                                                                    cgcTACGCTTATCAGGCATACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCAGACATTAttttt                                                                                >  2:788039/1‑150 (MQ=255)
                                                                    cccTACGCTTATCAGGCCTCCGTGGCTCATGCAATATATTGAATTTGCCCGATCTTGTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGATAAATGCCGGGTTTAATCCCGACATTAttttt                                                                                >  2:1211534/3‑150 (MQ=255)
                                                                            ttATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGcc                                                                        <  2:2484952/150‑1 (MQ=255)
                                                                            ttATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGcc                                                                        >  1:1709732/1‑150 (MQ=255)
                                                                              aTCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCcac                                                                      <  1:53413/150‑1 (MQ=255)
                                                                               tCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCcaca                                                                     <  2:1119986/150‑1 (MQ=255)
                                                                                cAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACaa                                                                    >  2:1993197/1‑150 (MQ=255)
                                                                                  ggCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACaaaa                                                                  >  1:1950592/1‑150 (MQ=255)
                                                                                   gCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAAc                                                                 <  2:1466839/150‑1 (MQ=255)
                                                                                     cTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAActc                                                               >  2:713164/1‑150 (MQ=255)
                                                                                       aCGTGGCTCACGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAActctc                                                             <  1:1001115/150‑1 (MQ=255)
                                                                                        cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAActctct                                                            <  1:2230594/150‑1 (MQ=255)
                                                                                          gggCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTc                                                          <  2:1952467/149‑1 (MQ=255)
                                                                                             cTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAac                                                       >  2:1336215/1‑150 (MQ=255)
                                                                                             cTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAac                                                       >  2:883151/1‑150 (MQ=255)
                                                                                                 tGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAacacac                                                   <  1:413285/150‑1 (MQ=255)
                                                                                                  gCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACt                                                  >  1:2389554/1‑150 (MQ=255)
                                                                                                  gCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACt                                                  >  1:350822/1‑150 (MQ=255)
                                                                                                   cAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACTCACTg                                                 <  2:1359953/150‑1 (MQ=255)
                                                                                                   cAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTg                                                 <  2:623868/150‑1 (MQ=255)
                                                                                                    aaTATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGa                                                >  2:1788086/1‑150 (MQ=255)
                                                                                                     atatatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGaa                                               >  2:985802/1‑150 (MQ=255)
                                                                                                      tatatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAAt                                              <  2:2309852/150‑1 (MQ=255)
                                                                                                         atTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATaaa                                           >  1:243919/1‑150 (MQ=255)
                                                                                                         atTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATaaa                                           <  2:558618/150‑1 (MQ=255)
                                                                                                             aaTTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAgtg                                       >  1:1414942/1‑150 (MQ=255)
                                                                                                               tttGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGtgtg                                     <  1:342425/150‑1 (MQ=255)
                                                                                                                 tGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCa                                   >  1:2010352/1‑150 (MQ=255)
                                                                                                                  gCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAActctct                                                            >  1:1235813/1‑124 (MQ=255)
                                                                                                                  gCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAActctct                                                            <  2:1235813/124‑1 (MQ=255)
                                                                                                                   cACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATc                                 <  1:833015/150‑1 (MQ=255)
                                                                                                                    aCGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCg                                >  1:2137195/1‑150 (MQ=255)
                                                                                                                       aTCTTGTAGGCCGGATAAGGCGTTCAAGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCg                             >  1:1597702/1‑150 (MQ=255)
                                                                                                                          ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGc                          <  1:2440102/150‑1 (MQ=255)
                                                                                                                          ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGc                          <  1:2032931/150‑1 (MQ=255)
                                                                                                                          ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGc                          <  2:2074743/150‑1 (MQ=255)
                                                                                                                           ttcatgCCGGATGCGGCGTGAACGCCTTATCCGGCATGATTAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACTCTGAATAAAAGTGTGCATCGCCGGGCt                         >  1:636786/6‑150 (MQ=38)
                                                                                                                           tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCt                         >  1:1819447/1‑150 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAAcc                     >  2:589941/1‑150 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAAcc                     >  2:1472234/1‑150 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAAcc                     <  1:577975/150‑1 (MQ=255)
                                                                                                                                 ccGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGc                   >  2:522253/1‑150 (MQ=255)
                                                                                                                                 ccGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGc                   >  2:2004684/1‑150 (MQ=255)
                                                                                                                                  ctgATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCt                  <  1:788039/148‑1 (MQ=255)
                                                                                                                                  cGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCt                  <  1:275799/150‑1 (MQ=255)
                                                                                                                                       aaGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTcc             <  2:526928/150‑1 (MQ=255)
                                                                                                                                         ggCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAg           >  2:1377977/1‑150 (MQ=255)
                                                                                                                                         ggCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATAGCCGGGCTAACCGCTTTTCCAg           >  1:157308/1‑150 (MQ=255)
                                                                                                                                          ttattCATGCCGGATGCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGc          >  2:1695615/4‑150 (MQ=255)
                                                                                                                                          gCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGc          <  1:402718/150‑1 (MQ=255)
                                                                                                                                           cGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGCa         >  2:1344051/1‑150 (MQ=255)
                                                                                                                                            gtgaaCGCCTTATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTTATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGCat        <  1:1048897/146‑1 (MQ=255)
                                                                                                                                                aCGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGCatgatg    <  1:1329727/150‑1 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGCATGAtgt   <  2:2389554/150‑1 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGCATGAtgt   >  1:998715/1‑150 (MQ=255)
                                                                                                                                                 cgccgcATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGCATGAtgt   >  1:1749993/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCCAGCATGATGTG  >  NC_000913/2714223‑2714518

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