Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,304,838 A→G intergenic (+445/+270) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,8380AG98.1% 426.4 / ‑4.8 159intergenic (+445/+270)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (77/77);  total (81/78)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

GAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAT  >  NC_000913/2304741‑2304897
                                                                                                 |                                                           
gAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCggg                                                           >  1:1020172/1‑100 (MQ=255)
 aaGTGCTCCACCCCGTAGGTTGGATAAGGCTGTTACGCAGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGc                                                          <  2:3293126/100‑1 (MQ=255)
  ggggCTCCCCCCCGAAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGcc                                                         <  2:4530645/97‑1 (MQ=14)
  aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGcc                                                         >  1:2800690/1‑100 (MQ=255)
  aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGcc                                                         <  2:656797/100‑1 (MQ=255)
  aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGcc                                                         <  1:5133953/100‑1 (MQ=255)
  aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGcc                                                         >  2:4770420/1‑100 (MQ=255)
  aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGcc                                                         <  2:2895933/100‑1 (MQ=255)
   gTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCt                                                        <  2:2931437/100‑1 (MQ=255)
     ncTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTac                                                      >  1:88462/2‑100 (MQ=255)
     gCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCGACGGCAAAAGCGGCTGGCTGAGGCAACGCGGGCGCGCCTGATCAGGccggc                                                      >  2:4118825/1‑97 (MQ=2)
     gCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTac                                                      >  1:1010408/1‑100 (MQ=255)
     gCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTac                                                      <  1:3811621/100‑1 (MQ=255)
      cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTaca                                                     >  2:4596390/1‑100 (MQ=255)
      cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTaca                                                     >  2:765635/1‑100 (MQ=255)
      cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGCCCTaca                                                     >  2:4945609/1‑100 (MQ=255)
       tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTacac                                                    >  2:3132400/1‑100 (MQ=255)
       tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTacac                                                    >  1:5125290/1‑100 (MQ=255)
        ccACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTacacn                                                   <  1:4104610/100‑2 (MQ=255)
        ccACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACAcc                                                   >  2:25043/1‑100 (MQ=255)
         cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCg                                                  <  1:558404/100‑1 (MQ=255)
         cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCg                                                  >  2:2852832/1‑100 (MQ=255)
         cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCg                                                  >  1:311714/1‑100 (MQ=255)
           ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACCCCGCt                                                >  2:2263738/1‑100 (MQ=255)
           ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCt                                                >  2:2811870/1‑100 (MQ=255)
           ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCt                                                >  1:3331291/1‑100 (MQ=255)
            cccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtg                                               <  1:2601821/100‑1 (MQ=255)
            cccGCAGGTTGGACAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtg                                               <  2:4260204/100‑1 (MQ=14)
              cGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtgtg                                             <  1:1728040/100‑1 (MQ=255)
              cGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtgtg                                             <  1:4369213/100‑1 (MQ=255)
               gTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGa                                            <  1:316217/100‑1 (MQ=255)
                tAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGaa                                           >  1:3264132/1‑100 (MQ=255)
                 aGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAg                                          <  1:2649734/100‑1 (MQ=255)
                 aGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAg                                          <  2:2564104/100‑1 (MQ=255)
                 aGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAg                                          <  1:694081/100‑1 (MQ=255)
                 aGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAg                                          <  1:2013076/100‑1 (MQ=255)
                  ggTTGGATAAGGCGGTTACGTCGCATCCGACATCCAACGTCCGAGCCGGTTGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGt                                         >  2:5066302/1‑100 (MQ=255)
                    ttGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGc                                       <  2:3949319/100‑1 (MQ=255)
                     tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTAGACCGCTGTGGAGTGcc                                      >  2:371755/1‑99 (MQ=14)
                     tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCt                                      <  1:2190806/100‑1 (MQ=255)
                      ngATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGGAGttccc                                     >  1:4122422/2‑96 (MQ=14)
                      ggATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCCCTGTGAAGGGCtg                                     >  2:2726381/1‑99 (MQ=14)
                        aTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCa                                   <  1:3951495/100‑1 (MQ=255)
                         tAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTCCCCCGCTGTGAAGTGCCCCAc                                  >  2:3154505/1‑100 (MQ=14)
                         tAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAc                                  <  2:2002988/100‑1 (MQ=255)
                          aaGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAcc                                 >  1:4180357/1‑100 (MQ=255)
                           aGGCGGTTACGCCGCTTTCGATAACCCACGGCCGAGCGGGTTGCCTGAGGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAccc                                <  2:106991/100‑1 (MQ=11)
                           aGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAccc                                >  1:1290774/1‑100 (MQ=255)
                           aGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAccc                                >  2:74530/1‑100 (MQ=255)
                           aGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAccc                                <  2:2870351/100‑1 (MQ=255)
                           aGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAccc                                <  1:21903/100‑1 (MQ=255)
                           aGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGCGCCCCAccc                                >  2:3416235/1‑100 (MQ=255)
                            ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAcccc                               >  1:1918347/1‑100 (MQ=255)
                            ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAcccc                               <  2:1284031/100‑1 (MQ=255)
                            ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAcccc                               >  1:705127/1‑100 (MQ=255)
                            ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAcccc                               >  1:253915/1‑100 (MQ=255)
                              cGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTCCACCGCTGTGACGTGCCCCACCCCGt                             >  2:3257163/1‑100 (MQ=11)
                              cGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGt                             >  2:1896255/1‑100 (MQ=37)
                              cGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGt                             <  1:1962054/100‑1 (MQ=255)
                              cGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGt                             >  2:2101967/1‑100 (MQ=255)
                              cGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGt                             <  1:1037164/100‑1 (MQ=255)
                               ggTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTa                            <  2:109629/100‑1 (MQ=255)
                               ggTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTCGCGCGTCTTATCGGGCCTACACCGCCGTGAAGTGCTCCACCCCGtc                            >  2:3649294/1‑99 (MQ=11)
                                 ttACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                          <  1:4983675/100‑1 (MQ=25)
                                    cgccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                       <  2:4657475/100‑1 (MQ=21)
                                    cgccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGTGCTCCACCCCGTAGGTCg                       >  2:2341147/1‑100 (MQ=18)
                                    cgccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTCATCGGGCCTCCACCGCTGCGAAGTGCTCCACCCCGTAGGTCg                       >  1:1906934/1‑100 (MQ=11)
                                    cgccgcATCCGACAGCCAACGTCCGAGCCGGTTGCCTGATGCGATGCTCGCGCGTCTTATCGGGCCTACACCCCTGTGAAGTGCTCCCCCCCGTAGGTCg                       >  1:4524635/1‑100 (MQ=11)
                                     gccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                      >  2:3775732/1‑100 (MQ=21)
                                      ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                     >  1:2079158/1‑100 (MQ=21)
                                      ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                     >  1:2458707/1‑100 (MQ=21)
                                      ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                     <  1:3506340/100‑1 (MQ=18)
                                        gcATCCAACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAACGTGCCTACGCCGCTGTGAAGTGCTCCATCCCGTAGGTCGCATa                   >  2:4842680/1‑100 (MQ=11)
                                          aTCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAg                 <  1:5055766/100‑1 (MQ=18)
                                            ccGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc               >  1:4774099/1‑100 (MQ=18)
                                            ccGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc               <  1:469541/100‑1 (MQ=18)
                                            ccGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc               <  1:4649770/100‑1 (MQ=21)
                                              gACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCGTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:2849676/100‑1 (MQ=255)
                                              gACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGTAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:4141249/100‑1 (MQ=255)
                                              gACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:3517036/100‑1 (MQ=18)
                                              gACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:2626350/100‑1 (MQ=18)
                                                cATCCAGCGTCCGAGCCGATTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt           >  1:4299734/1‑100 (MQ=11)
                                                 aTCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa          >  1:1526189/1‑100 (MQ=21)
                                                  tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc         <  1:1725533/100‑1 (MQ=18)
                                                  tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc         >  1:2829732/1‑100 (MQ=21)
                                                  tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGGTAc         >  2:2213315/1‑100 (MQ=255)
                                                  tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTACTCCACCCCGTAGGTCGGATAAGGCGGTTAc         >  2:2391094/1‑100 (MQ=255)
                                                  tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCCCTGTGAAGTACACCACCCCGTAGGTCGGATAAGGCGGTTAc         >  1:1794656/1‑100 (MQ=11)
                                                  tCCAACGTCCGAGCCGGTTGCCTGAAGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTATGTCGGATAAGGCGGTTAc         >  1:2352350/1‑100 (MQ=255)
                                                   ccAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGACGTGCTCCACCCCGTAGGTCGGATAAGACGGTTAcg        >  2:2618608/1‑100 (MQ=255)
                                                   ccAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcg        >  1:309526/1‑100 (MQ=18)
                                                    cAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgc       >  2:3608810/1‑100 (MQ=18)
                                                    cAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgc       >  2:991048/1‑100 (MQ=21)
                                                     aaCGTCCGTGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgcc      <  2:1557725/100‑1 (MQ=16)
                                                     aaCGTCCGAGCCGGTTGCCTGATGCGACGCTGTCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgcc      <  2:5083732/100‑1 (MQ=16)
                                                     aaCGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgcc      <  1:4352557/100‑1 (MQ=18)
                                                      aCGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg     <  2:1578983/100‑1 (MQ=18)
                                                      aCGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg     <  1:4058782/100‑1 (MQ=18)
                                                       cGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccgc    >  2:3313992/1‑100 (MQ=18)
                                                         tCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAt  <  1:4196971/100‑1 (MQ=14)
                                                                                                 |                                                           
GAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAT  >  NC_000913/2304741‑2304897

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