Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,409,174 T→G 5.3% intergenic (+79/‑101) rnr → / → rlmB exoribonuclease R, RNase R/23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,1740TG 5.3% 267.8 / 6.1 76intergenic (+79/‑101)rnr/rlmBexoribonuclease R, RNase R/23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base G (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGA  >  NC_000913/4409027‑4409321
                                                                                                                                                   |                                                                                                                                                   
gCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                    >  3:803826/1‑149 (MQ=255)
 cGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                   <  2:280681/149‑1 (MQ=255)
  gAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                                  >  1:621502/1‑149 (MQ=255)
      aCGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                              <  1:330355/149‑1 (MQ=255)
       cGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                             <  1:725717/149‑1 (MQ=39)
        gCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                            <  2:202842/149‑1 (MQ=255)
        gCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATAAGTCGTTCACGCCGCATCCGGca                                                                                                                                            <  2:69534/149‑2 (MQ=255)
         cAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                           <  4:436335/149‑1 (MQ=255)
                 aGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGGGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATACTc                                                                                                                                   <  4:176460/149‑1 (MQ=255)
                   cTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCAc                                                                                                                                 <  4:803826/149‑1 (MQ=255)
                   cTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCAc                                                                                                                                 <  1:341393/149‑1 (MQ=255)
                        gCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaaa                                                                                                                            <  2:282383/149‑1 (MQ=255)
                          gACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAtt                                                                                                                          <  1:718507/149‑1 (MQ=255)
                             cAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                                                       >  2:603760/1‑149 (MQ=255)
                             cAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTCCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                                                       <  1:782328/149‑1 (MQ=255)
                               aagcgaagcgTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                                                                                                                     >  3:445906/1‑149 (MQ=255)
                                    aagcgTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAACAGGTTTAGATTGCTGACAAAATGCGCTTTGTTCCTGCCGGATGCGGTGTGACCGTCTTATCCGGCCTACATAATCACGCAAATGCAATATATTg                                                                                                                >  2:208856/1‑149 (MQ=255)
                                        gTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCACTATATTGCaga                                                                                                            <  3:764780/149‑1 (MQ=255)
                                        gTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga                                                                                                            <  4:638043/149‑1 (MQ=255)
                                           cggcgAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGAt                                                                                                         <  4:364746/149‑1 (MQ=255)
                                               gaagaaAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg                                                                                                     <  3:804562/149‑1 (MQ=255)
                                                   aaaaaaGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAgg                                                                                                 >  1:493665/1‑149 (MQ=255)
                                                     aaaaGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc                                                                                               >  2:172811/1‑149 (MQ=255)
                                                              agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg                                                                                      >  3:441322/1‑149 (MQ=255)
                                                               gagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                         >  2:79831/1‑97 (MQ=255)
                                                               gagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                         <  1:79831/97‑1 (MQ=255)
                                                               gagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt                                                                                     >  3:725762/1‑149 (MQ=255)
                                                                agTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTa                                                                                    <  1:886521/149‑1 (MQ=255)
                                                                    aTCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGc                                                                                <  3:740872/149‑1 (MQ=255)
                                                                    aTCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGc                                                                                >  4:772275/1‑149 (MQ=255)
                                                                       aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATc                                                                             <  2:296237/149‑1 (MQ=255)
                                                                          aCCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAgg                                                                          >  4:262215/1‑149 (MQ=255)
                                                                            ccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCa                                                                        >  4:587502/1‑149 (MQ=255)
                                                                             ctctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCaa                                                                       >  2:869663/1‑149 (MQ=255)
                                                                              tctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaaa                                                                                                                            >  1:705723/1‑95 (MQ=255)
                                                                              tctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaaa                                                                                                                            <  2:705723/95‑1 (MQ=255)
                                                                                tttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCTTCAGGCAAAAc                                                                    >  3:572966/1‑149 (MQ=255)
                                                                                     aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGc                                                                                >  1:895895/1‑132 (MQ=255)
                                                                                     aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGc                                                                                <  2:895895/132‑1 (MQ=255)
                                                                                      agaagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaac                                                              <  2:800636/149‑1 (MQ=255)
                                                                                         agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaacaac                                                           >  2:712506/1‑149 (MQ=255)
                                                                                           agGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGa                                                         <  3:262215/149‑1 (MQ=255)
                                                                                           agGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGa                                                         <  3:234871/149‑1 (MQ=255)
                                                                                            gggTTAGATTGCTGACAAAGTGCGCTTTGATCATGCCGGATGCTGCGTGAACGGCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCACAGACCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAg                                                        <  1:208856/149‑1 (MQ=21)
                                                                                             ggTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                             <  1:626867/63‑1 (MQ=255)
                                                                                             ggTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                             >  2:626867/1‑63 (MQ=255)
                                                                                             ggTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGa                                                         <  1:449598/147‑1 (MQ=255)
                                                                                             ggTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGa                                                         >  2:449598/1‑147 (MQ=255)
                                                                                              gTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTa                                                      >  3:672150/1‑149 (MQ=255)
                                                                                              gTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTa                                                      >  1:266087/1‑149 (MQ=255)
                                                                                                tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACa                                                    >  2:435689/1‑149 (MQ=255)
                                                                                                       cTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGa                                             <  1:869663/149‑1 (MQ=255)
                                                                                                       cTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGa                                             >  2:702803/1‑149 (MQ=255)
                                                                                                         gACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTa                                                                                    <  2:828676/108‑1 (MQ=255)
                                                                                                         gACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTa                                                                                    >  1:828676/1‑108 (MQ=255)
                                                                                                         gACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGc                                           >  4:109259/1‑149 (MQ=255)
                                                                                                         gACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGc                                           <  3:531285/149‑1 (MQ=255)
                                                                                                          aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                                                       >  1:171997/1‑72 (MQ=255)
                                                                                                          aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                                                       <  2:171997/72‑1 (MQ=255)
                                                                                                          aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAg                                            >  1:414283/1‑147 (MQ=255)
                                                                                                          aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAg                                            <  2:414283/147‑1 (MQ=255)
                                                                                                               aTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATg                                     >  2:600375/1‑149 (MQ=255)
                                                                                                                     tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATAAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTAc                               >  1:598777/1‑149 (MQ=255)
                                                                                                                      ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACg                              <  1:12911/149‑1 (MQ=255)
                                                                                                                      ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACg                              <  2:266087/149‑1 (MQ=255)
                                                                                                                       tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACgg                             <  2:24000/149‑1 (MQ=255)
                                                                                                                       tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACgg                             <  1:295979/149‑1 (MQ=255)
                                                                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg                                                                           >  3:442605/3‑98 (MQ=18)
                                                                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACg                     >  1:876044/1‑149 (MQ=37)
                                                                                                                                              aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCctgc      >  3:438106/1‑149 (MQ=255)
                                                                                                                                               cGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCctgct     <  2:62027/149‑1 (MQ=255)
                                                                                                                                                gCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCctgctg    >  4:182560/1‑149 (MQ=255)
                                                                                                                                                  cTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGa  >  3:842146/1‑149 (MQ=255)
                                                                                                                                                  cTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGa  >  2:225477/1‑149 (MQ=255)
                                                                                                                                                  cTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGa  >  1:349023/1‑149 (MQ=255)
                                                                                                                                                  cTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGa  >  4:710813/1‑149 (MQ=255)
                                                                                                                                                  cTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGa  >  4:772872/1‑149 (MQ=255)
                                                                                                                                                   ttATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGATCGAAATGATTTACGGCATc                         >  2:390432/1‑125 (MQ=255)
                                                                                                                                                   ttATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGATCGAAATGATTTACGGCATc                         <  1:390432/125‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGA  >  NC_000913/4409027‑4409321

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