Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,263 G→T 9.3% intergenic (+272/+135) 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,2630GT 9.3% 125.8 / 3.5 65intergenic (+272/+135)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (34/24);  new base T (5/2);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCA  >  NC_000913/1552116‑1552411
                                                                                                                                                   |                                                                                                                                                    
aTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                                                                                                     >  1:144429/1‑149 (MQ=255)
    aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTATAATGAGTCATTTTATTATTGTATACAACGTAGGCTTTGTTCATTCCGGATCCGGCGTTAAccc                                                                                                                                                 >  2:275192/1‑147 (MQ=17)
     aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTGATATTGTAGACAAAGTAGGCTTTGTTCATGCCGGATGCGGCGTGaaagac                                                                                                                                                >  1:13021/1‑145 (MQ=21)
     aTAGATTGCACTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                <  2:445967/149‑1 (MQ=255)
      tAGATTTCAGTGAACGTGTAGGCCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                               <  2:1029684/149‑1 (MQ=255)
              aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCGTCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCTTATTGTAGAAAACGTCGGCTTTGTTCATACAGAATCCGGCTTTAACGCCATATGCGGc                                                                                                                                       >  1:145667/1‑149 (MQ=18)
               gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                      <  1:703019/149‑1 (MQ=255)
                tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                                                                                                     <  2:527198/149‑1 (MQ=255)
                   aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAATGCGTCATTTTAATATTGTAGACAACGAAGGCTTTTTTCATGCAGGTAGCGACTTGAACGCATAATCCGGCATGaa                                                                                                                                  >  2:338905/1‑149 (MQ=21)
                          gtcatGATAAGCGTAGCGAATCAGGCAATGTTTCGTTTTGGATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAGGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCtt                                                                                                                           <  2:96151/145‑1 (MQ=17)
                           gCCTGATAAGCGTAGCGCTTCAGGCAATGTTGCGATTGTCATCAGTTTCAAAAGGCGCTGTAAAAGCCGTCATTTTCATATTGTAGACAACGTAGGCTTTCCGCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                                                                                          <  1:391596/149‑1 (MQ=21)
                           gCCGGATAAGCGTAGCTCATCATTCAATGTTGCTTTTGTAATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                                                                                          <  2:241457/149‑1 (MQ=17)
                            ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                                                                                                         >  1:956361/1‑149 (MQ=255)
                              ttaTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaaa                                                                                                                       <  2:186552/147‑1 (MQ=255)
                                   aGCGTAGCGCATCAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                                                                                                  <  2:381485/149‑1 (MQ=255)
                                   aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCaa                                                                                                                  <  2:144429/149‑1 (MQ=255)
                                      ggaGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAg                                                                                                               <  1:477636/147‑1 (MQ=255)
                                      gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATATTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGAATCAAAACCCTTCAAAACCAATAg                                                                                                               >  1:651333/1‑149 (MQ=18)
                                      gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTAATGCCGGATGCGTCTTGAACGCCTTATACGGCATGAAAAAACTTCAAATCAAATat                                                                                                               >  1:442317/1‑148 (MQ=18)
                                      gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTACAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGAAGTCTTTGTTCATTCCGAATGCGGCGTGAACGCCATATACGGCATGAAAAACCTTCaaataaaatat                                                                                                               >  1:461920/1‑148 (MQ=37)
                                      gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAAGAGTCATTTTAATATTGTAGAAAAAGTAGGCATTTATCATGACGGAGGCGTAGTGAAAGCCTTAACCGGCATGAAAAACCGTCTATTCCAAttt                                                                                                               >  2:56797/1‑147 (MQ=12)
                                      gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGAAGAAAAAGTAGGCATTGATCATGCCCGATGCGTAGTGAAAGACTTTTCCGGCATGAAAAACCTTCAACTCCAATaa                                                                                                               >  2:141174/1‑148 (MQ=21)
                                       tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTATTATTGTATACAAAGTAGGCTTTTTTCATGCCGGATGCGCCGTGAACGCCTTAAACTGCATTAACAGACTTAAAATCAATTAGa                                                                                                              >  1:666904/1‑149 (MQ=21)
                                                  ttaaTGTTGATTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAc                                                                                                   <  2:604061/147‑1 (MQ=255)
                                                    aaTGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgt                                                                                                 >  1:972109/1‑149 (MQ=255)
                                                       gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTa                                                                                              <  2:956361/149‑1 (MQ=255)
                                                            gTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGACCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGACt                                                                                         >  2:1182359/1‑149 (MQ=255)
                                                                   atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                  >  2:300445/1‑149 (MQ=39)
                                                                       gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAg                                                                              >  1:1126885/1‑149 (MQ=35)
                                                                          tCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAACCGGCATGAAAACCCTTCAAATCCAAAAGAATGCAGTGAACGAGAAGGCCTGATAAGCGTAgcgc                                                                           >  2:168101/1‑149 (MQ=11)
                                                                              aTGGCGCTGTAAAAGGCGTCATTTTCATATTTTAGACAACGTAGGCTTTGTTCATTCCGGATGCGGCGTGAACGCCTTATCCTTCATGAAAAAACTTCAAATACAATAGATTGAAGTTAACGTGTATGCCTGATAAGCGTAGCGAATCa                                                                       >  2:776310/1‑149 (MQ=11)
                                                                                    cTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt                                                                 >  1:544050/1‑149 (MQ=32)
                                                                                      gTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGt                                                               >  2:266119/1‑149 (MQ=32)
                                                                                       tAAAAGGCGTCATTTTCATATTGTCGACAACGTAGGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                              <  2:957384/149‑1 (MQ=17)
                                                                                       tAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                              >  1:629143/1‑149 (MQ=32)
                                                                                        aaaaGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGTGTAGCGCATCTGGCAATGTTg                                                             >  1:476487/1‑149 (MQ=17)
                                                                                        aaaaGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAAGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTg                                                             >  2:729219/1‑149 (MQ=17)
                                                                                             gCGTCATTTTCATATTGTAGAGAAAGTAGGCTTTGTTCATGCCGGATTCTGCGTGAACGCCTTATCCGGCTTGAAAACGCTTCAAATCCATTAGATTGCAGTGAACTTTTAGGCCTTATAAGCGTAGCGCATCAGGCAATGTTGCGttt                                                        <  1:239415/149‑1 (MQ=11)
                                                                                             gCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCATTCAAATACAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGttt                                                        >  1:748373/1‑149 (MQ=17)
                                                                                              cGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACACTTCAAATCCAATAGATTGCAGTTAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTg                                                       <  1:168101/149‑1 (MQ=17)
                                                                                                            tGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTaaa                                         <  2:513233/149‑1 (MQ=33)
                                                                                                             gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTTCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAAt                                        >  2:600281/1‑149 (MQ=21)
                                                                                                             gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCTGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAAt                                        >  1:741434/1‑149 (MQ=21)
                                                                                                             gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAAt                                        >  1:1048179/1‑149 (MQ=34)
                                                                                                             gTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACACTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGACATCAGTTCTAAAt                                        >  2:580418/1‑149 (MQ=255)
                                                                                                               aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATgg                                      <  1:64047/149‑1 (MQ=34)
                                                                                                                    aCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGAAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCtt                                 <  1:776310/149‑1 (MQ=25)
                                                                                                                            tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAACCCGTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGAGTAGCGCATCAGGAAATGTTGCGTGTGTCATAAGTTCTAAATGGCGCTTTATAAAGt                         >  2:1027185/1‑149 (MQ=14)
                                                                                                                            tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                         <  1:300445/149‑1 (MQ=255)
                                                                                                                             ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTg                        <  2:1048179/149‑1 (MQ=255)
                                                                                                                                ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACGCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATATAGTGCCa                     >  1:918280/1‑149 (MQ=255)
                                                                                                                                 tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAt                    >  2:594453/1‑149 (MQ=255)
                                                                                                                                 tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAAAGGCGCTTTAAAAAGAGACAt                    >  2:405783/1‑149 (MQ=11)
                                                                                                                                  cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtt                   <  1:406232/149‑1 (MQ=255)
                                                                                                                                    tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATACAATGCACTGAACGTGTAGGCATGATAAGCGAAGCGCAACAGGCAATGTATCGTTTGTCATCAGATCTAAATGGAGCTTTATAAAGTTCCAtttt                 >  2:684420/1‑149 (MQ=14)
                                                                                                                                    tGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                 >  2:64648/1‑149 (MQ=14)
                                                                                                                                     gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                >  1:449210/1‑149 (MQ=255)
                                                                                                                                       cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtttagaa              <  2:1001866/149‑5 (MQ=255)
                                                                                                                                       cGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAAGAGTCAttttaat              >  2:44309/1‑146 (MQ=1)
                                                                                                                                        ggataaggcgTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa             <  1:64648/143‑1 (MQ=17)
                                                                                                                                        ggataaggcgTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa             >  2:183721/7‑149 (MQ=17)
                                                                                                                                         gataaggcgTTCACGCCGCATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTAt            >  1:183151/6‑149 (MQ=17)
                                                                                                                                           tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTg          <  1:975436/149‑1 (MQ=255)
                                                                                                                                               gcgTTAACGCATTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTa      <  2:442317/149‑1 (MQ=255)
                                                                                                                                               gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTa      >  1:1010500/1‑149 (MQ=255)
                                                                                                                                                   gAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAATTTTTTATTGCGTAAcaa  >  1:616270/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCA  >  NC_000913/1552116‑1552411

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