Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 668,354 G→T 5.6% S121* (TCG→TAG)  rlmH ← 23S rRNA m(3)Psi1915 pseudouridine methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913668,3540GT 5.6% 310.5 / 6.6 89S121* (TCG→TAG) rlmH23S rRNA m(3)Psi1915 pseudouridine methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (44/40);  new base T (3/2);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

TGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGCCTGGAATATCGAG  >  NC_000913/668207‑668502
                                                                                                                                                   |                                                                                                                                                    
tGTAGTTTCATCCGCTGCGTTTTCTACTCGAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGa                                                                                                                                                     >  3:793502/1‑149 (MQ=255)
          tCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGc                                                                                                                                           <  1:152557/149‑1 (MQ=255)
          tCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGc                                                                                                                                           >  2:238223/1‑149 (MQ=255)
               tGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAgc                                                                                                                                      >  2:185647/1‑149 (MQ=255)
                    tttCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGAT‑GGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCcgccgc                                                                                                                                >  2:601992/1‑149 (MQ=255)
                       cTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCtt                                                                                                                              <  4:287987/149‑1 (MQ=255)
                       cTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCtt                                                                                                                              <  1:466674/149‑1 (MQ=255)
                        tACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCttt                                                                                                                             <  1:884948/149‑1 (MQ=255)
                           tCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAca                                                                                                                          >  3:823371/1‑149 (MQ=255)
                           tCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAca                                                                                                                          <  4:735305/149‑1 (MQ=255)
                            ccaaccacccTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAcag                                                                                                                         >  2:545002/8‑149 (MQ=255)
                             aaaGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTTCGACCAGCACGCGAACCAGCGGATGGGGGGGGGTAAGCGCCGACAGCGACCAGCTCTGCCCAGCCGCCGCTTTAcagg                                                                                                                        >  4:193768/1‑149 (MQ=255)
                             aaaGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGAc                                                                                                                                                    <  4:211124/121‑1 (MQ=255)
                             aaaGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGAc                                                                                                                                                    >  3:211124/1‑121 (MQ=255)
                             aaaGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAcagg                                                                                                                        >  4:618706/1‑149 (MQ=255)
                             aaaGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAcagg                                                                                                                        >  2:132423/1‑149 (MQ=255)
                             aaaGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAcagg                                                                                                                        >  4:713434/1‑149 (MQ=255)
                                gCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAcaggcag                                                                                                                     <  4:139024/149‑1 (MQ=255)
                                 cTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAcaggcagg                                                                                                                    <  3:713434/149‑1 (MQ=255)
                                  tCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACaggcaggc                                                                                                                   >  1:665772/1‑149 (MQ=255)
                                     cTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGAc                                                                                                                >  1:413826/1‑149 (MQ=255)
                                     cTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGAc                                                                                                                >  1:557382/1‑149 (MQ=255)
                                     cTTATCACTCACGGTGATAAGGATGGTTGGTGGTGAGGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAACGCCGACAGCGACCAGCTCTGCTCAGCCGCCTCTTTACAGGCAGGCGAc                                                                                                                >  4:578712/1‑149 (MQ=255)
                                         tcactcacGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAAcc                                                                                                            <  4:837697/149‑1 (MQ=255)
                                             tcacGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTc                                                                                                        >  1:861139/1‑149 (MQ=255)
                                             tcacGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTc                                                                                                        <  3:532637/149‑1 (MQ=255)
                                             tcacGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTc                                                                                                        >  3:311914/1‑149 (MQ=255)
                                              cacGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCa                                                                                                       <  1:399388/149‑1 (MQ=255)
                                              cacGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTATGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCa                                                                                                       <  2:600625/149‑1 (MQ=255)
                                                 ggTGATAAGGATGGTTGTTGGTGATCCTCCACGCCCGGTACAGACTCTCTGCGACCAGCATGCGAACCAGCGTATGGGTGAGGGTAGGCGCCGACAGCGCCCAGCGCTGCTCAGCCGCCGCTTTATAGGCAGGTGACAACCCGTCAGGc                                                                                                    <  4:644345/149‑1 (MQ=255)
                                                  gTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAcaggcag                                                                                                                     <  4:507849/131‑1 (MQ=255)
                                                  gTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTAcaggcag                                                                                                                     >  3:507849/1‑131 (MQ=255)
                                                  gTGATAAGGATGGTTGGTGGTGATGCTCCACGCCAGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGcc                                                                                                   <  3:335684/149‑1 (MQ=255)
                                                     aTAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGTGACAACCCTTCAGGCCCGc                                                                                                <  1:468699/149‑1 (MQ=255)
                                                     aTAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGc                                                                                                >  1:920243/1‑149 (MQ=255)
                                                      tAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGcc                                                                                               >  1:370404/1‑149 (MQ=255)
                                                         ggATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAAt                                                                                            <  1:185647/149‑1 (MQ=255)
                                                           aTGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCa                                                                                          <  1:704716/149‑1 (MQ=255)
                                                             ggttggtGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGt                                                                                        <  3:553405/149‑1 (MQ=255)
                                                              gttggtGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTa                                                                                       <  1:238223/149‑1 (MQ=255)
                                                              gttggtGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTa                                                                                       <  1:505594/149‑1 (MQ=255)
                                                               ttggtGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAAACAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAg                                                                                      >  1:551083/1‑149 (MQ=255)
                                                                    ggtgATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGa                                                                                 >  2:309214/1‑149 (MQ=255)
                                                                    ggtgATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGa                                                                                 <  2:777978/149‑1 (MQ=255)
                                                                            tCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTc                                                                             <  1:738114/145‑1 (MQ=255)
                                                                            tCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTc                                                                             >  2:738114/1‑145 (MQ=255)
                                                                             ccACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGAc                                                                        >  1:588343/1‑149 (MQ=255)
                                                                              cACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGAcc                                                                       <  1:508937/149‑1 (MQ=255)
                                                                                   ccGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCa                                                                  <  2:370404/149‑1 (MQ=255)
                                                                                    cGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAg                                                                 <  1:456573/149‑1 (MQ=255)
                                                                                      gTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCt                                                               >  1:527432/1‑149 (MQ=255)
                                                                                          aGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCa                                                           >  1:401988/1‑149 (MQ=255)
                                                                                           ggctctctGGGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCTACCAGCTCTGCTCAGCCGCCGCTTTGCAGGCAGCGGACAACCCCTCAGCCCCGCCAATCAGTATACTCACGTCGCGACCATCCAGCTTCCAg                                                          >  4:795545/3‑149 (MQ=255)
                                                                                           gACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAg                                                          >  1:593041/1‑149 (MQ=255)
                                                                                              tctctGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATc                                                                    <  3:171225/136‑1 (MQ=255)
                                                                                              tctctGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATc                                                                    >  4:171225/1‑136 (MQ=255)
                                                                                              tctctGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGt                                                       >  1:440709/1‑149 (MQ=255)
                                                                                              tctctGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGt                                                       >  3:161233/1‑149 (MQ=255)
                                                                                              tctctGCGACCAGCACGCGAACCAGCGGAT‑‑GGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTc                                                     >  1:413349/1‑149 (MQ=255)
                                                                                                tctGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCTACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTc                                                     >  3:397826/1‑149 (MQ=255)
                                                                                                tctGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTc                                                     <  4:366041/149‑1 (MQ=255)
                                                                                                  tGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCa                                                   >  2:850691/1‑149 (MQ=255)
                                                                                                   gCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAg                                                  <  2:627825/149‑1 (MQ=255)
                                                                                                         aGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGc                                            <  1:231652/149‑1 (MQ=255)
                                                                                                          gCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCg                                           >  4:219614/1‑149 (MQ=255)
                                                                                                             cgcgAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCt                                        >  4:318192/1‑149 (MQ=255)
                                                                                                              gcgAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCTACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTa                                       >  1:528461/1‑149 (MQ=255)
                                                                                                               cgAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTaa                                      <  1:859902/149‑1 (MQ=255)
                                                                                                               cgAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTaa                                      <  3:382660/149‑1 (MQ=255)
                                                                                                                      gCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTgcggc                               <  2:588343/149‑1 (MQ=255)
                                                                                                                      gCGGATGCGGGAGGGTAAGCGCCGACAGCTACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTgcggc                               <  3:795545/149‑1 (MQ=255)
                                                                                                                        ggATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGt                             >  4:673684/1‑149 (MQ=255)
                                                                                                                               ggAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAg                      >  2:703244/1‑149 (MQ=255)
                                                                                                                                 aGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAggg                    >  1:545653/1‑149 (MQ=255)
                                                                                                                                 aGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAggg                    <  2:440709/149‑1 (MQ=255)
                                                                                                                                  gggTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGc                   <  1:865711/149‑1 (MQ=255)
                                                                                                                                    gTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCtt                 >  2:219631/1‑149 (MQ=255)
                                                                                                                                    gTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCtt                 >  4:675032/1‑149 (MQ=255)
                                                                                                                                    gTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCtt                 >  1:87330/1‑149 (MQ=255)
                                                                                                                                      aaGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGc               <  1:261614/149‑1 (MQ=255)
                                                                                                                                       aGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGcc              >  2:452537/1‑149 (MQ=255)
                                                                                                                                        gcgcCGACAGCTACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGCCt             <  4:397826/149‑1 (MQ=255)
                                                                                                                                        gcgcCGACAGCGGCCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGCCt             <  1:266546/149‑1 (MQ=255)
                                                                                                                                        gcgcCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGCCt             >  3:802959/1‑149 (MQ=255)
                                                                                                                                           ccGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCg                                           <  3:27053/116‑1 (MQ=255)
                                                                                                                                           ccGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCg                                           >  4:27053/1‑116 (MQ=255)
                                                                                                                                             gACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCtt                 >  4:617881/1‑140 (MQ=255)
                                                                                                                                             gACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCtt                 <  3:617881/140‑1 (MQ=255)
                                                                                                                                              aCAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGCCTGGAata       <  2:414536/149‑1 (MQ=255)
                                                                                                                                                   gACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGCCTGGAATATCGAg  <  2:551083/149‑1 (MQ=255)
                                                                                                                                                   gACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGCCTGGAATATCGAg  <  3:915839/149‑1 (MQ=255)
                                                                                                                                                   gACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGCCTGGAATATCGAg  <  3:675032/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGTAGTTTCATCCGCTGCGTTTTCTACTCAAAGCTCCCTTATCACTCACGGTGATAAGGATGGTTGGTGGTGATGCTCCACGCCCGGTACAGACTCTCTGCGACCAGCACGCGAACCAGCGGATGGGGGAGGGTAAGCGCCGACAGCGACCAGCTCTGCTCAGCCGCCGCTTTACAGGCAGGCGACAACCCTTCAGGCCCGCCAATCAGTAGACTGACGTCGCGACCATCCAGCTTCCAGCGTTCCAGCTCAGCGGCTAACTGCGGCGTATCCCAGGGCTTGCCTGGAATATCGAG  >  NC_000913/668207‑668502

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