Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,037 T→C 15.6% intergenic (+31/‑96) rlmB → / → yjfI 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,0370TC15.6% 258.9 / 40.4 124intergenic (+31/‑96)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (50/53);  new base C (8/11);  total (60/64)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-01

Only 100 of 126 total aligned reads displayed.

GTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATG  >  NC_000913/4409938‑4410135
                                                                                                   |                                                                                                  
gTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGAtt                                                                                                    >  1:1421490/1‑100 (MQ=255)
 tttCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGAttt                                                                                                   <  2:630296/100‑1 (MQ=255)
   tcttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTgg                                                                                                 >  2:1646831/1‑100 (MQ=255)
     ttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                               <  2:235838/100‑1 (MQ=255)
      tcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATg                                                                                              >  1:905072/1‑100 (MQ=255)
        ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATggc                                                                                            <  2:1687942/100‑2 (MQ=255)
        ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtt                                                                                           >  1:268687/1‑100 (MQ=255)
          tGAACGTGTCGGTTGCGACCGGTATTTGCTTATATGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGctt                                                                                          <  1:2408755/100‑1 (MQ=255)
            aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttt                                                                                         >  2:2433445/1‑99 (MQ=255)
             aCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATACAGATTTGGATGGtttttt                                                                                       >  1:270441/1‑100 (MQ=255)
              cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                      <  2:852548/100‑1 (MQ=255)
              cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                      >  2:1739927/1‑100 (MQ=255)
               gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                     <  2:1028889/100‑4 (MQ=37)
               gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAACTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                     <  2:727954/100‑4 (MQ=37)
                tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                    <  2:2248954/100‑2 (MQ=255)
                tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTAAGAATAATGAAAAACCATCCAAATCTGGATGGcttttc                                                                                       <  1:233140/97‑2 (MQ=37)
                  tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                  <  1:2352233/100‑4 (MQ=255)
                  tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                  <  2:2099963/100‑4 (MQ=255)
                  tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgt                                                                                   >  1:2180511/1‑97 (MQ=255)
                   cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                 >  2:2114790/1‑96 (MQ=255)
                   cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                 <  1:2561538/100‑5 (MQ=255)
                   cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataatt                                                                                 <  1:1862933/100‑8 (MQ=25)
                    ggTTGCGACCGGAATTTGCTTATTTGTAGCTTTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataattc                                                                                <  1:1730953/100‑9 (MQ=21)
                    ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                               >  1:110390/1‑100 (MQ=255)
                    ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                               >  2:1489052/1‑100 (MQ=255)
                    ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataattc                                                                                <  2:1946639/100‑9 (MQ=25)
                    ggTTGCGACCGGAATTTGCTTAGTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataattc                                                                                <  1:1774811/100‑9 (MQ=25)
                     gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtat                                                                              <  1:879542/100‑1 (MQ=255)
                     gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtat                                                                              <  1:2088872/100‑1 (MQ=255)
                      ttGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtata                                                                             <  2:1492072/100‑1 (MQ=255)
                       tGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtcttt                                                                               >  1:1973761/1‑92 (MQ=255)
                       tGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAg                                                                            >  1:1628861/1‑100 (MQ=255)
                        gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAg                                                                            >  2:2046506/1‑100 (MQ=255)
                        gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGc                                                                           >  1:2480064/1‑100 (MQ=255)
                         cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                          <  1:1046874/100‑1 (MQ=255)
                         cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                          <  2:1155993/100‑1 (MQ=255)
                         cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                          >  2:2517166/1‑100 (MQ=255)
                           aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTgg                                                                        >  1:1018369/1‑100 (MQ=255)
                           aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTg                                                                         >  1:709694/1‑99 (MQ=255)
                           aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGCATGGTTTTTTTTGTCTATAGCTg                                                                         >  1:2157929/1‑99 (MQ=255)
                             cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTa                                                                      <  2:1396882/100‑1 (MQ=255)
                              ggATTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTTTCAGAATTATGAAAAGTCATCCAGATTTGGATGGTTTTGTTTGTTTATAGCTGGTaa                                                                     <  1:1490776/100‑1 (MQ=255)
                               gAATTTTCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                     <  1:1885220/99‑1 (MQ=255)
                               gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                    <  1:1825361/100‑1 (MQ=255)
                               gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                    <  1:463873/100‑1 (MQ=255)
                                aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGa                                                                   <  2:1098729/100‑1 (MQ=255)
                                 aTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGAt                                                                  >  1:224775/1‑100 (MQ=255)
                                  tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTa                                                                      <  2:2341281/95‑1 (MQ=255)
                                  tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                 <  2:2400657/100‑1 (MQ=255)
                                  tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                 >  2:2279525/1‑100 (MQ=255)
                                  ttagCTGATTTGAAGCGGTGCGCCAGCGCAGTGAATTTCTCAGAATTATGAAAAGCCATCCAGATTTTGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                    <  1:1286645/94‑1 (MQ=255)
                                    tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                               >  1:1288100/1‑100 (MQ=255)
                                     gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                              <  2:1179905/100‑1 (MQ=255)
                                     gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                              >  2:704041/1‑100 (MQ=255)
                                      cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                             >  2:2042912/1‑100 (MQ=255)
                                        tATTTGAAGCGGTGCGCCAGCGAAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTATGATAAtt                                                              >  1:2050305/1‑97 (MQ=255)
                                          tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                         >  2:2386415/1‑100 (MQ=255)
                                           ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                        >  2:565219/1‑100 (MQ=255)
                                               aGCGTTGCGCCAACGCAACTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTTTAGCTGGTAAGGTAATTACGTAtt                                                       >  1:6616/1‑98 (MQ=255)
                                                 cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGTaaa                                                  <  2:2558940/100‑1 (MQ=255)
                                                 cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                  >  2:461620/1‑100 (MQ=255)
                                                    tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTTTCTAGCTGGTAAGATAATTACGCCTTGCAAatat                                               >  1:2462434/1‑100 (MQ=255)
                                                      cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTTCAAATATTc                                             >  1:1734251/1‑100 (MQ=255)
                                                        ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                           >  2:1923365/1‑100 (MQ=255)
                                                          aGCCCATCTAATTTCTCTTAATTATGAACAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                            <  1:1427796/97‑1 (MQ=255)
                                                             gcAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGAATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                      <  1:226628/100‑1 (MQ=255)
                                                             gcAGATAATTTCGCAGAAGTATGAAAAGCCATCGAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                      <  1:1248932/100‑1 (MQ=255)
                                                              cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                     >  2:1626318/1‑100 (MQ=255)
                                                              cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                     <  2:2315893/100‑1 (MQ=255)
                                                               aGCTGATTGCGCAGGATTATGAAAAGCCATCGAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                    <  2:1517929/100‑1 (MQ=255)
                                                               aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATTTTTTTTTTTTTATTTACGTATTGCAAATATTCATTttatta                                    >  1:2294716/1‑100 (MQ=255)
                                                                gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                   <  2:2230212/100‑1 (MQ=255)
                                                                  tAATTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAat                                  <  1:2561771/99‑1 (MQ=255)
                                                                  tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                 >  2:1054780/1‑100 (MQ=255)
                                                                  tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAat                                  <  1:1293795/99‑1 (MQ=255)
                                                                         tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                          >  2:663446/1‑100 (MQ=255)
                                                                         tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                          <  2:2118822/100‑1 (MQ=255)
                                                                         tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                          >  1:30211/1‑100 (MQ=255)
                                                                          cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                         <  2:35723/100‑1 (MQ=255)
                                                                          cAGAATTATGAAAAGCCATCCAGATATGGACGGCTTTTTTTAACTAGAGACAATAACAGAATTACGTATGGCAAATAGACAGGTTAGTAATAATTAAAc                          >  1:1120008/1‑99 (MQ=255)
                                                                           aGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                        >  1:368107/1‑100 (MQ=255)
                                                                            gAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                         <  1:1210026/98‑1 (MQ=255)
                                                                              aTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                     <  2:461567/100‑1 (MQ=255)
                                                                                caaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                   >  1:2433830/5‑100 (MQ=37)
                                                                                caaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATTTTATTAATAATTAAACTTAATCt                   >  1:887544/5‑100 (MQ=37)
                                                                                 aaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                  >  2:1172862/4‑100 (MQ=38)
                                                                                 aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTACACTTAATCta                  >  2:2474885/1‑100 (MQ=255)
                                                                                  tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTTTTTTTTAATAAGAAAATTAGTTATTGTATTTTTTCGTTTTATTGATTTTTAAACTTAATCttt                 >  1:636626/1‑98 (MQ=255)
                                                                                  tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTTTATTGCTGGTAAGATAATTACGTATTGCAAATATTCTTTTTATTAATATTTATACTTAATCt                   >  1:516704/1‑98 (MQ=255)
                                                                                  tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATACTTAGGTATTGCAAATAGTAATTCTCTTAATATTTACACTTAATCtat                 >  1:2122062/1‑100 (MQ=255)
                                                                                  aaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                 >  1:1477805/3‑100 (MQ=38)
                                                                                   gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                <  2:2297740/100‑1 (MQ=255)
                                                                                   aaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTTGTAAGATAATTACGAATTGCAAATATTCATTTTATTAATATTTATACTTAATCtata                >  1:1206940/2‑100 (MQ=37)
                                                                                    aaaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat               >  1:82609/1‑100 (MQ=255)
                                                                                    aaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATATCTGGTAAGATAATTACGTATTTCAAATATTCATTTTATTAATATTTAAACTTAATCTatat               >  1:910945/1‑100 (MQ=255)
                                                                                      aaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg             <  2:530181/100‑1 (MQ=255)
                                                                                          cctCCATATTTGGATGGTTTTTTTTGTCTATAGCTGGTAATATAATTACTTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa           <  1:1670935/96‑1 (MQ=255)
                                                                                             ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata      <  1:1319807/100‑1 (MQ=255)
                                                                                                gATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat   >  2:1010532/1‑100 (MQ=255)
                                                                                                 aTTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg  <  1:104601/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATG  >  NC_000913/4409938‑4410135

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