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)
        ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTT‑AAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtt                                                                                           >  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)
                    ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCta                                                                               >  1:110390/1‑100 (MQ=255)
                    ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCta                                                                               >  2:1489052/1‑100 (MQ=255)
                    ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataattc                                                                                <  2:1946639/100‑9 (MQ=25)
                    ggTTGCGACCGGAATTTGCTTAGTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataattc                                                                                <  1:1774811/100‑9 (MQ=25)
                     gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCtat                                                                              <  1:879542/100‑1 (MQ=255)
                     gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCtat                                                                              <  1:2088872/100‑1 (MQ=255)
                      ttGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCtata                                                                             <  2:1492072/100‑1 (MQ=255)
                       tGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtcttt                                                                               >  1:1973761/1‑92 (MQ=255)
                       tGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAg                                                                            >  1:1628861/1‑100 (MQ=255)
                        gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAg                                                                            >  2:2046506/1‑100 (MQ=255)
                        gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGc                                                                           >  1:2480064/1‑100 (MQ=255)
                         cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCt                                                                          <  1:1046874/100‑1 (MQ=255)
                         cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCt                                                                          <  2:1155993/100‑1 (MQ=255)
                         cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCt                                                                          >  2:2517166/1‑100 (MQ=255)
                           aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTgg                                                                        >  1:1018369/1‑100 (MQ=255)
                           aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTg                                                                         >  1:709694/1‑99 (MQ=255)
                           aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGCATGG‑TTTTTTTTGTCTATAGCTg                                                                         >  1:2157929/1‑99 (MQ=255)
                             cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTa                                                                      <  2:1396882/100‑1 (MQ=255)
                              ggATTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTTTCAGAATTATGAAAAGTCATCCAGATTTGGATGG‑TTTTGTTTGTTTATAGCTGGTaa                                                                     <  1:1490776/100‑1 (MQ=255)
                               gAATTTTCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTaa                                                                     <  1:1885220/99‑1 (MQ=255)
                               gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAg                                                                    <  1:1825361/100‑1 (MQ=255)
                               gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAg                                                                    <  1:463873/100‑1 (MQ=255)
                                aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGa                                                                   <  2:1098729/100‑1 (MQ=255)
                                 aTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGAt                                                                  >  1:224775/1‑100 (MQ=255)
                                  tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTa                                                                      <  2:2341281/95‑1 (MQ=255)
                                  tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                 <  2:2400657/100‑1 (MQ=255)
                                  tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                 >  2:2279525/1‑100 (MQ=255)
                                  ttagCTGATTTGAAGCGGTGCGCCAGCGCAGTGAATTTCTCAGAATTATGAAAAGCCATCCAGATTTTGATGG‑TTTTTTTTGTCTATAGCTGGTAAg                                                                    <  1:1286645/94‑1 (MQ=255)
                                    tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAAt                                                               >  1:1288100/1‑100 (MQ=255)
                                     gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                              <  2:1179905/100‑1 (MQ=255)
                                     gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                              >  2:704041/1‑100 (MQ=255)
                                      cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                             >  2:2042912/1‑100 (MQ=255)
                                        tATTTGAAGCGGTGCGCCAGCGAAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTATGATAAtt                                                              >  1:2050305/1‑97 (MQ=255)
                                          tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                         >  2:2386415/1‑100 (MQ=255)
                                           ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                        >  2:565219/1‑100 (MQ=255)
                                               aGCGTTGCGCCAACGCAACTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTTTAGCTGGTAAGGTAATTACGTAtt                                                       >  1:6616/1‑98 (MQ=255)
                                                 cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGTaaa                                                  <  2:2558940/100‑1 (MQ=255)
                                                 cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                  >  2:461620/1‑100 (MQ=255)
                                                    tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTTTCTAGCTGGTAAGATAATTACGCCTTGCAAatat                                               >  1:2462434/1‑100 (MQ=255)
                                                      cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTTCAAATATTc                                             >  1:1734251/1‑100 (MQ=255)
                                                        ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                           >  2:1923365/1‑100 (MQ=255)
                                                          aGCCCATCTAATTTCTCTTAATTATGAACAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                            <  1:1427796/97‑1 (MQ=255)
                                                             gcAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGAATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                      <  1:226628/100‑1 (MQ=255)
                                                             gcAGATAATTTCGCAGAAGTATGAAAAGCCATCGAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                      <  1:1248932/100‑1 (MQ=255)
                                                              cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                     >  2:1626318/1‑100 (MQ=255)
                                                              cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                     <  2:2315893/100‑1 (MQ=255)
                                                               aGCTGATTGCGCAGGATTATGAAAAGCCATCGAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                    <  2:1517929/100‑1 (MQ=255)
                                                               aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATTTTTTTTTTTTTATTTACGTATTGCAAATATTCATTttatta                                    >  1:2294716/1‑100 (MQ=255)
                                                                gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                   <  2:2230212/100‑1 (MQ=255)
                                                                  tAATTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAat                                  <  1:2561771/99‑1 (MQ=255)
                                                                  tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                 >  2:1054780/1‑100 (MQ=255)
                                                                  tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAat                                  <  1:1293795/99‑1 (MQ=255)
                                                                         tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                          >  2:663446/1‑100 (MQ=255)
                                                                         tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                          <  2:2118822/100‑1 (MQ=255)
                                                                         tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                          >  1:30211/1‑100 (MQ=255)
                                                                          cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                         <  2:35723/100‑1 (MQ=255)
                                                                          cAGAATTATGAAAAGCCATCCAGATATGGACGG‑CTTTTTTTAACTAGAGACAATAACAGAATTACGTATGGCAAATAGACAGGTTAGTAATAATTAAAc                          >  1:1120008/1‑99 (MQ=255)
                                                                           aGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                        >  1:368107/1‑100 (MQ=255)
                                                                            gAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                         <  1:1210026/98‑1 (MQ=255)
                                                                              aTTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                     <  2:461567/100‑1 (MQ=255)
                                                                                caaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                   >  1:2433830/5‑100 (MQ=37)
                                                                                caaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATTTTATTAATAATTAAACTTAATCt                   >  1:887544/5‑100 (MQ=37)
                                                                                 aaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                  >  2:1172862/4‑100 (MQ=38)
                                                                                 aTGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTACACTTAATCta                  >  2:2474885/1‑100 (MQ=255)
                                                                                  tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTT‑TTTTTTTTTAATAAGAAAATTAGTTATTGTATTTTTTCGTTTTATTGATTTTTAAACTTAATCttt                 >  1:636626/1‑98 (MQ=255)
                                                                                  tGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTTTATTGCTGGTAAGATAATTACGTATTGCAAATATTCTTTTTATTAATATTTATACTTAATCt                   >  1:516704/1‑98 (MQ=255)
                                                                                  tGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATACTTAGGTATTGCAAATAGTAATTCTCTTAATATTTACACTTAATCtat                 >  1:2122062/1‑100 (MQ=255)
                                                                                  aaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                 >  1:1477805/3‑100 (MQ=38)
                                                                                   gAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                <  2:2297740/100‑1 (MQ=255)
                                                                                   aaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTTGTAAGATAATTACGAATTGCAAATATTCATTTTATTAATATTTATACTTAATCtata                >  1:1206940/2‑100 (MQ=37)
                                                                                    aaaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat               >  1:82609/1‑100 (MQ=255)
                                                                                    aaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATATCTGGTAAGATAATTACGTATTTCAAATATTCATTTTATTAATATTTAAACTTAATCTatat               >  1:910945/1‑100 (MQ=255)
                                                                                      aaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg             <  2:530181/100‑1 (MQ=255)
                                                                                          cctCCATATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAATATAATTACTTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa           <  1:1670935/96‑1 (MQ=255)
                                                                                             ccAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata      <  1:1319807/100‑1 (MQ=255)
                                                                                                gATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat   >  2:1010532/1‑100 (MQ=255)
                                                                                                 aTTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg  <  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: ‑