Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,906,661 A→C F134V (TTT→GTT)  rlmA ← 23S rRNA m(1)G745 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,906,6610AC98.3% 617.7 / ‑2.2 175F134V (TTT→GTT) rlmA23S rRNA m(1)G745 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base C (87/84);  total (88/87)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

CGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGC  >  NC_000913/1906564‑1906757
                                                                                                 |                                                                                                
cGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAg                                                                                                <  2:5093037/100‑1 (MQ=255)
 gTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAgn                                                                                               <  1:1534545/100‑2 (MQ=255)
 gTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGt                                                                                               <  2:224715/100‑1 (MQ=255)
 gTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGt                                                                                               <  2:103091/100‑1 (MQ=255)
  tGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTg                                                                                              >  1:2025055/1‑100 (MQ=255)
    cTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGAc                                                                                            <  2:4127463/100‑1 (MQ=255)
     tAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGAcc                                                                                           <  2:1373446/100‑1 (MQ=255)
      aaTTCTTCTGCTTTACACTGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCt                                                                                          >  2:2331226/1‑100 (MQ=255)
      aaTTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCt                                                                                          >  1:2173140/1‑100 (MQ=255)
      aaTTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCt                                                                                          >  2:1561236/1‑100 (MQ=255)
        ttcttcTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGc                                                                                        <  2:29863/100‑1 (MQ=255)
        ttcttcTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGc                                                                                        >  1:4082994/1‑100 (MQ=255)
         tcttctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCg                                                                                       <  2:1458386/100‑1 (MQ=255)
         tcttctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCg                                                                                       >  1:4837905/1‑100 (MQ=255)
           ttctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGa                                                                                     >  1:5471581/1‑100 (MQ=255)
            tctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGAt                                                                                    <  1:3086610/100‑1 (MQ=255)
             ctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGACGCGACACAAACAGTGACCGGCGGaga                                                                                   >  1:3422164/1‑98 (MQ=255)
               gCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGc                                                                                 >  2:709984/1‑100 (MQ=255)
               gCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGc                                                                                 >  1:5013026/1‑100 (MQ=255)
               gCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGc                                                                                 >  2:637782/1‑100 (MQ=255)
                  ttACACGGCGCGTAAATGCGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGtt                                                                              >  1:678725/1‑100 (MQ=255)
                  ttACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGtt                                                                              >  2:1503783/1‑100 (MQ=255)
                    acacGGCGCGTAAATACGTATTATGGCGTCCATACCGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTc                                                                            <  1:1685486/100‑1 (MQ=255)
                     cacGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTcg                                                                           <  1:5364360/100‑1 (MQ=255)
                       cGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTcgcc                                                                         >  2:190693/1‑100 (MQ=255)
                          cgcgTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTcgccgcc                                                                      >  2:5356422/1‑100 (MQ=255)
                            cgTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTcgccgccgc                                                                    >  1:1873211/1‑100 (MQ=255)
                             gTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGGCACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCt                                                                   >  1:5150964/1‑100 (MQ=255)
                             gTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCt                                                                   >  2:4005929/1‑100 (MQ=255)
                                aaTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCtttt                                                                >  1:2810120/1‑100 (MQ=255)
                                 aTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTa                                                               >  1:4960391/1‑100 (MQ=255)
                                  tACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTAt                                                              <  2:2639268/100‑1 (MQ=255)
                                  tACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTAt                                                              <  2:5028266/100‑1 (MQ=255)
                                   aCGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATc                                                             <  1:3010915/100‑1 (MQ=255)
                                    cGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCg                                                            >  1:284907/1‑100 (MQ=255)
                                     gTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGc                                                           <  1:95489/100‑1 (MQ=255)
                                     gTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGc                                                           <  1:2398962/100‑1 (MQ=255)
                                       attatGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAACGAAAAAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTa                                                         <  2:3845537/100‑1 (MQ=255)
                                         tatGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTAcc                                                       <  1:1707681/100‑1 (MQ=255)
                                         tatGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTAcc                                                       <  1:1557491/100‑1 (MQ=255)
                                          atGGTGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACCAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCt                                                      <  1:4988811/100‑1 (MQ=255)
                                          atGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCt                                                      >  1:5517463/1‑100 (MQ=255)
                                            ggCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTc                                                    >  2:3830116/1‑100 (MQ=255)
                                            ggCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTc                                                    <  2:4881582/100‑1 (MQ=255)
                                              cGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGa                                                  >  1:926380/1‑100 (MQ=255)
                                               gTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGaa                                                 >  2:3447764/1‑100 (MQ=255)
                                               gTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGaa                                                 <  1:967613/100‑1 (MQ=255)
                                               gTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGaa                                                 <  1:192563/100‑1 (MQ=255)
                                               gTCAATACTGGTATCGAAAAACGGCAAACGGTCGCTGGAAGCGACACAAAAAGTGACCTGAGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGaa                                                 <  2:831121/100‑1 (MQ=255)
                                                tCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGaaa                                                <  2:3231786/100‑1 (MQ=255)
                                                tCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTACCGCTACCTTCGaaa                                                <  2:2892795/100‑1 (MQ=255)
                                                  cATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACa                                              <  1:1701812/100‑1 (MQ=255)
                                                      cTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACCTCCa                                          >  2:4062250/1‑100 (MQ=255)
                                                      cTGGCATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCa                                          <  1:5376681/100‑1 (MQ=255)
                                                           aTCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCa                                     >  2:207051/1‑100 (MQ=255)
                                                             cGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCaaa                                   >  2:2696209/1‑100 (MQ=255)
                                                             cGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCaaa                                   >  1:110474/1‑100 (MQ=255)
                                                             cGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCaaa                                   >  1:110138/1‑100 (MQ=255)
                                                              ggAAAACGGCAAACGGTGGCTGGAAGCTACACAAAGAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAAc                                  >  1:5387801/1‑100 (MQ=255)
                                                              ggAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAAc                                  >  2:1272482/1‑100 (MQ=255)
                                                              ggAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAAc                                  >  1:3223154/1‑100 (MQ=255)
                                                              ggAAAACGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACACCCAGACCAAAc                                  >  2:445480/1‑100 (MQ=255)
                                                                aaaaCGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGGTAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgt                                >  2:2846805/1‑100 (MQ=255)
                                                                aaaaCGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgt                                >  1:4337132/1‑100 (MQ=255)
                                                                 aaaCGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgtg                               >  2:3713156/1‑100 (MQ=255)
                                                                 aaaCGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgtg                               >  2:3359470/1‑100 (MQ=255)
                                                                 aaaCGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgtg                               >  1:5055879/1‑100 (MQ=255)
                                                                  aaCGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgtgg                              <  1:3658424/100‑1 (MQ=255)
                                                                   aCGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgtggt                             <  2:4201441/100‑1 (MQ=255)
                                                                   aCGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgtggt                             >  1:675876/1‑100 (MQ=255)
                                                                   aCGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgtggt                             <  1:433486/100‑1 (MQ=255)
                                                                    cGGCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgtggtg                            >  2:951554/1‑100 (MQ=255)
                                                                    cGGCAAACGGTGGCTGGAAGCGACACAAACAATGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACgtggtg                            <  2:4478263/100‑1 (MQ=255)
                                                                     ggCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGa                           <  2:1815669/100‑1 (MQ=255)
                                                                     ggCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGa                           <  1:2415507/100‑1 (MQ=255)
                                                                      gCAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGAt                          >  2:3051615/1‑100 (MQ=255)
                                                                       cAAACGGTTGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTGTATCGCTACCTTCGACACATCCAGACCAAACGTGGTGAtt                         <  2:1194746/100‑1 (MQ=255)
                                                                       cAAACGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGAtt                         >  2:64419/1‑100 (MQ=255)
                                                                          aCGGTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCg                      <  2:1692955/100‑1 (MQ=255)
                                                                            ggTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCggg                    >  1:84629/1‑100 (MQ=255)
                                                                             gTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGc                   <  1:5303370/100‑1 (MQ=255)
                                                                             gTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGc                   <  2:574848/100‑1 (MQ=255)
                                                                             gTGGCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGc                   <  2:4568946/100‑1 (MQ=255)
                                                                                gCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAAc                <  2:1987349/100‑1 (MQ=255)
                                                                                gCTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAAc                <  2:4741753/100‑1 (MQ=255)
                                                                                 cTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACg               >  2:1291619/1‑100 (MQ=255)
                                                                                 cTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACg               >  1:725842/1‑100 (MQ=255)
                                                                                 cTGGAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACg               <  1:1516029/100‑1 (MQ=255)
                                                                                   ggAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCa             <  1:2269345/100‑1 (MQ=255)
                                                                                    gAAGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCAt            >  1:154542/1‑100 (MQ=255)
                                                                                     aaGCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATc           <  2:3578320/100‑1 (MQ=255)
                                                                                       gCGACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCgg         >  1:176001/1‑100 (MQ=255)
                                                                                         gACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCa       <  1:3995231/100‑1 (MQ=255)
                                                                                         gACACAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCa       <  1:3995392/100‑1 (MQ=255)
                                                                                          acacAAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGTTACCTTCGAAACATCCAGTCCAAACGTGGTGATTTCGGGCAACGCATCGGCaa      <  1:5173884/100‑1 (MQ=255)
                                                                                           cacaAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCaaa     >  1:1564699/1‑100 (MQ=255)
                                                                                           cacaAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCaaa     >  1:4151608/1‑100 (MQ=255)
                                                                                            acaAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCaaan    <  1:4449400/100‑2 (MQ=255)
                                                                                             caAACAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATg   >  2:3489287/1‑100 (MQ=255)
                                                                                              aaaaaGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGc  >  2:2326223/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGC  >  NC_000913/1906564‑1906757

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