Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 899,184 G→A A223T (GCC→ACC)  rlmC → 23S rRNA m(5)U747 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913899,1840GA97.8% 289.8 / ‑2.7 90A223T (GCC→ACC) rlmC23S rRNA m(5)U747 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (48/40);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CAACAGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGAC  >  NC_000913/899088‑899282
                                                                                                |                                                                                                  
caacaGGCACTGGCGGAGCGTTTTAATGACGTACCGCCGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCCCCCGCCAGTTGTACACCa                                                                                                 >  2:1213157/1‑100 (MQ=255)
 aacaGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTAcaccac                                                                                                >  1:1310717/1‑100 (MQ=255)
 aacaGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTACCAGACTAATCCGGCGGTCGCCAGCCAGTTGTAcaccac                                                                                                >  2:1635002/1‑100 (MQ=255)
  acaGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACaccacc                                                                                               >  2:1711283/1‑100 (MQ=255)
  acaGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACaccacc                                                                                               >  1:1933280/1‑100 (MQ=255)
  acaGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACaccacc                                                                                               >  2:930995/1‑100 (MQ=255)
     ggCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACcgcg                                                                                            >  2:186247/1‑100 (MQ=255)
     ggCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACcgcg                                                                                            >  2:186258/1‑100 (MQ=255)
      gCACTGGCGGAGCGTTTTAATGACGGACCGCTGTGGATCCGTCCGCAGAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACcgcgc                                                                                           >  1:2365446/1‑100 (MQ=255)
      gCACTGGCGGAGCGTTTTAATGACGGACCGCTGTGGATCCGTCCGCAGAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACagcgc                                                                                           >  1:2365437/1‑100 (MQ=255)
       cACTGGCGGAGCGTTTTAATGACGTACCGTTGTGGTTCCGTCCGCAAAGTGTCTTCCAGACTAAGCCGGCGGTCGCCAGCCAGTTGTACACCACcgcgcg                                                                                          <  2:2189598/100‑1 (MQ=255)
           ggCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGAc                                                                                      >  2:1912093/1‑100 (MQ=255)
           ggCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGAc                                                                                      <  1:1550246/100‑1 (MQ=255)
           ggCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGAc                                                                                      <  2:1883718/100‑1 (MQ=255)
            gCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACt                                                                                     >  2:1717780/1‑100 (MQ=255)
              ggAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTgg                                                                                   <  1:1749915/100‑1 (MQ=255)
               gAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTggg                                                                                  <  1:1153882/100‑1 (MQ=255)
                aGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGt                                                                                 <  2:308703/100‑1 (MQ=255)
                 gCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTa                                                                                >  2:531501/1‑100 (MQ=255)
                  cGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTAc                                                                               >  2:337798/1‑100 (MQ=255)
                   gTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACg                                                                              >  2:1056171/1‑100 (MQ=255)
                   gTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACg                                                                              >  1:477551/1‑100 (MQ=255)
                   gTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACg                                                                              >  1:2194636/1‑100 (MQ=255)
                     tttAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGAc                                                                            >  1:553781/1‑100 (MQ=255)
                       tAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAg                                                                          <  1:501219/100‑1 (MQ=255)
                        aaTGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGc                                                                         >  1:814025/1‑100 (MQ=255)
                          tGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTg                                                                       <  1:1140060/100‑1 (MQ=255)
                           gACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGc                                                                      >  1:2233744/1‑100 (MQ=255)
                            aCGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGcc                                                                     <  1:425915/100‑1 (MQ=255)
                            aCGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGcc                                                                     <  1:1463707/100‑1 (MQ=255)
                             cGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCg                                                                    >  1:116688/1‑100 (MQ=255)
                             cGTACCGCTGTGGATCCGTCCGCAAAGGTTCTTCCAGACTAATCCGGCGGTCCCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCg                                                                    >  2:1744029/1‑100 (MQ=255)
                              ttACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACGAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCgg                                                                   <  2:1668685/99‑1 (MQ=255)
                              gTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCgg                                                                   <  1:2538200/100‑1 (MQ=255)
                               tACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGt                                                                  <  1:403890/100‑1 (MQ=255)
                               tACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGt                                                                  <  2:1732820/100‑1 (MQ=255)
                               tACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGt                                                                  <  2:98082/100‑1 (MQ=255)
                               tACCGCTGTGGATCCGTCCGCAAAGTTTCGTCCAGACTACTCCGGCGGTCGGCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGt                                                                  <  2:98104/100‑1 (MQ=255)
                                 ccGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTa                                                                >  2:71604/1‑100 (MQ=255)
                                    cTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAAc                                                             <  2:796637/100‑1 (MQ=255)
                                    cTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAAc                                                             <  1:632166/100‑1 (MQ=255)
                                    cTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAAc                                                             <  2:2017950/100‑1 (MQ=255)
                                     tgtgGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACa                                                            >  2:793383/1‑100 (MQ=255)
                                     tgtgGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACa                                                            <  2:458720/100‑1 (MQ=255)
                                      gtgGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACat                                                           <  2:783494/100‑1 (MQ=255)
                                         gATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATAtg                                                        >  1:2070311/1‑100 (MQ=255)
                                          aTCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGGCGGCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATAtgt                                                       <  2:101411/100‑1 (MQ=255)
                                            ccgtccgCTAAGTTTCTTCCAGACTAATCCGGCGGTCGTCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTgg                                                     <  1:2522888/100‑1 (MQ=255)
                                            ccgtccgCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTgg                                                     <  1:1196519/100‑1 (MQ=255)
                                             cgtccgCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTggg                                                    <  2:560470/100‑1 (MQ=255)
                                             cgtccgCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATACGTggg                                                    >  2:886079/1‑100 (MQ=255)
                                                ccgCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATc                                                 >  1:2161177/1‑100 (MQ=255)
                                                  gCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTg                                               >  2:989856/1‑100 (MQ=255)
                                                   cAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGt                                              >  1:699115/1‑100 (MQ=255)
                                                     aaGTTTCTTCCAGCCTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCTCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTc                                            <  1:2380870/100‑1 (MQ=255)
                                                       gTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTg                                          >  1:1906977/1‑100 (MQ=255)
                                                         ttcttcCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCg                                        >  1:26742/1‑100 (MQ=255)
                                                           cttcCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGgt                                      >  2:495188/1‑100 (MQ=255)
                                                           cttcCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGgt                                      >  1:38759/1‑100 (MQ=255)
                                                            ttcCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGgtg                                     >  1:1033004/1‑100 (MQ=255)
                                                             tcCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGgtgt                                    >  1:1292121/1‑100 (MQ=255)
                                                             gccAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGgtgt                                    >  1:2514286/2‑100 (MQ=255)
                                                              ccAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGtgtg                                   >  2:1223337/1‑100 (MQ=255)
                                                              ccAGACTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGtgtg                                   >  2:1439587/1‑100 (MQ=255)
                                                                   cTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTgggggg                              >  1:1899824/1‑100 (MQ=255)
                                                                   cTAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTgggggg                              >  2:624596/1‑100 (MQ=255)
                                                                    tAATCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGCGGGc                             <  2:121314/100‑1 (MQ=255)
                                                                     aaTCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCt                            >  2:1849871/1‑100 (MQ=255)
                                                                      aTCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCtt                           >  2:2227063/1‑100 (MQ=255)
                                                                      aTCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCtt                           >  2:1484313/1‑100 (MQ=255)
                                                                       tCCGGCGGTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCttt                          <  1:344379/100‑1 (MQ=255)
                                                                          ggcggTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGt                       <  2:2074764/100‑1 (MQ=255)
                                                                             ggTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTa                    <  1:1720478/100‑1 (MQ=255)
                                                                              gTCGGCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTac                   <  2:2207153/100‑1 (MQ=255)
                                                                              gTCGCCAGCCGGTTGTGCACCACCGCGTGCGAGTGGGTATGGGAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGGGGGGGGCTTTGGTTTac                   <  1:2375760/100‑1 (MQ=255)
                                                                              gTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTac                   <  2:1302774/100‑1 (MQ=255)
                                                                              gTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTac                   <  1:1491039/100‑1 (MQ=255)
                                                                              gTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTac                   <  2:352537/100‑1 (MQ=255)
                                                                              gTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTac                   <  2:504248/100‑1 (MQ=255)
                                                                              gTCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTac                   <  2:484590/100‑1 (MQ=255)
                                                                               tCGCCAGCCAGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTaca                  >  1:2014496/1‑100 (MQ=255)
                                                                                 gccagccaGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATGTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACt                >  1:2617542/1‑100 (MQ=255)
                                                                                 gccagccaGTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACt                >  1:1981312/1‑100 (MQ=255)
                                                                                  ccagccagTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTg               <  2:1089200/100‑1 (MQ=255)
                                                                                  ccagccagTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTg               >  1:1405929/1‑100 (MQ=255)
                                                                                   nagccagTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTgc              >  1:187710/2‑100 (MQ=255)
                                                                                    agccagTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTgcg             >  1:823145/1‑100 (MQ=255)
                                                                                     gccagTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTgcgc            >  2:1973243/1‑100 (MQ=255)
                                                                                       cagTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGa          >  2:1710210/1‑100 (MQ=255)
                                                                                         gTTGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACg        <  1:1255758/100‑1 (MQ=255)
                                                                                           tGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGcc      >  2:1242000/1‑100 (MQ=255)
                                                                                           tGTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGcc      <  2:1881810/100‑1 (MQ=255)
                                                                                            gTACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCt     >  2:389007/1‑100 (MQ=255)
                                                                                             tACACCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTg    >  1:2306682/1‑100 (MQ=255)
                                                                                               cGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGAc  >  2:2181215/1‑100 (MQ=255)
                                                                                                |                                                                                                  
CAACAGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGAC  >  NC_000913/899088‑899282

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