Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,083,214 T→C 6.1% N870S (AAC→AGC)  fdoG ← formate dehydrogenase‑O, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,083,2140TC 6.1% 229.1 / 4.3 81N870S (AAC→AGC) fdoGformate dehydrogenase‑O, large subunit
Reads supporting (aligned to +/- strand):  ref base T (36/40);  new base C (2/3);  total (38/43)
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.77e-01

CTGCTCTGGCTGCAAAATCGCGTTCAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCT  >  NC_000913/4083066‑4083362
                                                                                                                                                    |                                                                                                                                                    
cTGCTCTGGCTGCAAAATCGCGTTCAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTATAGATACGTGCAGCCGGGtt                                                                                                                                                     >  2:502955/1‑150 (MQ=255)
   ctctGGCTGCAAAATCGCGTTCAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCga                                                                                                                                                  <  2:367957/150‑1 (MQ=255)
      tgtcttcAAAATCTCGTTCAACTGCGCGTGTTTGGTCAAGTCTTGGAAGTGCTCGGTCTTACGATAGGGGGTTCCGTCATACGGGAACTTATCGTCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGAt                                                                                                                                               <  2:677525/144‑1 (MQ=255)
         cTGCAAAATCGCGTTCAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGCTCGAGATAAc                                                                                                                                            >  2:595367/1‑150 (MQ=255)
         cTGCAAAATCGCGTTCAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGCTCGAGATAAc                                                                                                                                            >  2:590122/1‑150 (MQ=255)
         cTGCAAAATCGCGTTCAACAGCGCGTGTTTGGTCCAGAAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAAc                                                                                                                                            <  1:764124/150‑1 (MQ=255)
           gCAAAATCGCGTTCAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGt                                                                                                                                          >  2:27093/1‑150 (MQ=255)
            cAAAATCGCGTTCAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGtt                                                                                                                                         <  2:942041/150‑1 (MQ=255)
             aaatTCGCGTTCAACAGGGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTGAAAGAAACGTGCAGCCGGGTTCGAGATAACGttt                                                                                                                                        <  2:64435/150‑1 (MQ=255)
              aaaTCGCGTTCAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTg                                                                                                                                       >  2:1031956/1‑150 (MQ=255)
                   gcgTTCAACAGCGCGTGTTTGGTCCAGTTGTGGAAGTGCTCGGTCAGTCGATAGGTGGTTCCGACATACGGGTACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGc                                                                                                                                  <  1:527514/150‑1 (MQ=255)
                   ccgTTCTACAGCGCGTGGTTGGTCCAGAAGGGGAAGTGCTCGGGCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGGCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGc                                                                                                                                  <  1:476720/149‑1 (MQ=255)
                        cAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCgg                                                                                                                             <  1:1094697/150‑1 (MQ=255)
                                       ggTCCTGTAGTGGAAGGGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAgcgg                                                                                                              <  1:155447/150‑1 (MQ=255)
                                          ccAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGt                                                                                                           <  1:994697/150‑1 (MQ=255)
                                           cAGTAGTGGAAGTGCTCGGGCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGtt                                                                                                          <  2:1234675/150‑1 (MQ=255)
                                            agtagtGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGttt                                                                                                         <  2:19815/150‑1 (MQ=255)
                                              tagtGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCa                                                                                                       <  2:408456/150‑1 (MQ=255)
                                                gtGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCaaa                                                                                                     >  2:1163067/1‑150 (MQ=255)
                                                gcggAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCaaa                                                                                                     <  2:1215800/148‑1 (MQ=255)
                                                 tGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAAc                                                                                                    >  1:622870/1‑150 (MQ=255)
                                                  ggAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACg                                                                                                   <  1:20462/150‑1 (MQ=255)
                                                  ggAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGCCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACg                                                                                                   <  1:672319/150‑1 (MQ=255)
                                                  ggAAGTGCGCGGGCCGACGCGAGGTGGTGCCGACATACGGGCACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACg                                                                                                   <  1:950265/150‑1 (MQ=255)
                                                       tGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCg                                                                                              >  1:146248/1‑150 (MQ=255)
                                                       tGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCg                                                                                              <  1:1308546/150‑1 (MQ=255)
                                                       tGCTCGGTCAGACGATAGGTGGGTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCg                                                                                              <  1:1088892/150‑1 (MQ=255)
                                                          tCGGTCAGACGATAGGTTGATCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAg                                                                                           <  2:1102236/150‑1 (MQ=255)
                                                            gggcAGACGATAGGCGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTGCGGCGTCGTCGGTAAAGATCCGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTGCCAGCGGCGTTTCAAACGGCTCGTAgtg                                                                                         <  1:655891/147‑1 (MQ=255)
                                                             gTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAgtgt                                                                                        >  2:6177/1‑150 (MQ=255)
                                                               cAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTGGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCATCGGCGTTTCAAACGGCTCGTAGTGATc                                                                                      >  2:55449/1‑150 (MQ=255)
                                                                 gACGATAGGTGGTTCCGACATCCGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGGTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCg                                                                                    <  1:351012/150‑1 (MQ=255)
                                                                  aCGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCgg                                                                                   <  1:929968/150‑1 (MQ=255)
                                                                   cGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCggg                                                                                  >  1:1009380/1‑150 (MQ=255)
                                                                   cGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACACAATGCTTGGGCGTCGTCTTTAAAGATACGTGCAGCCGGTTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCggg                                                                                  <  2:241653/150‑1 (MQ=255)
                                                                        ggtggtTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCTGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACgg                                                                             >  2:631059/1‑150 (MQ=255)
                                                                          tggtTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGAc                                                                           >  1:490800/1‑150 (MQ=255)
                                                                                cGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCg                                                                     >  1:616854/1‑150 (MQ=255)
                                                                                cGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCg                                                                     >  2:812982/1‑150 (MQ=255)
                                                                                 gACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGc                                                                    >  1:324326/1‑150 (MQ=255)
                                                                                 gACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGc                                                                    >  1:420977/1‑150 (MQ=255)
                                                                                 gACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGc                                                                    >  1:426907/1‑150 (MQ=255)
                                                                                      actggAACTTATCGGCTTTAGCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCt                                                               <  2:399875/147‑1 (MQ=255)
                                                                                        gggAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTa                                                             >  1:549836/1‑150 (MQ=255)
                                                                                           aaCTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCg                                                          >  1:406974/1‑150 (MQ=255)
                                                                                           aaCTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGCTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCg                                                          <  1:590122/150‑1 (MQ=255)
                                                                                           aaCTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGCTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCg                                                          <  1:595367/150‑1 (MQ=255)
                                                                                           aaCTTATCGGCATTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCg                                                          <  1:636943/150‑1 (MQ=255)
                                                                                             cTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAg                                                        >  2:860062/1‑150 (MQ=255)
                                                                                              ttATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAgg                                                       >  1:21523/1‑150 (MQ=255)
                                                                                                aTCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGc                                                     <  2:1122458/150‑1 (MQ=255)
                                                                                                aTCGGCTGTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGTACCTTCCGCCATCTTATCGAGGGc                                                     <  2:1016942/150‑1 (MQ=255)
                                                                                                   ggCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCaaa                                                  >  1:462431/1‑150 (MQ=255)
                                                                                                     cTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATTACGATAGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCTAACGCCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCATACa                                                >  2:1070302/1‑150 (MQ=37)
                                                                                                     cTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGCGGGCAAACa                                                >  2:112347/1‑150 (MQ=255)
                                                                                                       ttACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGa                                              <  1:776597/150‑1 (MQ=255)
                                                                                                         aCCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACg                                            <  1:502955/150‑1 (MQ=255)
                                                                                                                gCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATg                                     >  1:562882/1‑150 (MQ=255)
                                                                                                                 cTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGc                                    >  1:707231/1‑150 (MQ=255)
                                                                                                                  ttCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAACCAGACGCCCCATGcc                                   >  2:385437/1‑150 (MQ=255)
                                                                                                                  ttCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGcc                                   >  1:701010/1‑150 (MQ=255)
                                                                                                                     ggcggcGTCTTTAAATATACGTGCAGCCGGGTTCGAGATAACGTTTGGTAGCAGCGGGGTAGTTCCCAGCGGCGTTTCTAAAGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGCCGCCCCATGCCTTc                                <  1:53598/150‑1 (MQ=25)
                                                                                                                      gCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTcc                               <  2:141233/150‑1 (MQ=255)
                                                                                                                      gCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGACGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTcc                               <  2:802035/150‑1 (MQ=255)
                                                                                                                       cgtcgtCTTTAAAGATACGTGCAGCCGGGCTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCct                              <  1:1118995/150‑1 (MQ=255)
                                                                                                                        gtcgtcTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCctg                             <  1:583898/150‑1 (MQ=255)
                                                                                                                        ggcggcTTTAAAGTTACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCctg                             <  1:578791/145‑1 (MQ=255)
                                                                                                                          cgtcTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCctgct                           >  2:802670/1‑150 (MQ=255)
                                                                                                                            tcTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCtgctgc                         <  2:622870/150‑1 (MQ=255)
                                                                                                                             cTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCTTGCCGTCCTGCTGCa                        >  2:1094477/1‑150 (MQ=255)
                                                                                                                             cTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTAGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCa                        >  2:159136/1‑150 (MQ=255)
                                                                                                                               ttAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGCGGGCAAACAGCCGCCCCATGCCTTCCGGCTGCATg                      >  1:57056/1‑150 (MQ=255)
                                                                                                                               ttAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATg                      >  1:1178139/1‑150 (MQ=255)
                                                                                                                                 aaaGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGCGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGAt                    >  2:863158/1‑150 (MQ=255)
                                                                                                                                 aaaGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGAt                    >  2:453661/1‑150 (MQ=255)
                                                                                                                                  aaGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATa                   >  1:1038898/1‑150 (MQ=255)
                                                                                                                                         gTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGcc            <  2:283491/150‑1 (MQ=255)
                                                                                                                                           gCATCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCg          <  1:1031956/150‑1 (MQ=255)
                                                                                                                                             aGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGAc        >  1:653410/1‑150 (MQ=255)
                                                                                                                                                cGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGAcgcc     <  1:253424/150‑1 (MQ=255)
                                                                                                                                                  ggTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGAcgccgc   <  1:1215947/150‑1 (MQ=255)
                                                                                                                                                   gTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGACTTCCTGCTGCATGATAAACGGCCCGACGCCGCt  >  1:219895/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGCTCTGGCTGCAAAATCGCGTTCAACAGCGCGTGTTTGGTCCAGTAGTGGAAGTGCTCGGTCAGACGATAGGTGGTTCCGACATACGGGAACTTATCGGCTTTACCCAATGCTTCGGCGTCGTCTTTAAAGATACGTGCAGCCGGGTTCGAGATAACGTTTGGATGCAGCGGGTTAGTTCCCAGCGGCGTTTCAAACGGCTCGTAGTGTTCCGGGAACGGACCTTCCGCCATCTTATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCT  >  NC_000913/4083066‑4083362

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