Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,179,001 C→T S41S (AGC→AGT rplA → 50S ribosomal subunit protein L1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,179,0010CT96.4% 1662.6 / 32.3 504S41S (AGC→AGTrplA50S ribosomal subunit protein L1
Reads supporting (aligned to +/- strand):  ref base C (11/7);  new base T (273/212);  total (285/219)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 507 total aligned reads displayed.

GGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGC  >  NC_000913/4178854‑4179151
                                                                                                                                                   |                                                                                                                                                      
gggCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCGAAATTCGTAGAAAGtgtg                                                                                                                                                     >  4:1474068/1‑151 (MQ=255)
  gCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGa                                                                                                                                                   <  5:1527831/151‑1 (MQ=255)
    cTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACg                                                                                                                                                 >  7:1370826/1‑151 (MQ=255)
      ggTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGt                                                                                                                                                <  1:1435847/150‑1 (MQ=255)
       gTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAg                                                                                                                                              <  3:881094/151‑1 (MQ=255)
         aGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCt                                                                                                                                            >  3:1006025/1‑151 (MQ=255)
             gaggaCTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTa                                                                                                                                        <  1:1259050/151‑1 (MQ=255)
              aggaCTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTa                                                                                                                                        >  6:1295466/1‑150 (MQ=255)
              aggaCTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTa                                                                                                                                        >  4:128525/1‑150 (MQ=255)
                   tAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCg                                                                                                                                  >  3:689103/1‑151 (MQ=255)
                    aaGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTc                                                                                                                                   <  1:378901/149‑1 (MQ=255)
                    aaGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCgg                                                                                                                                 <  4:1215253/151‑1 (MQ=255)
                    aaGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCg                                                                                                                                  <  7:1281884/150‑1 (MQ=255)
                    aaGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTc                                                                                                                                   <  3:1376000/149‑1 (MQ=255)
                     aGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGc                                                                                                                                <  6:1535785/151‑1 (MQ=255)
                     aGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATACGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGc                                                                                                                                >  5:310431/1‑151 (MQ=255)
                       aaaTGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCa                                                                                                                               <  4:1644330/150‑1 (MQ=255)
                         aTGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCg                                                                                                                            >  8:569761/1‑151 (MQ=255)
                           ggtttaaCTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGa                                                                                                                           <  4:755402/145‑1 (MQ=255)
                           ggttAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGAc                                                                                                                          >  7:858547/4‑151 (MQ=255)
                           ggttAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGAc                                                                                                                          >  1:288723/4‑151 (MQ=255)
                             ttAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGc                                                                                                                        >  8:1326157/2‑151 (MQ=255)
                             cTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGc                                                                                                                        >  8:878949/1‑151 (MQ=255)
                               aaaCTGACCAAGCGCATGGGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTc                                                                                                                      <  2:1223646/151‑1 (MQ=255)
                                 aCTGCCCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCTCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACTCTCGt                                                                                                                    >  8:355614/1‑151 (MQ=255)
                                 aCTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGt                                                                                                                    <  8:1082755/151‑1 (MQ=255)
                                 aCTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGt                                                                                                                    >  6:480855/1‑151 (MQ=255)
                                  cTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGt                                                                                                                    >  8:1273959/1‑150 (MQ=255)
                                  cTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGt                                                                                                                    >  5:1691034/1‑150 (MQ=255)
                                   tGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTa                                                                                                                   >  7:566444/1‑150 (MQ=255)
                                    gACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTa                                                                                                                   <  7:1052165/149‑1 (MQ=255)
                                      ccAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTACCCTCGGCATCGACGCTCGTAAATc                                                                                                               >  3:770340/1‑151 (MQ=255)
                                       cAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCt                                                                                                              >  3:1093065/1‑151 (MQ=255)
                                        aaGCGCATGTGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTg                                                                                                             >  1:1554661/1‑151 (MQ=255)
                                        aaGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCtt                                                                                                             >  2:87317/1‑150 (MQ=255)
                                         aGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTg                                                                                                             >  1:859762/1‑150 (MQ=255)
                                          gcgcATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGAc                                                                                                           <  3:1627603/151‑1 (MQ=255)
                                           cgcATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGAcc                                                                                                          <  1:1166065/151‑1 (MQ=255)
                                             cATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAg                                                                                                        <  2:1390985/151‑1 (MQ=255)
                                                  gtgtTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAacgt                                                                                                   <  4:1076819/151‑1 (MQ=255)
                                                   tgtTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAacgt                                                                                                   >  8:628656/1‑150 (MQ=255)
                                                     ttATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAacgtacg                                                                                                >  7:1628006/1‑151 (MQ=255)
                                                      tATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAacgtacgt                                                                                               >  6:303278/1‑151 (MQ=255)
                                                      tATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAacgtacgt                                                                                               >  7:1215014/1‑151 (MQ=255)
                                                       aTCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACgtg                                                                                              >  4:647762/1‑151 (MQ=255)
                                                         ccGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACgtgg                                                                                             <  4:1000722/150‑1 (MQ=255)
                                                             gagaAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCa                                                                                         >  1:671362/1‑150 (MQ=255)
                                                              agaAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAAc                                                                                       >  1:1075041/1‑151 (MQ=255)
                                                                    ttGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACt                                                                                 >  1:1005275/1‑151 (MQ=255)
                                                                      gATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTg                                                                                <  8:694133/150‑1 (MQ=255)
                                                                       aTGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGAATCGACGCTAGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGcc                                                                              >  1:1389612/1‑151 (MQ=255)
                                                                         gCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCCGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCg                                                                             <  2:1643704/150‑1 (MQ=255)
                                                                               aaaCAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTa                                                                      <  3:573982/151‑1 (MQ=255)
                                                                                aaCAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTAc                                                                     >  2:1161091/1‑151 (MQ=255)
                                                                                 agaGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTATATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACt                                                                    <  5:656932/149‑1 (MQ=255)
                                                                                 aCAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACt                                                                    <  1:1629724/151‑1 (MQ=255)
                                                                                 aCAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACt                                                                    <  6:680326/151‑1 (MQ=255)
                                                                                   aGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTgg                                                                  >  2:564400/1‑151 (MQ=255)
                                                                                   aGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTgg                                                                  <  6:1083703/151‑1 (MQ=255)
                                                                                    gTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGc                                                                 >  1:908994/1‑151 (MQ=255)
                                                                                       cGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGccgt                                                              <  3:12440/151‑1 (MQ=255)
                                                                                           aTCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGccgttccg                                                          >  6:923631/1‑151 (MQ=255)
                                                                                                cGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTcgc                                                     >  4:1557650/1‑151 (MQ=255)
                                                                                                  aaGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGt                                                   <  5:136251/151‑1 (MQ=255)
                                                                                                  aaGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGt                                                   <  5:498808/151‑1 (MQ=255)
                                                                                                  aaGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGt                                                   <  4:392234/151‑1 (MQ=255)
                                                                                                  aaGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGt                                                   >  7:1575809/1‑151 (MQ=255)
                                                                                                    gCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCTCGGTACTGGCCGATCCGTTCGCGt                                                   >  4:236790/1‑149 (MQ=255)
                                                                                                     cTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGc                                                >  8:234727/1‑151 (MQ=255)
                                                                                                       aTCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGcc                                               <  6:1191059/150‑1 (MQ=255)
                                                                                                         cGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTACCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGt                                             <  7:346743/150‑1 (MQ=255)
                                                                                                                ctgAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTAccc                                     >  8:1683469/1‑151 (MQ=255)
                                                                                                                 tgAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGCCCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTAccc                                     >  3:1497403/1‑150 (MQ=255)
                                                                                                                  gAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTAccc                                     <  4:1376963/149‑1 (MQ=255)
                                                                                                                  gAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCa                                    >  7:471914/1‑150 (MQ=255)
                                                                                                                   aaaGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCa                                    >  5:1001078/1‑149 (MQ=255)
                                                                                                                      gagCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTg                               <  5:1086702/151‑1 (MQ=255)
                                                                                                                      gagCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTg                               <  6:857858/151‑1 (MQ=255)
                                                                                                                       agCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTACCCTCGGCATCGACGCTCGTAAATCTGACCACAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGc                              >  7:998040/1‑151 (MQ=255)
                                                                                                                       agCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGc                              >  6:582844/1‑151 (MQ=255)
                                                                                                                         cTGGCGACTCCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATAGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGACGCACGGTACTGGACGTTACGTTCGCGTAGCCGTATTTACCAAAGGTGCa                             >  5:1133045/1‑150 (MQ=255)
                                                                                                                          tGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCaa                            >  2:425371/1‑150 (MQ=255)
                                                                                                                              gACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTAGCGTAGCCGTATTTACCCAAGGTGCAAACGCt                       >  6:172845/1‑151 (MQ=255)
                                                                                                                                 tGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGa                     >  5:981961/1‑150 (MQ=255)
                                                                                                                                  gCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTNACCTCGGCATCGACGNTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTNTTTACCCAAGGTGCAAACGCTGAAg                   <  8:402824/151‑1 (MQ=255)
                                                                                                                                  gCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGa                     <  7:244027/149‑1 (MQ=255)
                                                                                                                                  gCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAg                   <  5:826432/151‑1 (MQ=255)
                                                                                                                                   cTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAgc                  >  8:806508/1‑151 (MQ=255)
                                                                                                                                    tAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAgct                 >  7:526226/1‑151 (MQ=255)
                                                                                                                                    tAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTCCAACTGTACTGCCGCACGGAACTGGCCGTTCCGATCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAgct                 >  7:409780/1‑151 (MQ=255)
                                                                                                                                      aaTTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAgctgc               >  7:1465378/1‑151 (MQ=255)
                                                                                                                                       aTTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAgctgct              >  5:164585/1‑151 (MQ=255)
                                                                                                                                       aTTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAgctgct              >  7:786725/1‑151 (MQ=255)
                                                                                                                                       aTTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAgctgct              >  5:1276524/1‑151 (MQ=255)
                                                                                                                                       aTTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACACTGAAgctgct              >  8:1615720/1‑151 (MQ=255)
                                                                                                                                        ttCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTCGCCGTATTTACCCAAGGTGCAAACGCTGAAgctgct              >  4:542976/1‑150 (MQ=255)
                                                                                                                                         tCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGAACGTGGTGCAACTGTACTGCCGCACGGTACTCGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTCCTa             >  2:1010012/1‑150 (MQ=255)
                                                                                                                                            tAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGc         >  1:477229/1‑151 (MQ=255)
                                                                                                                                            tAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTCAAGCTGATAAAGc         >  3:463752/1‑151 (MQ=255)
                                                                                                                                             aGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGc         <  3:263905/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGC  >  NC_000913/4178854‑4179151

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: