Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,253,781 T→C E171G (GAG→GGG)  yhaNM ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,253,7810TC99.1% 366.2 / ‑5.0 116E171G (GAG→GGG) yhaNMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (60/55);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 122 total aligned reads displayed.

TCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACAC  >  NC_000913/3253691‑3253870
                                                                                          |                                                                                         
tCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTg                                                                                        <  1:3992164/93‑1 (MQ=255)
tCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTg                                                                                        <  2:4250559/93‑1 (MQ=255)
tCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTg                                                                                        <  1:2766702/93‑1 (MQ=255)
  gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                                                                      >  1:140221/1‑93 (MQ=255)
  gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                                                                      <  2:2054139/93‑1 (MQ=255)
  gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                                                                      >  1:1548490/1‑93 (MQ=255)
  gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                                                                      <  2:2054135/93‑1 (MQ=255)
  gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                                                                      >  2:4485666/1‑93 (MQ=255)
   ccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCca                                                                                     >  1:3967410/1‑93 (MQ=255)
    cgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcac                                                                                    >  2:2683886/1‑93 (MQ=255)
     gcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcaca                                                                                   >  1:410734/1‑93 (MQ=255)
      cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacac                                                                                  <  2:2369874/93‑1 (MQ=255)
      cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacac                                                                                  >  1:2611276/1‑93 (MQ=255)
      cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTAAGCGGAGACTCTTGCTCGCCCCATGCcacac                                                                                  >  2:1508165/1‑93 (MQ=255)
       aaaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacaca                                                                                 <  2:3896139/93‑1 (MQ=255)
       aaaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacaca                                                                                 >  1:3333272/1‑93 (MQ=255)
        aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCacacac                                                                                >  1:2392508/1‑93 (MQ=255)
          cGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGc                                                                              >  2:2786429/1‑93 (MQ=255)
           gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGGTCGCCCCCTGCCACACACGcc                                                                             <  1:841750/93‑1 (MQ=255)
           gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGcc                                                                             >  1:129383/1‑93 (MQ=255)
            ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                                                            <  1:1181333/93‑1 (MQ=255)
            ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                                                            <  2:2657988/93‑1 (MQ=255)
            ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                                                            >  1:4349286/1‑93 (MQ=255)
            ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTGCCTCCCCCCTTCCCCACACCCctg                                                                           >  1:1499171/1‑93 (MQ=255)
             gACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctg                                                                           <  2:838210/93‑1 (MQ=255)
              aCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACATGCctgc                                                                          <  2:4608636/93‑1 (MQ=255)
              aCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgc                                                                          >  1:590316/1‑93 (MQ=255)
              aCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgc                                                                          <  2:1874687/93‑1 (MQ=255)
               cTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgct                                                                         >  2:3758484/1‑93 (MQ=255)
               cTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgct                                                                         >  2:4554551/1‑93 (MQ=255)
                ttCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgctg                                                                        >  1:1442137/1‑93 (MQ=255)
                ttCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgctg                                                                        <  2:2699665/93‑1 (MQ=255)
                 tCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTgg                                                                       >  1:2785261/1‑93 (MQ=255)
                 tCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTgg                                                                       >  1:3416340/1‑93 (MQ=255)
                      gACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTaaa                                                                  >  2:792024/1‑93 (MQ=255)
                        cGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAaca                                                                <  2:4472898/93‑1 (MQ=255)
                         gAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacac                                                               <  2:3413814/93‑1 (MQ=255)
                           aCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacacca                                                             >  1:2773904/1‑93 (MQ=255)
                           aCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacacca                                                             <  1:4630659/93‑1 (MQ=255)
                            cTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacaccac                                                            <  1:291255/93‑1 (MQ=255)
                              tCAGTATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacaccacac                                                          <  1:2697290/93‑1 (MQ=255)
                                aGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAAcaccacacca                                                        >  1:3024759/1‑93 (MQ=255)
                                  gATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATc                                                      >  1:638689/1‑93 (MQ=255)
                                  gATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATc                                                      <  2:1065774/93‑1 (MQ=255)
                                  gATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATc                                                      <  2:505335/93‑1 (MQ=255)
                                   aTCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg                                                     >  2:2329672/1‑93 (MQ=255)
                                   aTCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg                                                     <  1:297933/93‑1 (MQ=255)
                                   aTCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg                                                     >  1:2941810/1‑93 (MQ=255)
                                   aTCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg                                                     >  1:10681/1‑93 (MQ=255)
                                      tcAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGc                                                  >  2:4113297/1‑93 (MQ=255)
                                      ccagccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGc                                                  <  2:1432107/92‑1 (MQ=255)
                                         gccagcGTCGTTCTTGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtc                                               >  2:4256840/1‑93 (MQ=255)
                                         gccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtc                                               >  1:3011430/1‑93 (MQ=255)
                                         gccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtc                                               >  1:4442214/1‑93 (MQ=255)
                                         gccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtc                                               <  1:1737659/93‑1 (MQ=255)
                                           cagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtctc                                             >  1:451717/1‑93 (MQ=255)
                                           cagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtctc                                             <  1:365736/93‑1 (MQ=255)
                                             gcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCgg                                           >  1:226909/1‑92 (MQ=255)
                                             gcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa                                           <  1:3216469/93‑1 (MQ=255)
                                             gcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa                                           <  1:3087371/93‑1 (MQ=255)
                                             gcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa                                           >  1:226910/1‑93 (MQ=255)
                                                  cgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGGAAAAACCACACCATCGTGCGCCTCGTTATg                                       >  2:1709358/2‑93 (MQ=255)
                                                     tCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGcac                                    <  2:1716819/93‑1 (MQ=255)
                                                       tGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACaa                                  <  1:1109983/93‑1 (MQ=255)
                                                       tGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACaa                                  >  2:1295631/1‑93 (MQ=255)
                                                       tGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACaa                                  >  1:2005551/1‑93 (MQ=255)
                                                         gAAATCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATg                                <  1:2495035/93‑1 (MQ=255)
                                                         gAAAGCCCCGTCAGCGGAGACTCTTGCTCCCCCCCTGCCACACACGCCTGCTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATg                                <  1:2495036/93‑1 (MQ=255)
                                                           aaGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGtt                              <  1:2598557/93‑1 (MQ=255)
                                                            aGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTg                             >  2:3447109/1‑93 (MQ=255)
                                                            aGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTg                             <  2:4558273/93‑1 (MQ=255)
                                                               cccGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                          <  1:2730570/92‑1 (MQ=255)
                                                               aCCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                          >  2:4308103/1‑93 (MQ=255)
                                                                ccGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAt                         >  1:3827753/1‑93 (MQ=255)
                                                                  gTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGc                       >  1:3943794/1‑93 (MQ=255)
                                                                  gTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGc                       >  1:2266408/1‑93 (MQ=255)
                                                                   tgaGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      <  1:538936/91‑1 (MQ=255)
                                                                   tCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      >  1:930869/1‑93 (MQ=255)
                                                                      gCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccgc                   >  1:682908/1‑93 (MQ=255)
                                                                       cGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccgcc                  >  2:1114664/1‑93 (MQ=255)
                                                                        ggAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGCGCGTCTCGATATGCACAATGTTAGTATGTccgccg                 >  2:3454570/1‑93 (MQ=255)
                                                                         gagaCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcga                >  1:89306/1‑93 (MQ=255)
                                                                         gagaCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcga                <  2:213988/93‑1 (MQ=255)
                                                                          agaCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgac               <  1:1992594/93‑1 (MQ=255)
                                                                            aCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacga             <  2:2437443/93‑1 (MQ=255)
                                                                            aCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacga             <  1:3246082/93‑1 (MQ=255)
                                                                            aCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacga             <  2:379994/93‑1 (MQ=255)
                                                                            aCTCGTGCTCGCCCCCTGCCACACACGCCTGCTGTGTAAACACCACACCATCGTGCGTCTCGATATGCACAATTTTGGTATGCCCGCcgacga             <  1:1984850/93‑1 (MQ=255)
                                                                               ctTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATgg          <  2:197909/93‑1 (MQ=255)
                                                                               ctTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATgg          <  1:4133183/93‑1 (MQ=255)
                                                                                ttGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATggg         >  1:1975727/1‑92 (MQ=255)
                                                                                ttGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGt         >  2:3142665/1‑93 (MQ=255)
                                                                                  gCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGa       <  1:1290709/93‑1 (MQ=255)
                                                                                   cTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGac      >  2:4344668/1‑93 (MQ=255)
                                                                                   cTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGac      <  2:2498957/93‑1 (MQ=255)
                                                                                    tCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGaca     >  1:2976577/1‑93 (MQ=255)
                                                                                    tCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGaca     >  2:726689/1‑93 (MQ=255)
                                                                                    tCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGaca     <  2:2907476/93‑1 (MQ=255)
                                                                                       ccccCTGGCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGacacac  <  2:3983838/93‑1 (MQ=255)
                                                                                       ccccCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGacacac  >  1:2384678/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACAC  >  NC_000913/3253691‑3253870

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