Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,552,310 T→G intergenic (+319/+88) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,3100TG96.0% 305.1 / 0.8 102intergenic (+319/+88)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (51/46);  total (54/48)
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.80e-01

Only 100 of 102 total aligned reads displayed.

CTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGT  >  NC_000913/1552269‑1552406
                                         |                                                                                                
cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTc                                        <  1:4936751/100‑1 (MQ=14)
 ttATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCt                                       <  1:3233073/100‑1 (MQ=17)
  tATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTa                                      <  1:3445070/100‑1 (MQ=18)
    tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTaaa                                    >  1:1588620/1‑100 (MQ=25)
    tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTaaa                                    >  1:306649/1‑100 (MQ=25)
    tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTaaa                                    >  1:2503661/1‑100 (MQ=25)
    tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTaaa                                    >  1:2503792/1‑100 (MQ=25)
     ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAAt                                   >  2:4503930/1‑100 (MQ=255)
     ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAAt                                   <  1:4830238/100‑1 (MQ=255)
     ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAAt                                   <  2:3444995/100‑1 (MQ=255)
     ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAAt                                   >  1:2700652/1‑100 (MQ=255)
     acGGCATGAAACCCCCTCAAATCCCATGGCTTGCTGGGGATGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAAt                                   <  2:4077802/99‑1 (MQ=11)
      cGGCATTACCACTCTTCAAATCCACTAGATTTCGGGGTACGTCAGGGCCTGATAAGCGTAGCGCATCAGGCATGGTGGCGTTTGTCATCAGTTCTAAAtt                                  >  2:4319439/1‑99 (MQ=2)
      cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                  >  2:2256557/1‑100 (MQ=255)
      cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                  >  2:3075817/1‑100 (MQ=255)
      cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                  <  1:3732692/100‑1 (MQ=255)
      cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                  >  2:3335185/1‑100 (MQ=255)
      cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                  <  2:898352/100‑1 (MQ=255)
      cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATg                                  >  1:232975/1‑100 (MQ=255)
      cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTCTGTCATCAGTTCTAAATg                                  >  2:4664843/1‑100 (MQ=255)
       ggCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATgg                                 <  2:1778225/100‑1 (MQ=37)
       ggCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATgg                                 >  1:5263590/1‑100 (MQ=255)
       ggCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATgg                                 >  2:1684282/1‑100 (MQ=255)
       ggCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATgg                                 <  2:1915187/100‑1 (MQ=255)
       ggCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATgg                                 <  2:2138937/100‑1 (MQ=255)
       ggCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATgg                                 >  1:1253454/1‑100 (MQ=255)
        gCATGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGgc                                >  2:2667705/1‑100 (MQ=255)
          aTGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGgcgc                              <  1:4720471/100‑1 (MQ=37)
           tGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCt                             <  1:2793539/100‑1 (MQ=38)
           tGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCt                             <  1:3489604/100‑1 (MQ=255)
           tGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCt                             >  2:1214163/1‑100 (MQ=255)
           tGAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCt                             >  1:3901789/1‑100 (MQ=255)
            gAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCtt                            >  2:17429/1‑100 (MQ=255)
            gAAAACCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCtg                            >  2:1035182/1‑99 (MQ=25)
             aaaaCCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCGATGTTGCGTTTGTCATCAGTTCTAAATGGCACttt                           >  2:4637238/1‑100 (MQ=14)
             aaaaCCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCttt                           >  2:3472290/1‑100 (MQ=255)
             aaaaCCCTTCAAATCCAATAGATGGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCttt                           <  2:2997973/100‑1 (MQ=255)
              aaaCCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTta                          >  1:1921724/1‑100 (MQ=255)
              aaaCCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTta                          <  2:447005/100‑1 (MQ=255)
              aaaCCCTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTta                          <  1:1332911/100‑1 (MQ=255)
                 cccTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATaa                       <  1:1728401/100‑1 (MQ=255)
                 cccTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATaa                       >  1:1549953/1‑100 (MQ=255)
                 cccTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATaa                       >  1:2733581/1‑100 (MQ=255)
                 cccTTCAAATCCAATAGATTGCAGGGAACGTGTAGGACTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTATAAATGGCGCTTTCTaa                       >  1:4779363/1‑100 (MQ=11)
                   cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGACAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTTTATAAAg                     >  2:4229366/1‑100 (MQ=255)
                   cTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAg                     >  1:1032831/1‑100 (MQ=255)
                   cTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAg                     >  2:3685338/1‑100 (MQ=255)
                   cTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAg                     <  1:1909784/100‑1 (MQ=255)
                   cTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAg                     >  1:4633436/1‑100 (MQ=255)
                   cTTCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAg                     >  1:2535814/1‑100 (MQ=255)
                   cTTCAAATCCAATAGATTGCAGGGAACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGAAATAAGTTCTAAATGGCGCTTTATAAAg                     >  2:1128784/1‑100 (MQ=14)
                    ttCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTGTAAATGGCGCTTTATAAAGt                    <  1:4354618/100‑1 (MQ=255)
                    ttCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                    <  1:1143979/100‑1 (MQ=255)
                    ttCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                    <  1:3737155/100‑1 (MQ=255)
                    ttCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                    >  1:2800758/1‑100 (MQ=255)
                    ttCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                    >  1:4579271/1‑100 (MQ=255)
                    ttCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGATTTATAAAGt                    <  1:4354617/100‑1 (MQ=255)
                     tCAAATCGAATAGCTGGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTg                   <  2:3091372/100‑1 (MQ=17)
                     tCAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTg                   <  2:2753308/100‑1 (MQ=255)
                     tCAAAGCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGtn                   <  2:2577087/100‑2 (MQ=255)
                      cAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGc                  >  2:4161823/1‑100 (MQ=255)
                      cAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGc                  <  1:2126349/100‑1 (MQ=255)
                      cAAATCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGc                  >  1:2870995/1‑100 (MQ=255)
                       aaaTCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGcc                 <  1:2046601/100‑1 (MQ=255)
                        aaTCCAATAGATTGCAGTGAACGTGTAGGCCGGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTTAAATGGCGCTTTATAAAggtcgg                >  1:3934550/1‑95 (MQ=14)
                        aaTCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCa                >  2:3499028/1‑100 (MQ=255)
                         aTCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCCTTATAAGGTGCCAt               >  1:3899436/1‑100 (MQ=17)
                         aTCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAt               <  1:4696279/100‑1 (MQ=17)
                          tCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTGTAAAGTGCCAtt              >  2:3426888/1‑100 (MQ=255)
                          tCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtt              <  2:1028682/100‑1 (MQ=255)
                          tCCAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtt              >  2:1570399/1‑100 (MQ=255)
                           ccAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttt             >  2:4340214/1‑100 (MQ=255)
                            cAATAGATTGCAGGGACCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt            >  2:3940009/1‑100 (MQ=255)
                            cAATAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAttta            <  2:3953855/100‑2 (MQ=17)
                             aaTCGATTGAAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt           <  2:2929603/100‑1 (MQ=255)
                             aaTAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt           >  2:707501/1‑100 (MQ=255)
                             aaTAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttag           <  1:1136406/100‑3 (MQ=17)
                              aTAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt          <  2:1675513/100‑1 (MQ=255)
                              aTAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt          <  2:3756780/100‑1 (MQ=255)
                              aTAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt          >  2:2661393/1‑100 (MQ=255)
                              aTAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt          <  2:303770/100‑1 (MQ=255)
                              aTAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt          <  2:4900378/100‑1 (MQ=255)
                              aTAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttaga          <  2:795334/100‑4 (MQ=17)
                              aTAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttaga          <  2:2410166/100‑4 (MQ=17)
                               tAGATTGCAGGGAACGTGTAGGCCTGATAATCGTAGCGCATCAGGGAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCActgtttt         <  2:5177309/100‑1 (MQ=17)
                               tAGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttt         <  2:4019413/100‑1 (MQ=255)
                                aGATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa        >  1:5446049/1‑100 (MQ=255)
                                 gATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt       >  2:619387/1‑100 (MQ=255)
                                 gATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt       >  2:1268672/1‑100 (MQ=255)
                                 gATTGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt       <  2:2841870/100‑1 (MQ=255)
                                   ttGCAGGGAACGTGTNGGCCTGATAANCGNAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATNg     <  2:2141220/100‑1 (MQ=255)
                                   ttGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTg     <  1:1509324/100‑1 (MQ=255)
                                    tGCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGc    >  1:42094/1‑100 (MQ=255)
                                     gCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCGTCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg   >  1:2216719/1‑100 (MQ=255)
                                     gCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg   <  2:1847004/100‑1 (MQ=255)
                                     gCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg   <  1:4032853/100‑1 (MQ=255)
                                     gCAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg   <  1:1014947/100‑1 (MQ=255)
                                      cAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGt  >  1:4870202/1‑100 (MQ=255)
                                      cAGGGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGt  >  1:3602862/1‑100 (MQ=255)
                                      cAGCGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAAGGTTGCGCTTGTCATCAGTTCTAAACGGCGCTTTATAAAGTGCCATTTTTTTATTGcct  >  2:1495221/1‑98 (MQ=18)
                                         |                                                                                                
CTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGT  >  NC_000913/1552269‑1552406

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: