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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,4070TA 7.2% 137.8 / 2.5 67S45T (TCT→ACT) rlmB23S rRNA 2'‑O‑ribose G2251 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (25/34);  new base A (2/3);  total (28/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 = 8.72e-01

GAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTC  >  NC_000913/4409265‑4409546
                                                                                                                                              |                                                                                                                                            
gAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGt                                                                                                                                              >  1:580681/1‑143 (MQ=255)
gAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAAAGTTATAAAGAAGTCTTAATATTAAAAGGCCGTGAAGATAAACGTCTGTTACAGCTGATTCACGCCCTTGAGt                                                                                                                                              >  1:24779/1‑143 (MQ=255)
   tACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGCTTTCAGGAAGTCCTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGCGcctc                                                                                                                                           >  2:220223/1‑143 (MQ=255)
         aaTGAGCGAAATGATTTACGGCATCCACGGAGTGCAGGACGGGCGTTTGCGTGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTGAGACTCAAGGGCg                                                                                                                                    <  2:402008/143‑1 (MQ=255)
            gAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTgg                                                                                                                                 <  2:124337/143‑1 (MQ=255)
            gAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGGCTGTTACCGCTGATTCACGCCCTGAGACTCAAGGGCGTga                                                                                                                                 <  1:287972/143‑2 (MQ=255)
            gAGCGAAATGATTTACGGCAGCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTgg                                                                                                                                 <  1:511856/143‑1 (MQ=255)
             aGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCGAGAGACACAGGGCGTGGt                                                                                                                                >  2:298317/1‑143 (MQ=255)
              gCGAAATTATTTCCTGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATGAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGtt                                                                                                                               <  1:554475/143‑1 (MQ=255)
              gCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGtt                                                                                                                               >  2:261807/1‑143 (MQ=255)
                gAAATGATTTACGGCATCCACGCAGTGCAGGTCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCAAAATTTTGGTAGGCCATGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTAt                                                                                                                             <  1:101601/143‑1 (MQ=255)
                 aaaTGATTTACGGCATCCACGCAGTGCATGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGTCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATc                                                                                                                            <  2:226632/143‑1 (MQ=255)
                 aaaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATc                                                                                                                            <  2:570163/143‑1 (MQ=255)
                      aTTTACGGCACCCACGCAGTGCAGGCCCTGCTGGAGCGCGGCTCTGAACGTGTGCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCGGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGtt                                                                                                                       <  2:234958/143‑1 (MQ=255)
                       tttACGGCATCCACGCAGTGCAGGCCCTACAACAAAGCACCCATGCTTGATCCCAGGAAGTCTTTATTCTAAAAGGCCGTGAAGATAAAAGTCTGTTAACGCTGATTCACGCACTTGAATCTCAGGGCGTGGTTATCCAGTTg                                                                                                                      >  1:234681/1‑143 (MQ=255)
                        ttACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTgg                                                                                                                     <  1:570353/143‑1 (MQ=255)
                         tACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCCTTATTTTAAAAGGCCGTAAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGc                                                                                                                    >  1:447569/1‑143 (MQ=255)
                           cGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAATCTCAGGGCGTGGTAAACCAGTCGGCaa                                                                                                                  >  2:476808/1‑143 (MQ=255)
                            ggCATCCACACAGTGCAGGCCCTGCAGGAGCGCGCCCCTGAACGTTTTCAGGAAGACTTTATCTTAAAAGGCCGTGAAGATAAACGTCTGCTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGCTATCCAGTTGGCaaa                                                                                                                 >  2:254016/1‑143 (MQ=255)
                             gCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCATGGCGTGGTTATCCAGTTGGCAAAc                                                                                                                <  1:492153/143‑1 (MQ=255)
                              cATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTCGCAAAcc                                                                                                               >  2:424697/1‑143 (MQ=255)
                               aTCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCCGGGCGTGGTTATCCAGTTGGCAAACCg                                                                                                              >  2:519706/1‑143 (MQ=255)
                                  cACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCGGAACGTTTTCAGGAAGTCTTTATGTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGCGCGTGGTTATTCAGTTGGCcaatcagca                                                                                                           >  1:558318/1‑137 (MQ=255)
                                        gTGCTGGCCCTGCTGGGGGGCTCCACTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAtct                                                                                                     <  2:334845/143‑1 (MQ=255)
                                                         gcgcgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGACAAACCGCCAATATCTCGACGAGAAAAGCGACg                                                                                    >  1:455819/1‑143 (MQ=255)
                                                         gcgcgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCAAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACg                                                                                    >  1:520313/1‑143 (MQ=255)
                                                           gctccccTGCACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCTCGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                                                                                  <  2:281514/140‑1 (MQ=255)
                                                                 ctcaaCGTTTTCAGTAAGTGTTTATTTTAAAAGGGCGGGAAGAAAACCGGCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTg                                                                            <  1:227225/140‑1 (MQ=255)
                                                                 ccgacagTTTTCAGGAACAGTTGATTCTAAAAGGCCGTGAAGAGAAACGTCTGTTACCGCTGATTCACGCCTGTAAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTg                                                                            <  1:187213/137‑1 (MQ=255)
                                                                  tGATCGTTTTCAGTATGTCTTTGTTTTAAACGGTCGTGAGAGTAAACGTCTGTTACCGCTGATTCGCGCCCTTGAGTCTTAGGGCGTGGTTATCCAGTTGTCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGc                                                                           <  2:560248/143‑1 (MQ=255)
                                                                        ttttCAGGAAGTCTTTATTTTAAAAGGCCGTGAATATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                     <  2:422897/143‑1 (MQ=255)
                                                                        gtcgcAGGGAGTGTTTGTGTGAAAGGGGCGTTAAGATAAACGTCTTGTACCGCTGATTCACGCCCTTGATGCTCAGGGCGTGGTTATCCAGTCGGCTAAGTGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                     <  2:247946/139‑1 (MQ=255)
                                                                         tttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCAAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAggg                                                                    >  1:520399/1‑143 (MQ=255)
                                                                                aaGTCTTTATTTTAAAAGTCCGTGAAGATAATCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGCGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATc                                                             <  1:293775/143‑1 (MQ=255)
                                                                                     tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGAATCAGAGCAGTATCGccag                                                        >  1:349624/1‑141 (MQ=255)
                                                                                      ttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGTGTCTCAGGGCGTGGTTATCCAGTTGGCAAACAACAAATATATCGAAGAGAAAAACGAAGGTGCCGTGCATCAAGGCATTATCGccccc                                                       >  1:144030/1‑141 (MQ=255)
                                                                                       gatgttAAAAGGCCGTGAAGATAAAAGTCTGTAACCGCTGATCCACGACATTGAGACGCAGAGCGTGGACATCACGTTGTCAAAATGCAGACAACTGGAGCAGAAACGAGACGGTGCGGAGCCTAAGGGCATAATCGccccca                                                      >  2:371752/5‑140 (MQ=255)
                                                                                          tttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGCAAAGCGACGGGGCCGTGCAGCAGGGCATTCTCGCCCGCGTgc                                                   >  2:239960/1‑142 (MQ=255)
                                                                                          tttAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACACCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTACACGAGAAAAGCGACGGTGCCGTGCATCAGGGAAGTATCGCACGCGTGa                                                   >  2:552423/1‑143 (MQ=255)
                                                                                            tAAAAGGCCGTGAAGATATTCGTAGGCAACCCCTGGTTTCCGCACAAAAAAGTCAGGGCGTGGTTGTCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAg                                                 <  1:253129/143‑1 (MQ=255)
                                                                                              aaaGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGcc                                               <  2:382545/143‑1 (MQ=255)
                                                                                                aggtcttgatgaTAAACGTCTGTTACCGCTGATTCACGCCCGTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAg                                             <  1:560411/137‑1 (MQ=255)
                                                                                                aGGCCGTGAATATATTCGTCTGTTATCGTTGTTTCACGTCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAg                                             <  1:467977/143‑1 (MQ=255)
                                                                                                aGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAg                                             <  1:379576/143‑1 (MQ=255)
                                                                                                     gTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAAGGCATTGTCGCCACCGTTCAGCCAAGACga                                        >  2:334262/1‑142 (MQ=255)
                                                                                                       gAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCa                                      <  2:447569/143‑1 (MQ=255)
                                                                                                         aGATAAACGTCTGGTGCCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATAGCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGt                                    <  2:395078/143‑1 (MQ=255)
                                                                                                           aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTAt                                  >  2:550191/1‑143 (MQ=255)
                                                                                                            tAAACGTCTGTTACCGCTGAGTCACGCCCTTGGGGCGCAGGGCGTGGGTATCCAGTTGGCAAACCGCCAATATCGCGACGAGAAAAGCGACGGTGTCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATc                                 <  1:553182/143‑1 (MQ=255)
                                                                                                            tAAACGTCTGTTACCGATGATTGACGCCCTTGAGTCTCAGGGCGTGGATTTCCGGGTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCATGACGTCAGTATc                                 >  1:127129/1‑143 (MQ=255)
                                                                                                              accGTCTGTTACCGCTGCTTCACGCCCTTGAGTTTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAg                               <  1:261807/141‑1 (MQ=255)
                                                                                                                cGTCTGGTAGCGCTGATGCACGCGCTTGATTGTCAGGGCGTGGTGATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGGCCGCGTGAAGCCAGGACGTCAGTATCAGGa                             <  2:435743/143‑1 (MQ=255)
                                                                                                                  tCTGTTACCGCTGATTCACGCTTTTTAGTCTCAGGGCGTGTATAACCAGTTGGCAAACCGCCAATATCACGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaa                           <  1:174670/143‑1 (MQ=255)
                                                                                                                   cTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGTCGTGGTTATCCAGTTGGCAATCCGCTAATATCTCGACGAGAAAAGCGACGGTGCCGTGCAGCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaaa                          <  1:461294/143‑1 (MQ=255)
                                                                                                                      ttACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATCTCTCGCCGAGCAAAGCGACGGTGCCGTGCATCAGATCATTATTGCACACGTAAAAAAAAGACGTCAGTATCAGGAAAACGa                       >  1:166381/1‑143 (MQ=255)
                                                                                                                      ttACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCTTGAAGCCAGGACGTCAGTATCAGGAAAACGa                       >  1:326486/1‑143 (MQ=255)
                                                                                                                           gCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGCGAAAAGAGACGGTGCTGTGCATCAGGGCATTCTCGCCCGCGTGCAGCCAGCGCGTCAGTATCAGGAATCCGCTCTGc                  >  2:220240/1‑143 (MQ=255)
                                                                                                                           gCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGCGAAGCCAGGACGTCAGTATCAGGAAAACGAACTGc                  >  2:494485/1‑143 (MQ=255)
                                                                                                                           gCTGATTCACGCCCTTGAGCCTCAGGGCGTGGTTATCCAGTTGGCACACCGCCAATATTGCGACGCGCAAAGCGGCGCTGCCGTGCACCAGGTCATCATCGCACGCGTGAAGCCAGGACGTCAGTATACGGCCAACGAACTGc                  >  2:132955/1‑143 (MQ=255)
                                                                                                                             tGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCg                <  1:301041/143‑1 (MQ=255)
                                                                                                                                ttaaCGCCCTCGGGTCTCAGAGCGTGGTTATCCAGTTGGCAAACCGCCACTATCTCGACGAGCAACGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt             <  2:129437/140‑1 (MQ=255)
                                                                                                                                ttCACGGCCTTGAGTCTCAGGTTGTGGTTATCTAGTTGACAAACCGCCAAGATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCACGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt             <  1:197608/143‑1 (MQ=255)
                                                                                                                                ttCACGCCCTTGAGTCTCAGGGCGTGGTTATTCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt             <  1:220223/143‑1 (MQ=255)
                                                                                                                                ttCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt             <  1:84178/143‑1 (MQ=255)
                                                                                                                                 tCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCATGACGTCAGTATCAGGAAAACGATCTGCCGGATc            <  2:127129/143‑1 (MQ=255)
                                                                                                                                   aCGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCACGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTg          >  1:328590/1‑143 (MQ=255)
                                                                                                                                   aCGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACAGTGCCGTGCATCAGAGCATTAACGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTg          >  1:485965/1‑143 (MQ=255)
                                                                                                                                    cGCCCTCGAGTCTCACGGCGTCGCTCTCCAGTTGGCAAACCTCCAATATCTCGACGAGAAAAGTGACGGTGTCGTGCATCAGGGCATTACCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGa         <  1:18147/143‑1 (MQ=255)
                                                                                                                                         ttGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCGAGGCGGCCAGCACCAGGAAAACGATCTGCCCGCTCTGCtcgct    >  2:561455/1‑143 (MQ=255)
                                                                                                                                           gAGTCTCAGGGCGTTGTTATTCAGTTGGCAAACTGCTAATATTTTGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttc  <  1:242416/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTC  >  NC_000913/4409265‑4409546

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