Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,289,526 A→G K94K (AAA→AAG nrfB → nitrite reductase, formate‑dependent, penta‑heme cytochrome c

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,289,5260AG98.0% 501.2 / ‑2.3 151K94K (AAA→AAGnrfBnitrite reductase, formate‑dependent, penta‑heme cytochrome c
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (77/68);  total (81/70)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 151 total aligned reads displayed.

AAACATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTT  >  NC_000913/4289431‑4289621
                                                                                               |                                                                                                
aaaCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtg                                                                                              <  1:3274810/100‑1 (MQ=255)
 aaCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtgt                                                                                            >  2:2531312/1‑100 (MQ=255)
  aCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtgtg                                                                                           >  1:943320/1‑100 (MQ=255)
  aCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtgtg                                                                                           <  1:2053670/100‑1 (MQ=255)
   cctGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGa                                                                                          <  1:2451755/98‑1 (MQ=255)
   cATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGa                                                                                          >  1:763380/1‑100 (MQ=255)
   cATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGa                                                                                          >  1:1264004/1‑100 (MQ=255)
   cATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCACGAAGGGATGAAGGaagaaa                                                                                           >  2:911993/1‑95 (MQ=255)
    aTGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGAt                                                                                         <  2:943332/100‑1 (MQ=255)
    aTGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGAt                                                                                         <  1:1428186/100‑1 (MQ=255)
      gCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATgc                                                                                       >  1:311049/1‑100 (MQ=255)
       cTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATgcg                                                                                      >  2:762603/1‑100 (MQ=255)
          ccGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGCTGCGCtt                                                                                   >  1:1872057/1‑100 (MQ=255)
          ccGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCtt                                                                                   <  1:657429/100‑1 (MQ=255)
          ccGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCtt                                                                                   >  1:2911306/1‑100 (MQ=255)
           cGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCttt                                                                                  <  2:2643024/100‑1 (MQ=255)
           cGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCttt                                                                                  >  1:1314590/1‑100 (MQ=255)
              tATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAAc                                                                               >  2:3409140/1‑100 (MQ=255)
              tATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAAc                                                                               <  2:752542/100‑1 (MQ=255)
              tATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAAc                                                                               <  2:1527906/100‑1 (MQ=255)
               aTCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACg                                                                              <  1:2728012/100‑1 (MQ=255)
               aTCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACg                                                                              >  1:2631141/1‑100 (MQ=255)
                  aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGACGGGGTGAAGGATGTGATGCGCTTTAACGAGc                                                                           >  1:2166458/1‑100 (MQ=255)
                  aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGc                                                                           >  1:293854/1‑100 (MQ=255)
                  aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGc                                                                           <  2:1931792/100‑1 (MQ=255)
                  aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGc                                                                           >  1:1067848/1‑100 (MQ=255)
                   aCCCGAATAACAAACTGCCTGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGcc                                                                          <  1:3140207/100‑1 (MQ=255)
                    cccGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCg                                                                         <  2:3354880/100‑1 (MQ=255)
                    cccGAATAACAAACTGCCCGTCACCTGCACCAACTCTCCCCGCCACCCATCACCGCACCACCCCGTCTGCGGGAACGATCTGATCCGCTTTAACCACCCg                                                                         >  1:3021649/1‑100 (MQ=255)
                     ccGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCg                                                                         <  1:3286334/99‑1 (MQ=255)
                      cGAATAACAAACTACCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGAt                                                                       >  1:2153164/1‑100 (MQ=255)
                       gAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATg                                                                      >  2:1867409/1‑100 (MQ=255)
                        aaTAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGt                                                                     >  1:1961783/1‑100 (MQ=255)
                          tAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTAc                                                                   >  1:2690840/1‑100 (MQ=255)
                           accaacagcccGGACACCTGCACCAATGGACAGGGCCGGCCATCACGCCCACTCCGCGAAGGGGTACAGGTGTTGATGCGCTTTAACGAGCCGATGTACa                                                                  <  1:2085876/91‑1 (MQ=255)
                           aacaacCTCCCGGACACCTCCACAAATTGGGTCGGCCAGCACTCACCCCAACACCGCGAAGGGGTGAAGGGTGTGATGCGCTTTAACGAGCCGATGAACa                                                                  <  1:318486/100‑1 (MQ=255)
                           aaCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACa                                                                  <  2:2697652/100‑1 (MQ=255)
                            aCAAACTGCCGGTCACCTGCACCACCTGTCACCGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACa                                                                  <  2:735452/99‑1 (MQ=255)
                            aCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACaa                                                                 >  1:1059554/1‑100 (MQ=255)
                            aCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACaa                                                                 >  2:152675/1‑100 (MQ=255)
                            aCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACaa                                                                 >  1:2101209/1‑100 (MQ=255)
                             cAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAg                                                                <  2:1059568/100‑1 (MQ=255)
                             cAAACTGACGGACCCCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAg                                                                >  1:1820953/1‑100 (MQ=255)
                               aaCTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGt                                                              <  2:1431894/100‑1 (MQ=255)
                                 cTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTg                                                            >  2:1277985/1‑100 (MQ=255)
                                 cTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTg                                                            <  2:3231170/100‑1 (MQ=255)
                                   gCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTgg                                                           <  1:1009171/99‑1 (MQ=255)
                                   gCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGga                                                          >  1:1909413/1‑100 (MQ=255)
                                   gCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGga                                                          >  2:1108604/1‑100 (MQ=255)
                                    ccGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGgag                                                         <  1:382111/100‑1 (MQ=255)
                                     cGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGgaga                                                        >  1:593775/1‑100 (MQ=255)
                                      ggTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGagag                                                       <  1:2963913/100‑1 (MQ=255)
                                       gTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGc                                                      <  1:1152456/100‑1 (MQ=255)
                                         cACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAg                                                    >  2:3246549/1‑100 (MQ=255)
                                         cACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAg                                                    >  2:636643/1‑100 (MQ=255)
                                          aCCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGa                                                   >  1:830797/1‑100 (MQ=255)
                                           ccTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGaa                                                  >  2:2433302/1‑100 (MQ=255)
                                           ccTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGaa                                                  >  2:2188020/1‑100 (MQ=255)
                                             tGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACa                                                >  2:516315/1‑100 (MQ=255)
                                              gCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGCAGGGGTGGAGGATGTGATGCGCTTTACCGAGCCGATGTACAAGGTTGGAGAGCAGAACAg                                               >  1:1556343/1‑100 (MQ=255)
                                              gCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAg                                               >  1:2540451/1‑100 (MQ=255)
                                                aCCAACTGTCACGGCCAGCCATCACCGCAACACCGCGCAGGGGGGAAGGATGTGATGCGCTGTAACGAGCCGCCGTACAAGGTTGGAGAGCAGAACAGCg                                             >  1:1229032/1‑100 (MQ=255)
                                                aCCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCg                                             >  2:2688826/1‑100 (MQ=255)
                                                      tGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTa                                       <  1:746351/100‑1 (MQ=255)
                                                      tGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTa                                       <  2:2607972/100‑1 (MQ=255)
                                                       gTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGGGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTa                                       <  1:536696/99‑1 (MQ=255)
                                                        tCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATg                                     <  2:1872084/100‑1 (MQ=255)
                                                        tCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATg                                     >  2:2372330/1‑100 (MQ=255)
                                                         cACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGt                                    >  1:1464640/1‑100 (MQ=255)
                                                          aCGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTc                                   <  1:2245275/100‑1 (MQ=255)
                                                           cGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGGGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTcc                                  <  1:3050454/100‑1 (MQ=255)
                                                           cGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTcc                                  >  2:1670092/1‑100 (MQ=255)
                                                           cGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTcc                                  >  1:3218480/1‑100 (MQ=255)
                                                            gtccagccaTCACAGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCt                                 >  2:2356398/3‑100 (MQ=255)
                                                             gccagccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTg                                >  1:1456247/1‑100 (MQ=255)
                                                              ccagccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGt                               <  2:1215459/100‑1 (MQ=255)
                                                                agccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAGCAGCGTCTGTATGTCCTGTCa                             <  2:2944468/100‑1 (MQ=255)
                                                                agccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCa                             >  2:3010985/1‑100 (MQ=255)
                                                                agccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCa                             <  2:2153195/100‑1 (MQ=255)
                                                                 gccaTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCAt                            >  2:2398198/1‑100 (MQ=255)
                                                                   caTCACCGCACCACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGTGCAGAACAGCGTCTGTATGTCCTGTCATCt                          <  1:2279606/100‑1 (MQ=255)
                                                                   caTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTACGTCCTGTCATCt                          <  2:2282378/100‑1 (MQ=255)
                                                                    aTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTg                         >  2:2540890/1‑100 (MQ=255)
                                                                    aTCACCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTg                         >  2:1469279/1‑100 (MQ=255)
                                                                       aCCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCg                      <  2:2504277/100‑1 (MQ=255)
                                                                       aCCGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCg                      <  1:1566374/100‑1 (MQ=255)
                                                                        ccGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCgg                     <  1:3409088/100‑1 (MQ=255)
                                                                         cGCAACACCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGa                    >  1:2584910/1‑100 (MQ=255)
                                                                           cAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAAc                  <  1:1807216/100‑1 (MQ=255)
                                                                             acacCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAg                <  2:2653354/100‑1 (MQ=255)
                                                                              cacCGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGt               >  2:1717813/1‑100 (MQ=255)
                                                                                ctgcgAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTg             <  2:9383/98‑1 (MQ=255)
                                                                                ccGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTg             >  1:639443/1‑100 (MQ=255)
                                                                                ccGCGAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTg             <  2:1056745/100‑1 (MQ=255)
                                                                                  gcgAAGGGGTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTgag           >  1:9383/1‑98 (MQ=255)
                                                                                       ggggTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAg      <  2:1901065/100‑1 (MQ=255)
                                                                                       ggggTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAg      <  2:2259036/100‑1 (MQ=255)
                                                                                        gggTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGc     >  1:2396619/1‑100 (MQ=255)
                                                                                        gggTGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGc     <  2:125407/100‑1 (MQ=255)
                                                                                           tGAAGGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGtt  >  2:808585/1‑100 (MQ=255)
                                                                                               |                                                                                                
AAACATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATGTGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTT  >  NC_000913/4289431‑4289621

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