Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,644,260 A→G 5.7% P8P (CCT→CCC rlmN ← dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,644,2600AG 5.7% 221.3 / 2.5 70P8P (CCT→CCCrlmNdual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (29/37);  new base G (2/2);  total (31/39)
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.40e-01

AAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCAC  >  NC_000913/2644113‑2644406
                                                                                                                                                   |                                                                                                                                                  
aaaGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCTGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGTCGGTTGAGATCCAGCATGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAg                                                                                                                                                   <  8:46437/149‑1 (MQ=255)
    ttGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGgtgt                                                                                                                                               <  2:212177/149‑1 (MQ=255)
       tCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCGGGTGTGAc                                                                                                                                            >  7:87152/1‑149 (MQ=255)
       tCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCGGGTGTGAc                                                                                                                                            >  7:152981/1‑149 (MQ=255)
          cagcaATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa                                                                                                                                         <  6:118835/149‑1 (MQ=255)
            gcaATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAAtt                                                                                                                                       <  2:173006/149‑1 (MQ=255)
               aTAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTATCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGtt                                                                                                                                    <  6:244358/149‑1 (MQ=255)
                  gTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGCTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGGTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCa                                                                                                                                  <  2:144984/148‑1 (MQ=255)
                   tGATACATCCAATTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCAGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGa                                                                                                                                <  5:55538/149‑1 (MQ=255)
                    gATACATCCACTTCATCACCTGATCGGCGCGGAAGTGTTTTTCCCCTAAATCTTTAAAAACCTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGAc                                                                                                                               <  4:189628/149‑1 (MQ=255)
                     aTACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACa                                                                                                                              >  6:244234/1‑149 (MQ=255)
                         aTCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGGGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATg                                                                                                                          <  2:110588/149‑1 (MQ=255)
                             acgtcatcaCCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCGGCCGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTAt                                                                                                                      >  3:114238/4‑149 (MQ=255)
                             aCTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTAt                                                                                                                      >  6:357945/1‑149 (MQ=255)
                               ttCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTc                                                                                                                    >  5:256625/1‑149 (MQ=255)
                                 catcaCCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCg                                                                                                                  >  1:54227/1‑149 (MQ=255)
                                  atcaCCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCgg                                                                                                                 >  1:81437/1‑149 (MQ=255)
                                   tcaCCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGc                                                                                                                <  5:357945/149‑1 (MQ=255)
                                      ccTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTc                                                                                                             <  5:244234/149‑1 (MQ=255)
                                       cTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCg                                                                                                            >  3:81820/1‑149 (MQ=255)
                                          aTCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGtta                                                                                                         <  4:131601/149‑1 (MQ=255)
                                                  ggAAGGGTTTTTCACCAAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGt                                                                                                 >  2:128027/1‑149 (MQ=255)
                                                   gAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGtt                                                                                                <  8:213245/149‑1 (MQ=255)
                                                   gAAGGGTTATTCACCTAAATCATAAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTGGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGtt                                                                                                <  4:114238/149‑1 (MQ=255)
                                                    aaGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTa                                                                                               <  2:81437/149‑1 (MQ=255)
                                                    aaGGGTGTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTa                                                                                               <  4:205327/149‑1 (MQ=255)
                                                     aGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTAt                                                                                              <  4:81820/149‑1 (MQ=255)
                                                           tttCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGccc                                                                                         <  1:14799/148‑1 (MQ=255)
                                                            ttCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACTTTATGGCCACt                                                                                       >  6:122398/1‑149 (MQ=255)
                                                             tCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCt                                                                                       <  8:286304/148‑1 (MQ=255)
                                                               aCCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTAAGGCCCCTg                                                                                      >  7:73474/1‑147 (MQ=255)
                                                                 cTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGAc                                                                                                                               <  7:30747/104‑1 (MQ=255)
                                                                 cTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGAc                                                                                                                               >  8:30747/1‑104 (MQ=255)
                                                                  tAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAggg                                                                                 >  6:265684/1‑149 (MQ=255)
                                                                   aaaTCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGc                                                                                                                <  8:218147/117‑1 (MQ=255)
                                                                   aaaTCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGc                                                                                                                >  7:218147/1‑117 (MQ=255)
                                                                      tCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGttatta                                                                                                      >  4:310641/1‑124 (MQ=255)
                                                                      tCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGttatta                                                                                                      <  3:310641/124‑1 (MQ=255)
                                                                            aaaaaCTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCGGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGaaaaaa                                                                        <  8:152981/148‑1 (MQ=255)
                                                                            aaaaaCTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCGGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGaaaaaa                                                                        <  8:87152/148‑1 (MQ=255)
                                                                            aaaaaCTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGaaaa                                                                          >  8:378966/1‑146 (MQ=255)
                                                                             aaaaCTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAg                                                                       >  7:137196/1‑148 (MQ=255)
                                                                               aaCTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGaaa                                                                    <  3:95390/149‑1 (MQ=255)
                                                                                      gCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCg                                                             <  5:154932/149‑1 (MQ=255)
                                                                                        aTCTGCTGACGTTTGCGATCCTGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGAATAATTTTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTTAAAAAAGAAACGCCCCGGt                                                           <  8:137196/149‑1 (MQ=255)
                                                                                        aTCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGt                                                           <  1:128027/149‑1 (MQ=255)
                                                                                             ctgACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCg                                                      <  1:367752/149‑1 (MQ=255)
                                                                                             ctgACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCg                                                      >  2:29934/1‑149 (MQ=255)
                                                                                                  ggTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGgctgc                                                 >  2:203945/1‑149 (MQ=255)
                                                                                                    ttGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTACAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGttat                                                                                                        >  1:330625/1‑92 (MQ=255)
                                                                                                    ttGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTACAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGttat                                                                                                        <  2:330625/92‑1 (MQ=255)
                                                                                                        gaTCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTcc                                           <  2:81471/149‑1 (MQ=255)
                                                                                                          tCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCgg                                         <  6:256625/149‑1 (MQ=255)
                                                                                                             agcagGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCg                                              >  5:327975/1‑141 (MQ=255)
                                                                                                             agcagGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCg                                              <  6:327975/141‑1 (MQ=255)
                                                                                                                 ggTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCt                                  <  6:178517/149‑1 (MQ=255)
                                                                                                                   ttGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGgctgct                                                <  1:343051/133‑1 (MQ=255)
                                                                                                                   ttGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGgctgct                                                >  2:343051/1‑133 (MQ=255)
                                                                                                                   ttGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTg                                 >  7:301044/1‑148 (MQ=255)
                                                                                                                      aTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCAtt                             <  1:203945/149‑1 (MQ=255)
                                                                                                                          ttCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTa                          <  7:378966/148‑1 (MQ=255)
                                                                                                                             cATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACaaa                      >  8:366205/1‑149 (MQ=255)
                                                                                                                               tCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTACGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAAtt                    >  2:245533/1‑149 (MQ=255)
                                                                                                                                      tggtgACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGgcgc             <  7:366205/149‑1 (MQ=255)
                                                                                                                                        gtgACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGGCGCAc           >  8:148296/1‑149 (MQ=255)
                                                                                                                                        gtgACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGGCGCAc           >  8:285336/1‑149 (MQ=255)
                                                                                                                                         tgACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGGCGCACg          >  8:91498/1‑149 (MQ=255)
                                                                                                                                         tgACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGGCGCACg          >  8:394995/1‑149 (MQ=255)
                                                                                                                                             gTTTTCAGTTGTGACTAATTGTTCCGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGATACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTAGAAATTCTGGCGCACGGATg      <  8:73474/149‑1 (MQ=255)
                                                                                                                                                 tCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCAc  >  4:322787/1‑149 (MQ=255)
                                                                                                                                                   aGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGgc               >  3:263007/1‑134 (MQ=255)
                                                                                                                                                   aGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGgc               <  4:263007/134‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAGTTGTCGCAGCAATAGTGATACATCCACTTCATCACCTGATCGGCGCGGAAGGGTTTTTCACCTAAATCTTTAAAAAACTCCCGCATCTGCTGACGGTTGAGATCCAGCAGGTTGATTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAATTGTTCAGACATATGCTATTCCGGCCTCGTTATTACACGTTATGGCCCCTGGAGGGTTGAAAAAAGAAACGCCCCGGTGAGCGGCTGCTCGTCCGGGGGCGCTGCATTGTACAAATTCTGGCGCACGGATGCCAC  >  NC_000913/2644113‑2644406

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