Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 206,483 A→G Y453C (TAC→TGC)  dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913206,4830AG98.4% 383.4 / ‑3.9 124Y453C (TAC→TGC) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (63/59);  total (63/61)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 131 total aligned reads displayed.

GGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTACGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGCTGCCG  >  NC_000913/206392‑206571
                                                                                           |                                                                                        
ggTCGTGATGCGGTATCGCAGATCATCACCTTCGGGACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGGCGCGTGCTGGGGCATCCGTAc                                                                                         >  1:277597/1‑93 (MQ=255)
  tCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCgg                                                                                       >  1:1732224/1‑93 (MQ=255)
   cGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGc                                                                                      >  1:4185029/1‑93 (MQ=255)
    gTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCt                                                                                     <  1:847172/93‑1 (MQ=255)
     tGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCtt                                                                                    <  1:2836312/93‑1 (MQ=255)
     ggATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCtt                                                                                    <  2:3671387/92‑1 (MQ=255)
      gATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCttt                                                                                   >  2:201613/1‑93 (MQ=255)
      gATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCttt                                                                                   <  1:3860358/93‑1 (MQ=255)
       aTGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTg                                                                                  >  2:4312601/1‑93 (MQ=255)
       aTGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTg                                                                                  >  2:3142830/1‑93 (MQ=255)
        tGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGACGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGt                                                                                 <  2:3099118/93‑1 (MQ=255)
         gCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTc                                                                                >  1:559685/1‑93 (MQ=255)
         gCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTc                                                                                <  2:3162907/93‑1 (MQ=255)
         gCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTc                                                                                >  2:103838/1‑93 (MQ=255)
           ggTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGa                                                                              <  1:3882373/93‑1 (MQ=255)
           ggTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGa                                                                              <  2:227118/93‑1 (MQ=255)
            gTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGAt                                                                             >  2:3504632/1‑93 (MQ=255)
            gTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGAt                                                                             >  1:29340/1‑93 (MQ=255)
             tATCGCCGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATc                                                                            <  1:1541113/93‑1 (MQ=255)
             tATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATc                                                                            >  1:378866/1‑93 (MQ=255)
              aTCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCg                                                                           >  1:1793519/1‑93 (MQ=255)
               tCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGt                                                                          >  1:1613805/1‑93 (MQ=255)
                 gCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTAt                                                                        >  2:1385626/1‑93 (MQ=255)
                  cAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTAtc                                                                       <  1:4055761/93‑1 (MQ=255)
                  cAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTACGGCTTTGTCGATCGTAtc                                                                       <  1:4366573/93‑1 (MQ=255)
                   aGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTAtct                                                                      >  2:2589172/1‑93 (MQ=255)
                   aGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTAtct                                                                      >  1:146038/1‑93 (MQ=255)
                    gATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTAtctc                                                                     <  1:102284/93‑1 (MQ=255)
                      tcatcaCCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGa                                                                   >  1:1429518/1‑93 (MQ=255)
                          caCCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACt                                                               >  2:4172107/1‑93 (MQ=255)
                          caCCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACt                                                               >  2:2056134/1‑93 (MQ=255)
                           cccTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTg                                                              <  1:2543650/92‑1 (MQ=255)
                           aCCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCGTCCGTGCGGCTTTGTCGATCGGATCTCGAAACTg                                                              >  2:1328127/1‑93 (MQ=255)
                           aCCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTg                                                              >  2:87355/1‑93 (MQ=255)
                             cTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGAt                                                            >  2:4314208/1‑93 (MQ=255)
                              ttCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATc                                                           <  2:284186/93‑1 (MQ=255)
                               tCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGCAACTGATcc                                                          >  2:3673087/1‑93 (MQ=255)
                               tCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATcc                                                          <  1:398267/93‑1 (MQ=255)
                               tCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATcc                                                          <  1:2922/93‑1 (MQ=255)
                               tCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATcc                                                          <  1:1277556/93‑1 (MQ=255)
                                cGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATccc                                                         <  1:234063/93‑1 (MQ=255)
                                cGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATccc                                                         <  1:891162/93‑1 (MQ=255)
                                  gTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATcccgc                                                       >  2:1210778/1‑93 (MQ=255)
                                  gTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATcccgc                                                       >  2:4745958/1‑93 (MQ=255)
                                   tACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATcccgcc                                                      <  2:2020950/93‑1 (MQ=255)
                                   tACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATcccgcc                                                      <  1:904939/93‑1 (MQ=255)
                                    aCAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATcccgccc                                                     <  1:1028104/93‑1 (MQ=255)
                                    aCAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATcccgccc                                                     >  2:2814998/1‑93 (MQ=255)
                                      aaTGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGa                                                   >  1:4569544/1‑93 (MQ=255)
                                       aTGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGAt                                                  >  2:3920684/1‑93 (MQ=255)
                                       aTGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGAt                                                  >  1:2216580/1‑93 (MQ=255)
                                         ggcggcGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATcc                                                <  2:1284522/93‑1 (MQ=255)
                                         ggcggcGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATATCGAAACTGATCCCGCCCGATcc                                                >  2:590775/1‑93 (MQ=255)
                                         ggcggcGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGATTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATcc                                                >  1:1041225/1‑93 (MQ=255)
                                          gcggcgAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCg                                               <  1:4106021/93‑1 (MQ=255)
                                            ggcgAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCggg                                             <  2:907978/93‑1 (MQ=255)
                                              cgAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGa                                           <  2:940147/93‑1 (MQ=255)
                                              cgAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGa                                           >  2:2561386/1‑93 (MQ=255)
                                               gAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGAt                                          >  1:3045552/1‑93 (MQ=255)
                                               gAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGAt                                          <  1:1126392/93‑1 (MQ=255)
                                                 aaGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGa                                        <  2:885763/93‑1 (MQ=255)
                                                  aGCAGTGATCCGCGCCGAAGACCGCGTGCTGGGGCATCCGTGCGGCTTTGTCAATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGAc                                       <  2:4077132/93‑1 (MQ=255)
                                                    cGGTGAGCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGc                                     <  1:3450393/93‑1 (MQ=255)
                                                       tggtCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTACGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTgg                                  <  2:3829351/90‑1 (MQ=255)
                                                       tGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTgg                                  >  2:1230945/1‑93 (MQ=255)
                                                        gATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGc                                 <  1:2905791/93‑1 (MQ=255)
                                                         aTCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCg                                <  2:4290186/93‑1 (MQ=255)
                                                         aTCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCg                                >  1:1451394/1‑93 (MQ=255)
                                                            cgcgACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGaaa                             <  2:4648706/93‑1 (MQ=255)
                                                             gcgACGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAg                            >  2:3200585/1‑93 (MQ=255)
                                                                aCGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGt                         >  2:256188/1‑93 (MQ=255)
                                                                 cGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGtt                        <  2:886933/93‑1 (MQ=255)
                                                                 cGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGtt                        <  2:3571035/93‑1 (MQ=255)
                                                                 cGTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGtt                        <  1:1119537/93‑1 (MQ=255)
                                                                  gTAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGttt                       <  2:2245188/93‑1 (MQ=255)
                                                                   tAGGCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTg                      <  2:3390637/93‑1 (MQ=255)
                                                                     ggCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGaa                    >  1:2738224/1‑93 (MQ=255)
                                                                     ggCCGCGTGCTGGGGCAGCCGTGCGTCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGaa                    <  1:1562194/93‑1 (MQ=255)
                                                                      gCCGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAag                   <  1:4007200/93‑1 (MQ=255)
                                                                       ccGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAagc                  <  1:4280934/93‑1 (MQ=255)
                                                                       ccGCGTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAagc                  <  1:2576841/93‑1 (MQ=255)
                                                                        cgcgTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAagcc                 <  2:931996/93‑1 (MQ=255)
                                                                         gcgTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGCAAGCGTTTGAagccg                >  2:2237220/1‑93 (MQ=255)
                                                                         gcgTGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAagccg                <  1:1057067/93‑1 (MQ=255)
                                                                            tGCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAagccgagc             <  1:4393764/93‑1 (MQ=255)
                                                                             gCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAagccgagcc            >  1:604444/1‑93 (MQ=255)
                                                                             gCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAagccgagcc            >  1:156372/1‑93 (MQ=255)
                                                                             gCTGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAagccgagcc            <  1:1735304/93‑1 (MQ=255)
                                                                               tGGGGCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGc          >  1:2536194/1‑93 (MQ=255)
                                                                                ggggCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCa         >  1:2428670/1‑93 (MQ=255)
                                                                                  ggCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGc       <  1:4403019/93‑1 (MQ=255)
                                                                                  ggCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGGTTGAAGCCGAGCCGCAGc       >  1:2476669/1‑93 (MQ=255)
                                                                                   gCATCCGTGCGGTTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGCt      <  1:3162195/93‑1 (MQ=255)
                                                                                   gCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGCt      >  2:763639/1‑93 (MQ=255)
                                                                                   gCATCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGCt      <  2:2915994/93‑1 (MQ=255)
                                                                                     aTCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGCTGc    >  2:803607/1‑93 (MQ=255)
                                                                                      tCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTTAAGCCGAGCCTCAGGTGc  >  2:3504551/1‑92 (MQ=255)
                                                                                      tCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGCTGcc   <  2:4572958/93‑1 (MQ=255)
                                                                                      tCCGTGCGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGCTGcc   >  2:1605356/1‑93 (MQ=255)
                                                                                       ccGTGCGGCTTTGTCGATCGTATCTCTAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGCTGCCg  <  1:4019618/93‑1 (MQ=255)
                                                                                           |                                                                                        
GGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGCATCCGTACGGCTTTGTCGATCGTATCTCGAAACTGATCCCGCCCGATCCGGGGATGACGCTGGCGAAAGCGTTTGAAGCCGAGCCGCAGCTGCCG  >  NC_000913/206392‑206571

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