Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,283,932 A→G N110S (AAC→AGC)  narH → nitrate reductase A subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,283,9320AG88.1% 978.6 / 99.4 310N110S (AAC→AGC) narHnitrate reductase A subunit beta
Reads supporting (aligned to +/- strand):  ref base A (17/20);  new base G (123/150);  total (140/170)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

CATGGGTAACCGTGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAGGGCCGA  >  NC_000913/1283834‑1284030
                                                                                                  |                                                                                                  
 aTGGGTAACCGTGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCt                                                                                                 >  2:2359137/1‑101 (MQ=255)
      tAACCGTGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATa                                                                                            <  2:2358801/101‑1 (MQ=255)
      tAACCGTGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATa                                                                                            <  1:2358640/101‑1 (MQ=255)
       aaCCGTGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATAc                                                                                           <  1:2358509/101‑1 (MQ=255)
       aaCCGTACCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATAc                                                                                           >  2:2359144/1‑101 (MQ=255)
          cGTGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACcgc                                                                                        <  2:2358857/101‑1 (MQ=255)
           gTGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACcgcg                                                                                       <  1:2358847/101‑1 (MQ=255)
           gTGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACcgcg                                                                                       >  2:2359149/1‑101 (MQ=255)
            tGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCgcgc                                                                                      >  2:2359152/1‑101 (MQ=255)
              ccATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCg                                                                                    <  1:2358713/101‑1 (MQ=255)
               cATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCgg                                                                                   >  2:2359157/1‑101 (MQ=255)
                aTGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGa                                                                                  >  2:2359160/1‑101 (MQ=255)
                 tgctgcTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGaa                                                                                 >  1:2359162/1‑101 (MQ=255)
                 tgctgcTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGaa                                                                                 <  1:2358751/101‑1 (MQ=255)
                  gctgctGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAg                                                                                <  1:2358761/101‑1 (MQ=255)
                     gctgGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGgca                                                                             <  2:2358568/101‑1 (MQ=255)
                     gctgGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGgca                                                                             >  2:2359167/1‑101 (MQ=255)
                     gctgGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACAATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGaa                                                                                 <  1:2358652/97‑1 (MQ=255)
                      ctgGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGgcag                                                                            <  1:2358592/101‑1 (MQ=255)
                       tgGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGgcagc                                                                           >  1:2359168/1‑101 (MQ=255)
                       tgGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGgcagc                                                                           <  2:2358672/101‑1 (MQ=255)
                        gggTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGgcag                                                                            <  2:2358872/99‑1 (MQ=255)
                          gTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACAATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCaaa                                                                        <  2:2358557/101‑1 (MQ=255)
                           tAAAACCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAAt                                                                       <  2:2358719/101‑1 (MQ=255)
                            aaaaTCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATc                                                                      >  2:2359174/1‑101 (MQ=255)
                             aaaTCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGAAGCAAATCg                                                                     >  1:2359176/1‑101 (MQ=255)
                              aaTCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGc                                                                    >  2:2359178/1‑101 (MQ=255)
                                 cTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGc                                                                 >  2:2359182/1‑101 (MQ=255)
                                 cTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGc                                                                 <  2:2358420/101‑1 (MQ=255)
                                  ttCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGcc                                                                >  1:2359185/1‑101 (MQ=255)
                                   tCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCg                                                               <  1:2358736/101‑1 (MQ=255)
                                    cGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGa                                                              >  2:2359187/1‑101 (MQ=255)
                                        aaCCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGc                                                          >  1:2359191/1‑101 (MQ=255)
                                        aaCCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGc                                                          <  2:2358753/101‑1 (MQ=255)
                                        aaCCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGc                                                          <  1:2358468/101‑1 (MQ=255)
                                        aaCCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGc                                                          <  1:2358610/101‑1 (MQ=255)
                                         aCCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGcc                                                         >  2:2359194/1‑101 (MQ=255)
                                          cccGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGccc                                                        <  2:2359061/101‑1 (MQ=255)
                                             gCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTc                                                     <  1:2358520/101‑1 (MQ=255)
                                              cATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTcc                                                    <  1:2359149/101‑1 (MQ=255)
                                              cATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTcc                                                    <  1:2358576/101‑1 (MQ=255)
                                                tCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTc                                                     >  2:2359198/1‑98 (MQ=255)
                                                 cTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCcgcg                                                 <  2:2358694/101‑1 (MQ=255)
                                                  tGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGgcag                                                                            >  2:2359203/1‑73 (MQ=255)
                                                   gCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTcc                                                    >  1:2359204/1‑96 (MQ=255)
                                                     cGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCg                                             >  1:2359209/1‑101 (MQ=255)
                                                     cGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCg                                             >  2:2359211/1‑101 (MQ=255)
                                                      ggggATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCcgc                                                  <  1:2358663/95‑1 (MQ=255)
                                                        ggATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTg                                          <  2:2358599/101‑1 (MQ=255)
                                                         gATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGa                                         <  1:2358731/101‑1 (MQ=255)
                                                          aTCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTg                                          <  1:2359143/99‑1 (MQ=255)
                                                          aTCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGAt                                        <  1:2358369/101‑1 (MQ=255)
                                                           tCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGAtt                                       <  1:2358634/101‑1 (MQ=255)
                                                            cgacgaTTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTa                                      >  1:2359215/1‑101 (MQ=255)
                                                             gacgaTTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTAc                                     >  2:2359219/1‑101 (MQ=255)
                                                              acgaTTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCa                                                                   <  2:2358778/70‑1 (MQ=255)
                                                               cgaTTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCg                                   >  2:2359222/1‑101 (MQ=255)
                                                                gaTTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCgg                                  >  1:2359223/1‑101 (MQ=255)
                                                                gaTTATTACGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCgg                                  >  1:2359226/1‑101 (MQ=255)
                                                                  ttattaCGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCgggg                                >  1:2359230/1‑101 (MQ=255)
                                                                    attcCGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTc                                              <  2:2358644/85‑1 (MQ=255)
                                                                     ttaCGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCcg                                                   <  2:2358508/79‑1 (MQ=255)
                                                                     ttaCGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAAc                             <  1:2358676/101‑1 (MQ=255)
                                                                      taCGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGcc                                                                >  1:2359234/1‑65 (MQ=255)
                                                                       aCGAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACgg                           <  1:2358546/101‑1 (MQ=255)
                                                                        cGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGa                          >  2:2359236/1‑101 (MQ=255)
                                                                         gAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGaa                              <  2:2358582/96‑1 (MQ=255)
                                                                         gAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGAt                         <  1:2358679/101‑1 (MQ=255)
                                                                         gAGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGAt                         >  2:2359239/1‑101 (MQ=255)
                                                                          aGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAAt                                                                       >  2:2359243/1‑54 (MQ=255)
                                                                          aGCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATg                        <  2:2358925/101‑1 (MQ=255)
                                                                           gCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATgg                       <  1:2358693/101‑1 (MQ=255)
                                                                           gCCGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATgg                       >  2:2359246/1‑101 (MQ=255)
                                                                             cGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGAtt                                       <  2:2358825/83‑1 (MQ=255)
                                                                             cGTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCg                     >  2:2359250/1‑101 (MQ=255)
                                                                              gTTCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGa                    >  1:2359252/1‑101 (MQ=255)
                                                                               ttCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGaa                   <  1:2358934/101‑1 (MQ=255)
                                                                                tCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGaaa                  <  2:2358993/101‑1 (MQ=255)
                                                                                tCGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGaaa                  <  1:2358704/101‑1 (MQ=255)
                                                                                 cGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGaaaa                 <  2:2358688/101‑1 (MQ=255)
                                                                                 cGATTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGaaaa                 >  2:2359255/1‑101 (MQ=255)
                                                                                  gATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAAt                <  1:2358836/101‑1 (MQ=255)
                                                                                   aTTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATc               <  1:2358912/101‑1 (MQ=255)
                                                                                   aTTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATc               >  1:2359258/1‑101 (MQ=255)
                                                                                   aTTTTGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATc               <  2:2358618/101‑1 (MQ=255)
                                                                                    ttttGACTATCAGAGCCTGCATCCCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCt                                           <  2:2358528/72‑1 (MQ=255)
                                                                                    ttttGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCg              <  2:2358741/101‑1 (MQ=255)
                                                                                     tttGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGa             >  2:2359261/1‑101 (MQ=255)
                                                                                      ttGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGaa            >  2:2359264/1‑101 (MQ=255)
                                                                                       tGACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGaaa           <  1:2359030/101‑1 (MQ=255)
                                                                                        gACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGaaaa          <  1:2358776/101‑1 (MQ=255)
                                                                                        gACTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGaaaa          <  2:2358635/101‑1 (MQ=255)
                                                                                         aCTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGaaaaa         <  1:2358886/101‑1 (MQ=255)
                                                                                          cTATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAg        <  1:2358815/101‑1 (MQ=255)
                                                                                           tATCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAgg       >  2:2359268/1‑101 (MQ=255)
                                                                                            aTCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAggg      >  1:2359272/1‑101 (MQ=255)
                                                                                            aTCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAggg      >  2:2359271/1‑101 (MQ=255)
                                                                                             tCAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAGGGc     <  2:2358733/101‑1 (MQ=255)
                                                                                              cAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAGGGcc    >  1:2359277/1‑101 (MQ=255)
                                                                                              cAGAGCCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAGGGcc    >  2:2359278/1‑101 (MQ=255)
                                                                                                  |                                                                                                  
CATGGGTAACCGTGCCATGCTGCTGGGTAAAATCTTCGCTAACCCGCATCTGCCGGGGATCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAGGGCCGA  >  NC_000913/1283834‑1284030

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