Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,914,929 G→A 5.1% L24L (CTC→CTT)  rlmD ← 23S rRNA m(5)U1939 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,914,9290GA 5.1% 256.0 / 2.3 79L24L (CTC→CTT) rlmD23S rRNA m(5)U1939 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base A (2/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

TATCGCTTAACCGGCGTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGC  >  NC_000913/2914784‑2915077
                                                                                                                                                 |                                                                                                                                                    
tATCGCTTAACCGGCGTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAgg                                                                                                                                                   >  7:95421/1‑149 (MQ=255)
  tCGCTTAACCGGCGTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGt                                                                                                                                                  <  2:125437/148‑1 (MQ=255)
    gCTTAACCGGCGTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCg                                                                                                                                                <  7:228782/148‑1 (MQ=255)
     cTTAACCGGCGTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCAAGGTCGtt                                                                                                                                              >  8:219297/1‑149 (MQ=255)
     cTTAACCGGCGTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCAAGGTCGtt                                                                                                                                              >  8:219293/1‑149 (MQ=255)
           cGGCGTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTg                                                                                                                                         >  6:1885/1‑148 (MQ=255)
           cGGCGTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTg                                                                                                                                         >  7:173693/1‑148 (MQ=255)
               gTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACg                                                                                                                                    >  7:327041/1‑149 (MQ=255)
                tACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACgg                                                                                                                                   >  3:199542/1‑149 (MQ=255)
                tACGACTTTAGCGCGGGCATACTGGTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCGGCGGCAATAATCCGGGGTTAATTAGGGTTTTGCCGTTATGTCGCGCCACGCCCTGCCCAAAAGAGGCTAGGTCGTTGCCGGAGACgg                                                                                                                                   >  8:144865/1‑149 (MQ=255)
                    aCTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTAt                                                                                                                               <  5:1885/149‑1 (MQ=255)
                    aCTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGGCGAGGTCGTTGACTGAAACGGTTAt                                                                                                                               <  7:354085/149‑1 (MQ=255)
                     cTTTAGCGCGGGCATACTGTTTTTTACCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATg                                                                                                                              <  7:122572/149‑1 (MQ=255)
                        tAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGag                                                                                                                            <  5:247821/148‑2 (MQ=255)
                        tAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGa                                                                                                                             <  7:17688/147‑1 (MQ=255)
                          gcgcGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCTGCAATAATCCGGGTATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCt                                                                                                                          <  1:38160/148‑1 (MQ=255)
                                    aCTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGcgtcgtc                                                                                                               >  8:251006/1‑149 (MQ=255)
                                     cTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGcgtcgtcg                                                                                                              >  1:348719/1‑149 (MQ=255)
                                        ttttttATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGc                                                                                                           >  8:231107/1‑149 (MQ=255)
                                                cTTCAGTAACAGTACCTTCCGCGTTTTCCTGCGGCATTAATCCGGGGATAAATAGCGTGTTGCCGTTATGTCGCGCCCCGCCCTGACCACAAGAGTCGAGGTCGGTGCCTGCAACGGTTATGCTCTGACGCGGCGTCGGGCGTCGtttt                                                                                                   >  6:264902/1‑149 (MQ=255)
                                                   cagtaacagtaaCTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTg                                                                                                  <  8:95421/147‑1 (MQ=255)
                                                     gtaacagtaacTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCaga                                                                                              <  2:348719/149‑1 (MQ=255)
                                                     gtaacagtaacTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCaga                                                                                              <  8:327041/149‑1 (MQ=255)
                                                      taacagtaacTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCagag                                                                                             >  7:266938/1‑149 (MQ=255)
                                                       aacagtaacTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCagag                                                                                             >  4:268266/1‑148 (MQ=255)
                                                        acagtaacTTCCGCGGTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTa                                                                                           <  3:246008/149‑1 (MQ=255)
                                                         cagtaacTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCagag                                                                                             >  1:252272/1‑146 (MQ=255)
                                                            taacTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGaaa                                                                                                                                      <  6:58110/102‑1 (MQ=255)
                                                            taacTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGaaa                                                                                                                                      >  5:58110/1‑102 (MQ=255)
                                                              acTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAAtt                                                                                      <  4:154118/148‑1 (MQ=255)
                                                                ttCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTgcg                                                                                   >  1:235249/1‑149 (MQ=255)
                                                                ttCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTgcg                                                                                   >  7:263152/1‑149 (MQ=255)
                                                                    gcgTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCAt                                                                               >  4:14016/1‑149 (MQ=255)
                                                                       ttttCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTgg                                                                            <  7:251006/149‑1 (MQ=255)
                                                                         ttCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGc                                                                           <  2:116534/148‑1 (MQ=255)
                                                                           ccTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCAAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGcaa                                                                         <  7:219297/148‑3 (MQ=255)
                                                                           ccTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCAAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCga                                                                         <  7:219293/148‑1 (MQ=255)
                                                                            cTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCgaga                                                                       >  1:49094/1‑149 (MQ=255)
                                                                              gCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACt                                                                     >  6:340156/1‑149 (MQ=255)
                                                                               cGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGAc                                                                      <  6:123489/147‑1 (MQ=255)
                                                                                 gCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTtc                                                                  >  7:58949/1‑149 (MQ=255)
                                                                                 gCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTtc                                                                  <  6:242498/149‑1 (MQ=255)
                                                                                       atCCTGGTATAAATAGCGGTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAAt                                                             <  7:157827/148‑1 (MQ=255)
                                                                                       atCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAAtt                                                            <  5:340156/149‑1 (MQ=255)
                                                                                          cGGGGATACCTAGCGTTTTGCCGTTATCTCGCGCCCCGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCCTTGGCGAGACTTTCTCAATTTaa                                                         <  8:36670/149‑1 (MQ=255)
                                                                                          cGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTaa                                                         <  8:49365/149‑1 (MQ=255)
                                                                                           ggggATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTaa                                                         >  1:302036/1‑148 (MQ=255)
                                                                                           ggggATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAAc                                                        <  8:173693/149‑1 (MQ=255)
                                                                                            gggATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTTACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCAATTGGCGAGACTTTCTCAATTTAACa                                                       <  2:252272/149‑1 (MQ=255)
                                                                                               aTAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAgtg                                                    >  3:239087/1‑149 (MQ=255)
                                                                                                   aTAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATc                                                                                                                           <  4:119243/74‑1 (MQ=255)
                                                                                                   aTAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATc                                                                                                                           >  3:119243/1‑74 (MQ=255)
                                                                                                   aTAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGAc                                                <  2:302036/149‑1 (MQ=255)
                                                                                                       cGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTa                                            <  7:231107/149‑1 (MQ=255)
                                                                                                        gTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTaa                                           <  7:162163/149‑1 (MQ=255)
                                                                                                        gTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGGCGATGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCTTGCGTCGTTTTGCAGAGTAGAATTGAGACATTGTCGATACTTTCTCAATTTAACAGTGTGACCTTaa                                           >  6:159171/1‑149 (MQ=255)
                                                                                                              ccGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTcc                                     <  2:235249/149‑1 (MQ=255)
                                                                                                              ccGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTc                                      <  4:239087/148‑1 (MQ=255)
                                                                                                               cGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGGTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTccc                                    <  4:63316/149‑1 (MQ=255)
                                                                                                                 ttATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCa                                   >  2:115516/1‑148 (MQ=255)
                                                                                                                   aTGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCa                                                                                                >  6:244337/1‑85 (MQ=255)
                                                                                                                   aTGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGGCGTGCGTCGTTTTGCa                                                                                                <  5:244337/85‑1 (MQ=255)
                                                                                                                      tCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACgg                             <  8:58949/149‑1 (MQ=255)
                                                                                                                      tCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACg                              >  7:339046/1‑148 (MQ=255)
                                                                                                                        gcgcCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCACAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCATAATTGTCCCATAACGGaa                           >  7:23745/1‑149 (MQ=255)
                                                                                                                            cACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTcc                       <  7:353921/149‑1 (MQ=255)
                                                                                                                            cACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTc                        >  8:353921/1‑148 (MQ=255)
                                                                                                                              cGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCAt                     <  3:268266/149‑1 (MQ=255)
                                                                                                                              cGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCAt                     <  7:144865/149‑1 (MQ=255)
                                                                                                                               gCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAATCTGTTATGATCTGACGAGTCGTCGTGCGTCGTTTTGCTGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATaaa                               <  6:347892/138‑2 (MQ=255)
                                                                                                                                 ccTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGAc                  >  4:91121/1‑149 (MQ=255)
                                                                                                                                   tGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGCACTCCATGAcc                 >  5:46395/1‑148 (MQ=255)
                                                                                                                                   tGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGAcc                 >  6:296892/1‑148 (MQ=255)
                                                                                                                                       cAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTa            <  1:35957/149‑1 (MQ=255)
                                                                                                                                        aaaaGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACGTTCTGAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTAc           >  4:280566/1‑149 (MQ=255)
                                                                                                                                            gagTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGcc       >  1:29423/1‑149 (MQ=255)
                                                                                                                                               tCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTtct                                                                 <  6:21657/88‑1 (MQ=255)
                                                                                                                                               tCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTtct                                                                 >  5:21657/1‑88 (MQ=255)
                                                                                                                                               tCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTgc    >  8:89063/1‑149 (MQ=255)
                                                                                                                                                 gAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTgcgc  >  3:242098/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TATCGCTTAACCGGCGTACGACTTTAGCGCGGGCATACTGTTTTTTATCTTCAGTAACAGTAACTTCCGCGTTTTCCTGCGGCAATAATCCGGGGATAAATAGCGTTTTGCCGTTATGTCGCGCCACGCCCTGACCAAAAGAGTCGAGGTCGTTGACTGAAACGGTTATGATCTGACGCGTCGTCGTGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGC  >  NC_000913/2914784‑2915077

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