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→AAG) nrfBnitrite 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.

AAACATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTT  >  NC_000913/4289431‑4289621
                                                                                               |                                                                                                
aaaCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtg                                                                                              <  1:3274810/100‑1 (MQ=255)
 aaCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtg‑t                                                                                            >  2:2531312/1‑100 (MQ=255)
  aCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtg‑tg                                                                                           >  1:943320/1‑100 (MQ=255)
  aCATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGAtg‑tg                                                                                           <  1:2053670/100‑1 (MQ=255)
   cctGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGa                                                                                          <  1:2451755/98‑1 (MQ=255)
   cATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGa                                                                                          >  1:763380/1‑100 (MQ=255)
   cATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGa                                                                                          >  1:1264004/1‑100 (MQ=255)
   cATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCACGAAGGGATGAAGGaagaaa                                                                                           >  2:911993/1‑95 (MQ=255)
    aTGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGAt                                                                                         <  2:943332/100‑1 (MQ=255)
    aTGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGAt                                                                                         <  1:1428186/100‑1 (MQ=255)
      gCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATgc                                                                                       >  1:311049/1‑100 (MQ=255)
       cTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATgcg                                                                                      >  2:762603/1‑100 (MQ=255)
          ccGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGCTGCGCtt                                                                                   >  1:1872057/1‑100 (MQ=255)
          ccGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCtt                                                                                   <  1:657429/100‑1 (MQ=255)
          ccGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCtt                                                                                   >  1:2911306/1‑100 (MQ=255)
           cGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCttt                                                                                  <  2:2643024/100‑1 (MQ=255)
           cGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCttt                                                                                  >  1:1314590/1‑100 (MQ=255)
              tATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAAc                                                                               >  2:3409140/1‑100 (MQ=255)
              tATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAAc                                                                               <  2:752542/100‑1 (MQ=255)
              tATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAAc                                                                               <  2:1527906/100‑1 (MQ=255)
               aTCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACg                                                                              <  1:2728012/100‑1 (MQ=255)
               aTCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACg                                                                              >  1:2631141/1‑100 (MQ=255)
                  aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGACGGGGTGAAGGATG‑TGATGCGCTTTAACGAGc                                                                           >  1:2166458/1‑100 (MQ=255)
                  aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGc                                                                           >  1:293854/1‑100 (MQ=255)
                  aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGc                                                                           <  2:1931792/100‑1 (MQ=255)
                  aaCCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGc                                                                           >  1:1067848/1‑100 (MQ=255)
                   aCCCGAATAACAAACTGCCTGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGcc                                                                          <  1:3140207/100‑1 (MQ=255)
                    cccGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCg                                                                         <  2:3354880/100‑1 (MQ=255)
                    cccGAATAACAAACTGCCCGTCACCTGCACCAACTCTCCCCGCCACCCATCACCGCACCACCCCGTCTGCGGGAACGATC‑TGATCCGCTTTAACCACCCg                                                                         >  1:3021649/1‑100 (MQ=255)
                     ccGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCg                                                                         <  1:3286334/99‑1 (MQ=255)
                      cGAATAACAAACTACCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGAt                                                                       >  1:2153164/1‑100 (MQ=255)
                       gAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATg                                                                      >  2:1867409/1‑100 (MQ=255)
                        aaTAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGt                                                                     >  1:1961783/1‑100 (MQ=255)
                          tAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTAc                                                                   >  1:2690840/1‑100 (MQ=255)
                           accaacagcccGGACACCTGCACCAATGGACAGGGCCGGCCATCACGCCCACTCCGCGAAGGGGT‑ACAGGTGTTGATGCGCTTTAACGAGCCGATGTACa                                                                  <  1:2085876/91‑1 (MQ=255)
                           aacaacCTCCCGGACACCTCCACAAATTGGGTCGGCCAGCACTCACCCCAACACCGCGAAGGGGTGAAGGGTG‑TGATGCGCTTTAACGAGCCGATGAACa                                                                  <  1:318486/100‑1 (MQ=255)
                           aaCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACa                                                                  <  2:2697652/100‑1 (MQ=255)
                            aCAAACTGCCGGTCACCTGCACCACCTGTCACCGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACa                                                                  <  2:735452/99‑1 (MQ=255)
                            aCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACaa                                                                 >  1:1059554/1‑100 (MQ=255)
                            aCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACaa                                                                 >  2:152675/1‑100 (MQ=255)
                            aCAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACaa                                                                 >  1:2101209/1‑100 (MQ=255)
                             cAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAg                                                                <  2:1059568/100‑1 (MQ=255)
                             cAAACTGACGGACCCCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAg                                                                >  1:1820953/1‑100 (MQ=255)
                               aaCTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGt                                                              <  2:1431894/100‑1 (MQ=255)
                                 cTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTg                                                            >  2:1277985/1‑100 (MQ=255)
                                 cTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTg                                                            <  2:3231170/100‑1 (MQ=255)
                                   gCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTgg                                                           <  1:1009171/99‑1 (MQ=255)
                                   gCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGga                                                          >  1:1909413/1‑100 (MQ=255)
                                   gCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGga                                                          >  2:1108604/1‑100 (MQ=255)
                                    ccGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGgag                                                         <  1:382111/100‑1 (MQ=255)
                                     cGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGgaga                                                        >  1:593775/1‑100 (MQ=255)
                                      ggTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGagag                                                       <  1:2963913/100‑1 (MQ=255)
                                       gTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGc                                                      <  1:1152456/100‑1 (MQ=255)
                                         cACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAg                                                    >  2:3246549/1‑100 (MQ=255)
                                         cACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAg                                                    >  2:636643/1‑100 (MQ=255)
                                          aCCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGa                                                   >  1:830797/1‑100 (MQ=255)
                                           ccTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGaa                                                  >  2:2433302/1‑100 (MQ=255)
                                           ccTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGaa                                                  >  2:2188020/1‑100 (MQ=255)
                                             tGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACa                                                >  2:516315/1‑100 (MQ=255)
                                              gCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGCAGGGGTGGAGGATG‑TGATGCGCTTTACCGAGCCGATGTACAAGGTTGGAGAGCAGAACAg                                               >  1:1556343/1‑100 (MQ=255)
                                              gCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAg                                               >  1:2540451/1‑100 (MQ=255)
                                                aCCAACTGTCACGGCCAGCCATCACCGCAACACCGCGCAGGGGGGAAGGATG‑TGATGCGCTGTAACGAGCCGCCGTACAAGGTTGGAGAGCAGAACAGCg                                             >  1:1229032/1‑100 (MQ=255)
                                                aCCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCg                                             >  2:2688826/1‑100 (MQ=255)
                                                      tGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTa                                       <  1:746351/100‑1 (MQ=255)
                                                      tGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTa                                       <  2:2607972/100‑1 (MQ=255)
                                                       gTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGGGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTa                                       <  1:536696/99‑1 (MQ=255)
                                                        tCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATg                                     <  2:1872084/100‑1 (MQ=255)
                                                        tCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATg                                     >  2:2372330/1‑100 (MQ=255)
                                                         cACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGt                                    >  1:1464640/1‑100 (MQ=255)
                                                          aCGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTc                                   <  1:2245275/100‑1 (MQ=255)
                                                           cGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGGGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTcc                                  <  1:3050454/100‑1 (MQ=255)
                                                           cGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTcc                                  >  2:1670092/1‑100 (MQ=255)
                                                           cGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTcc                                  >  1:3218480/1‑100 (MQ=255)
                                                            gtccagccaTCACAGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCt                                 >  2:2356398/3‑100 (MQ=255)
                                                             gccagccaTCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTg                                >  1:1456247/1‑100 (MQ=255)
                                                              ccagccaTCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGt                               <  2:1215459/100‑1 (MQ=255)
                                                                agccaTCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAGCAGCGTCTGTATGTCCTGTCa                             <  2:2944468/100‑1 (MQ=255)
                                                                agccaTCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCa                             >  2:3010985/1‑100 (MQ=255)
                                                                agccaTCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCa                             <  2:2153195/100‑1 (MQ=255)
                                                                 gccaTCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCAt                            >  2:2398198/1‑100 (MQ=255)
                                                                   caTCACCGCACCACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGTGCAGAACAGCGTCTGTATGTCCTGTCATCt                          <  1:2279606/100‑1 (MQ=255)
                                                                   caTCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTACGTCCTGTCATCt                          <  2:2282378/100‑1 (MQ=255)
                                                                    aTCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTg                         >  2:2540890/1‑100 (MQ=255)
                                                                    aTCACCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTg                         >  2:1469279/1‑100 (MQ=255)
                                                                       aCCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCg                      <  2:2504277/100‑1 (MQ=255)
                                                                       aCCGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCg                      <  1:1566374/100‑1 (MQ=255)
                                                                        ccGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCgg                     <  1:3409088/100‑1 (MQ=255)
                                                                         cGCAACACCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGa                    >  1:2584910/1‑100 (MQ=255)
                                                                           cAACACCGCGAAGGGGTGAAAGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAAc                  <  1:1807216/100‑1 (MQ=255)
                                                                             acacCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAg                <  2:2653354/100‑1 (MQ=255)
                                                                              cacCGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGt               >  2:1717813/1‑100 (MQ=255)
                                                                                ctgcgAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTg             <  2:9383/98‑1 (MQ=255)
                                                                                ccGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTg             >  1:639443/1‑100 (MQ=255)
                                                                                ccGCGAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTg             <  2:1056745/100‑1 (MQ=255)
                                                                                  gcgAAGGGGTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTgag           >  1:9383/1‑98 (MQ=255)
                                                                                       ggggTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAg      <  2:1901065/100‑1 (MQ=255)
                                                                                       ggggTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAg      <  2:2259036/100‑1 (MQ=255)
                                                                                        gggTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGc     >  1:2396619/1‑100 (MQ=255)
                                                                                        gggTGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGc     <  2:125407/100‑1 (MQ=255)
                                                                                           tGAAGGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGtt  >  2:808585/1‑100 (MQ=255)
                                                                                               |                                                                                                
AAACATGCTTCCGTTATCAACCCGAATAACAAACTGCCGGTCACCTGCACCAACTGTCACGGCCAGCCATCACCGCAACACCGCGAAGGGGTGAAAGATG‑TGATGCGCTTTAACGAGCCGATGTACAAGGTTGGAGAGCAGAACAGCGTCTGTATGTCCTGTCATCTGCCGGAACAGTTGCAAAAAGCGTT  >  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: ‑