Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,409,407 T→A 7.8% S45T (TCT→ACT)  rlmB → 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,4070TA 7.8% 452.1 / 19.3 138S45T (TCT→ACT) rlmB23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (61/66);  new base A (6/5);  total (67/71)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 141 total aligned reads displayed.

AACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCA  >  NC_000913/4409262‑4409555
                                                                                                                                                 |                                                                                                                                                     
 acGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑a                                                                                                                                                  <  1:367858/149‑1 (MQ=255)
   gAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑Agg                                                                                                                                                <  3:463785/149‑1 (MQ=255)
        cATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                           <  4:508230/149‑1 (MQ=255)
         aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                          <  2:934447/149‑1 (MQ=255)
         aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                          <  3:751343/149‑1 (MQ=255)
         aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGTga                                                                                                                                          <  4:215455/149‑2 (MQ=255)
            aaTGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTa                                                                                                                                       >  4:867340/1‑149 (MQ=255)
             aTGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGTgaatca                                                                                                                                      <  2:910155/149‑6 (MQ=255)
               gAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATcc                                                                                                                                    >  3:896304/1‑149 (MQ=255)
                  cGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGt                                                                                                                                 >  3:884640/1‑149 (MQ=255)
                    aaaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTg                                                                                                                               >  1:375871/1‑149 (MQ=255)
                             aCGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGc                                                                                                                      <  4:662069/149‑1 (MQ=255)
                               ggCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCa                                                                                                                    >  2:39898/1‑149 (MQ=255)
                                  aTCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAata                                                                                                                 >  4:744804/1‑149 (MQ=255)
                                   tCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAatat                                                                                                                >  2:342756/1‑149 (MQ=255)
                                      aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGt                                                                                                                                            <  4:674744/118‑1 (MQ=255)
                                      aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGt                                                                                                                                            >  3:674744/1‑118 (MQ=255)
                                      aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGACTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtctc                                                                                                             >  1:644647/1‑149 (MQ=255)
                                       cGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcg                                                                                                            >  2:669246/1‑149 (MQ=255)
                                       cGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcg                                                                                                            >  1:246712/1‑149 (MQ=255)
                                         cAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                                                                                          <  3:435044/149‑1 (MQ=255)
                                              cAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaa                                                                                                     <  3:213354/149‑1 (MQ=255)
                                                ggCCCTGCTGGTGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTAAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaaa                                                                                                   <  1:788825/149‑1 (MQ=255)
                                                ggCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaaa                                                                                                   <  4:890183/149‑1 (MQ=255)
                                                  cccTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGc                                                                                                 <  3:759215/149‑1 (MQ=255)
                                                        tgGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                                                                                           >  2:125878/1‑149 (MQ=255)
                                                        tgGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                                                                                           <  4:896304/149‑1 (MQ=255)
                                                         ggAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTg                                                                                          <  1:342756/149‑1 (MQ=255)
                                                             cgcgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGt                                                                                      >  4:143740/1‑149 (MQ=255)
                                                                  cccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                                 >  4:399783/1‑149 (MQ=255)
                                                                  cccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                                 <  1:86684/149‑1 (MQ=255)
                                                                   ccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGa                                                                                               <  2:583969/134‑1 (MQ=255)
                                                                   ccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCa                                                                                <  2:375871/149‑1 (MQ=255)
                                                                    cTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAg                                                                               >  4:503417/1‑149 (MQ=255)
                                                                     tGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                              <  4:370627/149‑1 (MQ=255)
                                                                     tGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                              >  4:467501/1‑149 (MQ=255)
                                                                     tGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATCCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                              >  1:604838/1‑149 (MQ=255)
                                                                         cGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCAt                                                                          <  4:120855/149‑1 (MQ=255)
                                                                           ttttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTa                                                                        >  4:485042/1‑149 (MQ=255)
                                                                            tttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTAt                                                                       >  2:397710/1‑149 (MQ=255)
                                                                                aGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGcc                                                                   >  1:357557/1‑149 (MQ=255)
                                                                                 ggAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGccc                                                                  >  3:865494/1‑149 (MQ=255)
                                                                                   aaGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCcgc                                                                <  4:865494/149‑1 (MQ=255)
                                                                                        tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACTGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                           >  3:449006/1‑149 (MQ=255)
                                                                                        tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTgg                                                                                                                              <  1:684929/82‑1 (MQ=255)
                                                                                        tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTgg                                                                                                                              >  2:684929/1‑82 (MQ=255)
                                                                                        tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAg                                                                                                  >  2:718054/1‑110 (MQ=255)
                                                                                        tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTgca                                                           >  3:561716/1‑147 (MQ=255)
                                                                                        tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                           <  2:501062/149‑1 (MQ=255)
                                                                                        tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                           >  4:647974/1‑149 (MQ=255)
                                                                                           atgttAAAGGGCCGGGAAGATAAACGTCTGTTACCGCCGAAACGCCGCCTTGAGTCTC‑AGAGCGTGGTTATCGAGTTGGCAAACCGCCAATATCTCGACGACAAACGCGACGGTGCCGTGCCTCAGGGCATTATCGTCCGCGTGAAAcc                                                        <  2:70352/146‑1 (MQ=255)
                                                                                            ttttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCa                                                       <  1:226674/149‑1 (MQ=255)
                                                                                             tttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAg                                                      >  3:477322/1‑149 (MQ=255)
                                                                                               tAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTTAAGCCAGGa                                                    <  1:716499/149‑1 (MQ=255)
                                                                                                aaaaGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGAc                                                   <  1:8059/149‑1 (MQ=255)
                                                                                                 aaaGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACg                                                  >  1:241063/1‑149 (MQ=255)
                                                                                                  aaGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGt                                                 <  1:39898/149‑1 (MQ=255)
                                                                                                   aGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                                >  2:735961/1‑149 (MQ=255)
                                                                                                       cGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTa                                            <  1:785477/149‑1 (MQ=255)
                                                                                                          gAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATc                                                                      <  1:573429/120‑1 (MQ=255)
                                                                                                            aGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAgg                                       <  2:162783/149‑1 (MQ=255)
                                                                                                             gATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGa                                      >  3:647163/1‑149 (MQ=255)
                                                                                                                aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaaa                                   <  3:143740/149‑1 (MQ=255)
                                                                                                                aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaaa                                   <  2:155164/149‑1 (MQ=255)
                                                                                                                   cGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGa                                <  4:299688/149‑1 (MQ=255)
                                                                                                                    gTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGAt                               <  3:399783/149‑1 (MQ=255)
                                                                                                                       tGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTg                            >  1:138082/1‑149 (MQ=255)
                                                                                                                       tGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTg                            >  2:120653/1‑149 (MQ=255)
                                                                                                                       tGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTg                            >  4:178069/1‑149 (MQ=255)
                                                                                                                        gTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGc                           <  3:612338/149‑1 (MQ=255)
                                                                                                                         ttACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                                 <  1:886686/94‑1 (MQ=255)
                                                                                                                         ttACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                           >  2:76213/1‑116 (MQ=255)
                                                                                                                         ttACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGcc                          >  1:16607/1‑149 (MQ=255)
                                                                                                                          tACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCg                         <  2:407774/149‑1 (MQ=255)
                                                                                                                             cgttGATTCCCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAGGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAATCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                      <  1:41419/146‑1 (MQ=255)
                                                                                                                             cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                              >  3:229958/1‑93 (MQ=255)
                                                                                                                             cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                      <  1:154966/149‑1 (MQ=255)
                                                                                                                             cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                      >  3:97252/1‑149 (MQ=255)
                                                                                                                             cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTGATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                              <  4:229958/93‑1 (MQ=255)
                                                                                                                              gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCa                                               >  2:645204/1‑123 (MQ=255)
                                                                                                                              gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCa                                               <  1:645204/123‑1 (MQ=255)
                                                                                                                              gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                     >  4:650125/1‑149 (MQ=255)
                                                                                                                               cTGATCCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCt                    <  2:604838/149‑1 (MQ=255)
                                                                                                                                  aTTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAt                 <  1:747001/149‑1 (MQ=255)
                                                                                                                                   accACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtc                >  2:841194/3‑149 (MQ=255)
                                                                                                                                    ccACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCCGTATCAGGAAAACGATCTGCCGGATCTGAtcg               >  3:580144/2‑149 (MQ=255)
                                                                                                                                    accacgccTTGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACTGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtc                >  3:852459/7‑149 (MQ=255)
                                                                                                                                     cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgc              >  1:260031/1‑149 (MQ=255)
                                                                                                                                      aCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                     >  2:261332/1‑141 (MQ=255)
                                                                                                                                      aCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                     <  1:261332/141‑1 (MQ=255)
                                                                                                                                       cGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgctt            <  2:357557/149‑1 (MQ=255)
                                                                                                                                       cGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgctt            >  2:724570/1‑149 (MQ=255)
                                                                                                                                        gCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttc           <  2:241063/149‑1 (MQ=255)
                                                                                                                                        gCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttc           <  1:735961/149‑1 (MQ=255)
                                                                                                                                           cctGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttcgc         >  2:287392/3‑149 (MQ=255)
                                                                                                                                            ttGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTc       >  4:169122/1‑149 (MQ=255)
                                                                                                                                            ttGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTc       >  2:880769/1‑149 (MQ=255)
                                                                                                                                              gAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGa     <  3:867340/149‑1 (MQ=255)
                                                                                                                                                 tctc‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATTTGCCGGATCTGATCGCTTCGCTCGATCa  >  2:110095/1‑149 (MQ=255)
                                                                                                                                                 tctc‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCa  <  4:511056/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCA  >  NC_000913/4409262‑4409555

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