Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,874,594 A→G intergenic (‑14/+101) yidA ← / ← b3698 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,874,5940AG83.6% 629.9 / 112.4 256intergenic (‑14/+101)yidA/b3698hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/17);  new base G (121/93);  total (146/110)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

Only 100 of 261 total aligned reads displayed.

AACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCC  >  NC_000913/3874506‑3874686
                                                                                        |                                                                                            
aaCGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGAtt                                                                                          <  1:5342169/93‑1 (MQ=255)
aaCGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGAtt                                                                                          <  1:3435171/93‑1 (MQ=255)
aaCGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGAtt                                                                                          <  1:3688449/93‑1 (MQ=255)
     cggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTc                                                                                     <  1:6604847/93‑1 (MQ=255)
     cggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTc                                                                                     <  2:4881536/93‑1 (MQ=255)
     cggGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTc                                                                                     >  2:355290/1‑93 (MQ=255)
      gtgTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTcc                                                                                    <  2:1541788/91‑1 (MQ=255)
      gggTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTcc                                                                                    <  2:5541500/93‑1 (MQ=255)
       ggTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCt                                                                                   <  2:5010182/93‑1 (MQ=255)
        gTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTg                                                                                  >  2:4457418/1‑93 (MQ=255)
           aaaTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGAcc                                                                               >  1:2030263/1‑93 (MQ=255)
            aaTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGAccc                                                                              >  2:1727974/1‑93 (MQ=255)
             aTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCa                                                                             <  1:1166646/93‑1 (MQ=255)
             aTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCa                                                                             <  2:3259671/93‑1 (MQ=255)
             aTGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCa                                                                             <  2:6088686/93‑1 (MQ=255)
              tGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCAt                                                                            <  1:4518377/93‑1 (MQ=255)
              tGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCAt                                                                            <  2:2497190/93‑1 (MQ=255)
               ggTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATg                                                                           >  1:6486187/1‑93 (MQ=255)
                 tgtgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCt                                                                         <  1:4563706/93‑1 (MQ=255)
                  gtgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTa                                                                        >  2:3487947/1‑93 (MQ=255)
                   tgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTaa                                                                       >  2:4651532/1‑93 (MQ=255)
                   tgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTaa                                                                       >  1:772026/1‑93 (MQ=255)
                   tgATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTaa                                                                       >  1:3286143/1‑93 (MQ=255)
                    gATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAAc                                                                      >  1:1486705/1‑93 (MQ=255)
                     aTCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAAcg                                                                     >  1:3617180/1‑93 (MQ=255)
                       cGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAAcgcg                                                                   >  2:2390766/1‑93 (MQ=255)
                       cGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAAcgcg                                                                   >  1:2763404/1‑93 (MQ=255)
                       cGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAAcgcg                                                                   <  1:445222/93‑1 (MQ=255)
                         ggCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGAt                                                                 >  2:620730/1‑93 (MQ=255)
                           cagcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTc                                                               <  2:5215895/93‑1 (MQ=255)
                             gcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCg                                                             >  2:6423644/1‑93 (MQ=255)
                             gcagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCg                                                             >  1:1747437/1‑93 (MQ=255)
                              cagAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGc                                                            <  1:3970700/93‑1 (MQ=255)
                                gAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTc                                                          >  1:4190486/1‑93 (MQ=255)
                                gAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTc                                                          <  2:186272/93‑1 (MQ=255)
                                gAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTc                                                          <  2:3167794/93‑1 (MQ=255)
                                gAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGCTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTc                                                          >  2:1050478/1‑93 (MQ=255)
                                   gggTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCaaa                                                       >  2:1255950/1‑93 (MQ=255)
                                   gggTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCaaa                                                       >  2:1133468/1‑93 (MQ=255)
                                     gTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCaaaaa                                                     <  1:987008/93‑1 (MQ=255)
                                       gCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATc                                                   >  2:5308855/1‑93 (MQ=255)
                                        ccatccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCa                                                  <  1:3866414/93‑1 (MQ=255)
                                        ccatccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCa                                                  <  2:3590948/93‑1 (MQ=255)
                                         catccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAACATCAg                                                 >  2:2732489/1‑93 (MQ=255)
                                         catccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAg                                                 <  1:1757929/93‑1 (MQ=255)
                                          atccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGt                                                >  2:2569490/1‑93 (MQ=255)
                                           tccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTa                                               <  2:4286064/93‑1 (MQ=255)
                                            ccatATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTac                                              >  2:3731253/1‑93 (MQ=255)
                                             catATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTaca                                             >  1:1945973/1‑93 (MQ=255)
                                                atCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTacaaca                                          >  2:1947144/1‑93 (MQ=255)
                                                 tCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAacac                                         <  2:5164023/93‑1 (MQ=255)
                                                  cGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACAcc                                        <  1:5574143/93‑1 (MQ=255)
                                                   gATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACAccc                                       >  2:2217585/1‑93 (MQ=255)
                                                   gATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACAccc                                       >  1:5055407/1‑93 (MQ=255)
                                                    aTAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCg                                      >  1:215455/1‑93 (MQ=255)
                                                     tAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGa                                     >  1:6377825/1‑93 (MQ=255)
                                                     tAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGa                                     >  1:2231377/1‑93 (MQ=255)
                                                      aGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAg                                    <  2:4397585/93‑1 (MQ=255)
                                                      aGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAg                                    >  2:4962003/1‑93 (MQ=255)
                                                      aGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAg                                    <  2:3976412/93‑1 (MQ=255)
                                                       gCAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAgg                                   >  2:6330527/1‑93 (MQ=255)
                                                        cAATGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAggg                                  >  1:4740920/1‑93 (MQ=255)
                                                         aaTGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGa                                 >  1:5447674/1‑93 (MQ=255)
                                                          aTGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAACAATCAGTACAACACCCGAGGGaa                                >  2:6012238/1‑93 (MQ=255)
                                                           tGAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaa                               >  2:2895293/1‑93 (MQ=255)
                                                            gAGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaaa                              >  1:2673202/1‑93 (MQ=255)
                                                             aGTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAg                             <  2:40360/93‑1 (MQ=255)
                                                              gTTTAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgg                            >  1:809045/1‑93 (MQ=255)
                                                               tttAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACTCCCGAGGGAAAAggg                           <  1:5830226/93‑1 (MQ=255)
                                                               tttAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggg                           <  1:4910846/93‑1 (MQ=255)
                                                                ttAATAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggag                          >  1:6259588/1‑91 (MQ=255)
                                                                 tAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCATTACACCCCCCGAGGGAAAAggggg                         <  1:2177518/93‑1 (MQ=255)
                                                                  aaTAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGa                        >  2:4233163/1‑93 (MQ=255)
                                                                  aaTAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGa                        <  1:4255235/93‑1 (MQ=255)
                                                                  aaTAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGa                        >  2:5394583/1‑93 (MQ=255)
                                                                  aaTAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGa                        >  1:1028403/1‑93 (MQ=255)
                                                                   aTAGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGAt                       >  1:4800894/1‑93 (MQ=255)
                                                                   aTAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGAt                       <  2:34102/93‑1 (MQ=255)
                                                                     aGCCATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGc                     >  1:1588221/1‑93 (MQ=255)
                                                                      gCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCa                    >  1:1299229/1‑93 (MQ=255)
                                                                        cATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAAc                  >  1:5760968/1‑93 (MQ=255)
                                                                        cATGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAAc                  >  1:4969658/1‑93 (MQ=255)
                                                                        cATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAAc                  >  2:1639340/1‑93 (MQ=255)
                                                                        cATGAGATCCCCAGGGGGATTGGTTCCTGACCCATGCTAACGCGATTCAGCTCAAAAATCAGTACACCACCCGAGGGAAAAGGGGGATGCAAc                  >  1:3162726/1‑93 (MQ=255)
                                                                          tGAGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAAcgc                >  2:217047/1‑93 (MQ=255)
                                                                           gAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAAcgcg               >  2:1888950/1‑93 (MQ=255)
                                                                            aGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAAcgcgc              >  1:6472118/1‑93 (MQ=255)
                                                                            aGTTCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAAcgcgc              <  1:3542564/93‑1 (MQ=255)
                                                                             gTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAAcgcgcg             >  1:1550267/1‑93 (MQ=255)
                                                                              ttCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGt            <  1:6121385/93‑1 (MQ=255)
                                                                               tCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGtt           <  2:5871382/93‑1 (MQ=255)
                                                                               tCCCCAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGtt           <  2:1489472/93‑1 (MQ=255)
                                                                                ccccAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTg          <  2:3115121/93‑1 (MQ=255)
                                                                                 cccAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGc         >  2:1981686/1‑93 (MQ=255)
                                                                                  ccAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCg        >  1:2512609/1‑93 (MQ=255)
                                                                                   cAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGt       >  2:6102612/1‑93 (MQ=255)
                                                                                   cAGGTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGt       <  2:670809/93‑1 (MQ=255)
                                                                                   cAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGt       <  1:1329645/93‑1 (MQ=255)
                                                                                     ggTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGc     >  2:3497766/1‑93 (MQ=255)
                                                                                     ggTGGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGc     >  2:2971773/1‑93 (MQ=255)
                                                                                        |                                                                                            
AACGGCGGGTGAAATGGTGTGATCGGGCAGCAGAAGGGTGCCATCCATATCGATAGCAATGAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCC  >  NC_000913/3874506‑3874686

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