Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,551,710 T→C F125L (TTC→CTC)  fdnI → formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,551,7100TC96.9% 316.6 / ‑7.0 96F125L (TTC→CTC) fdnIformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base C (50/43);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CATCACAACATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGC  >  NC_000913/1551566‑1551855
                                                                                                                                                |                                                                                                                                                 
cATCACAACATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGt                                                                                                                                               >  2:824636/1‑149 (MQ=255)
cATCACAACATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAAAGCCGGGCAAAAGATGGTGTTCTGGGCGATCATGAGCATGGTTCTCGt                                                                                                                                               >  4:162577/1‑149 (MQ=255)
 aTCACAACATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGtgct                                                                                                                                            >  1:1087334/1‑151 (MQ=255)
    acaacaTCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTgctgct                                                                                                                                         <  1:1100074/151‑1 (MQ=255)
    acaacaTCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTgctgct                                                                                                                                         <  2:670756/151‑1 (MQ=255)
     caacaTCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTTGTCGATCATGAGCATGATTCTCGGgctgct                                                                                                                                         >  2:725585/1‑150 (MQ=255)
     caacaTCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCAGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGctgctg                                                                                                                                        <  2:128317/151‑1 (MQ=255)
       acaTCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGGACAAAGCCGGGGAAAAGATTATGTTCTGGGCGATCATGAGCATGATTCTCGTGCTGCTgg                                                                                                                                       >  2:766141/1‑150 (MQ=255)
        caTCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGGAAGTACAACGCCGGGCCAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGGGCTGCTGGTg                                                                                                                                     >  1:406714/1‑151 (MQ=255)
          tCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGa                                                                                                                                    <  2:1041256/150‑1 (MQ=255)
             cGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCgg                                                                                                                                <  3:228708/151‑1 (MQ=255)
             cGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCgg                                                                                                                                >  4:47051/1‑151 (MQ=255)
              ggATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCggg                                                                                                                               <  2:659398/151‑1 (MQ=255)
              ggATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCggg                                                                                                                               >  1:755725/1‑151 (MQ=255)
               gATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCgggg                                                                                                                              <  4:20536/151‑1 (MQ=255)
                aTAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCgggg                                                                                                                              >  4:400555/1‑150 (MQ=255)
                aTAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGt                                                                                                                             <  3:179255/151‑1 (MQ=255)
                    gaaagaTATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGAtt                                                                                                                         <  1:162830/151‑1 (MQ=255)
                         atatTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                                                                                                    >  1:166407/1‑151 (MQ=255)
                           atTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGc                                                                                                                  <  1:215967/151‑1 (MQ=255)
                            ttCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGGACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGGGACCGGGGGGATTATCTGGCg                                                                                                                 >  2:61876/1‑151 (MQ=255)
                             accGTGGATGTTGAACATTGTCGAAGTATTGAAAGGCCATGAGCATAAAGTGGCGGATGGCGGCGCGGACAACGACGGGCGAAAAATCATCTTCTTGTTGATCTTAAACACGATTCTCTTGCTCTTGGTGACCGGGGTGGTTTTCTgg                                                                                                                   >  2:739943/2‑148 (MQ=255)
                              ccGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGGTGTCGGTAAGGACAAAGCCGGGCAAAAGATGATGTTCTGGTCGATTATGAGCATGAATCTCGGGCTGCTGGGGACCGGGGGGATTATTTGGCg                                                                                                                 >  2:166587/1‑149 (MQ=255)
                                 tGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGGGGTCGGTATGGGTAACCCCGGGCAAAAAGTGAAGTTTTGGTCGATACTCAGCCTGGTTCTCGTGCGGTTGGTGACCGGGGTGGTTATTTGGCGTcc                                                                                                              >  2:1103015/1‑149 (MQ=255)
                                  ggCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGt                                                                                                            <  4:77571/150‑1 (MQ=255)
                                  ggCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTa                                                                                                           <  1:1247705/151‑1 (MQ=255)
                                  ggCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTATGGTCGAACCAGAGAATGATTTTCGAGCTGCTCGTGACCGGGGTGATTATCTGGCGTCCg                                                                                                             >  1:954286/1‑149 (MQ=255)
                                   gCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTAc                                                                                                          >  1:422618/1‑151 (MQ=255)
                                   gCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATAAATCTCGTGCTGCTGGGGACCGGGGAGATTATCTGGCGTCCg                                                                                                             >  1:189034/1‑148 (MQ=255)
                                      gTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACtt                                                                                                        <  4:16228/150‑1 (MQ=255)
                                      gTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTCGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGCCCGGGGTGATTCTCTAGCGTCCGTACt                                                                                                         <  2:498482/149‑1 (MQ=255)
                                      gTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGGAAGTAAAACGCCCGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGGGCTGCTGGTGACCGGGGTGGTTATCTGGCGTCCGTAc                                                                                                          >  2:46208/1‑148 (MQ=255)
                                       ttGAACATTGTCGAAGTATTGAAAGGAAATGAGCATAAAGTGGCGGATGCCGGTAAGTACAATGCCGGGGAAAAAATGAAGATCTGGGCGATCATTAGCATGATTTTCGAGCTGATGGTGACCGGGGTGATTATCTGGTGTCCGTACttt                                                                                                       >  4:185945/1‑150 (MQ=255)
                                          aaCATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTgcgc                                                                                                   >  1:471455/1‑151 (MQ=255)
                                          aaCATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTgcgc                                                                                                   <  1:876196/151‑1 (MQ=255)
                                           aCATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGGAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGGGATTATCTGGCGTCCGTACTTTgcgc                                                                                                   >  2:1108098/1‑150 (MQ=255)
                                                 tCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGGCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACtt                                                                                            >  1:188719/1‑151 (MQ=255)
                                                 tCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGACGGTAAGAAAAAAGCCGGGCAAAAGATGATGTTACGGTCGGTCATGAGCATGATTTTCGTGCTGCCGGTGACCGGGGTGATTATTTGGCGTCCGTAATTTGCGCAGTACtt                                                                                            >  2:123770/1‑151 (MQ=255)
                                                     aGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGGGATTATCTGGCGTCCGTACTTTGCGCAGTACTTccc                                                                                         >  2:826821/1‑150 (MQ=255)
                                                      gTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTccc                                                                                         >  2:105468/1‑149 (MQ=255)
                                                      gTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGGGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCg                                                                                        >  3:411013/1‑150 (MQ=255)
                                                        aTTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGAt                                                                                      <  4:255558/150‑1 (MQ=255)
                                                        aTTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATg                                                                                     >  4:399090/1‑151 (MQ=255)
                                                        aTTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGACGGGCAAAAGGTGATGGTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGGGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGAt                                                                                      >  2:730166/1‑150 (MQ=255)
                                                              aGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGgtt                                                                               <  2:1057814/151‑1 (MQ=255)
                                                                  aaTGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATTAGCATGATTCTCGTGCTGCTGGTGACCGGGGGGATTATCTGGCGTCCGGACTTTGCGCAGTAATTCCCGATGGCGgttgtt                                                                            >  3:386237/1‑150 (MQ=255)
                                                                  aaTGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTc                                                                           <  2:226064/151‑1 (MQ=255)
                                                                   aTGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCg                                                                          <  2:331269/151‑1 (MQ=255)
                                                                   aTGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGGCGATCATGAGCATGATTCTCGGGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGgttgtt                                                                            >  3:176161/1‑149 (MQ=255)
                                                                    tGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGGGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCg                                                                          >  1:602293/1‑150 (MQ=255)
                                                                          tAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACa                                                                     >  2:715766/1‑149 (MQ=255)
                                                                            aaGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCct                                                                 <  1:655060/151‑1 (MQ=255)
                                                                                ggCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGTATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCctgct                                                              >  1:58077/1‑150 (MQ=255)
                                                                                    gATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGGTcc                                                         <  2:1121469/151‑1 (MQ=255)
                                                                                       gTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCAcg                                                      <  1:653364/151‑1 (MQ=255)
                                                                                         cGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCAcgc                                                     <  3:427881/150‑1 (MQ=255)
                                                                                         cGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCAcgc                                                     <  1:1181404/150‑1 (MQ=255)
                                                                                          ggTAAGGACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGGGACCGGGGGGATTATCTGGCGTCCGTACTTTGCGCAGGACTTTCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCgg                                                   >  3:432978/1‑151 (MQ=255)
                                                                                           gTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGc                                                  <  3:2615/151‑1 (MQ=255)
                                                                                           gTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGc                                                  >  4:265443/1‑151 (MQ=255)
                                                                                            tAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCt                                                 >  1:149971/1‑151 (MQ=255)
                                                                                             aaGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTg                                                >  2:1115260/1‑151 (MQ=255)
                                                                                               gTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCg                                              <  4:8718/151‑1 (MQ=255)
                                                                                                    aCGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTat                                         <  1:833380/151‑1 (MQ=255)
                                                                                                    aCGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTat                                         >  4:448343/1‑151 (MQ=255)
                                                                                                     cGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTatc                                        >  1:97070/1‑151 (MQ=255)
                                                                                                      gCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTatc                                        >  2:776667/1‑150 (MQ=255)
                                                                                                      gCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTatc                                        >  1:614554/1‑150 (MQ=255)
                                                                                                        cGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTatcatc                                     >  1:1135832/1‑151 (MQ=255)
                                                                                                          ggCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCt                                   <  2:774832/151‑1 (MQ=255)
                                                                                                             aaaaGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGAt                                <  4:27676/151‑1 (MQ=255)
                                                                                                             aaaaGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGAt                                >  4:334181/1‑151 (MQ=255)
                                                                                                              aaaGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATc                               <  2:1036270/151‑1 (MQ=255)
                                                                                                               aaGATGATGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATcc                              >  4:359800/1‑151 (MQ=255)
                                                                                                                 gatgatGTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGGGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCAc                            >  1:1019205/1‑151 (MQ=255)
                                                                                                                    gatgTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGcc                         >  2:192900/1‑151 (MQ=255)
                                                                                                                     atgTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGcc                         >  4:252482/1‑150 (MQ=255)
                                                                                                                     atgTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGcc                         >  4:249281/1‑150 (MQ=255)
                                                                                                                       gTTCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATc                      >  1:610788/1‑151 (MQ=255)
                                                                                                                         tCTGGTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCt                    >  4:22230/1‑151 (MQ=255)
                                                                                                                            ggTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTg                   <  2:1110918/149‑1 (MQ=255)
                                                                                                                            ggTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGa                  <  3:173155/150‑1 (MQ=255)
                                                                                                                            ggTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGa                  <  1:381408/150‑1 (MQ=255)
                                                                                                                            ggTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGa                  <  4:325985/150‑1 (MQ=255)
                                                                                                                            ggTCGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGa                  <  2:591767/150‑1 (MQ=255)
                                                                                                                               cGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATcc               >  3:348881/1‑150 (MQ=255)
                                                                                                                               cGATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCa              <  1:544729/151‑1 (MQ=255)
                                                                                                                                gATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCat             <  1:753297/151‑1 (MQ=255)
                                                                                                                                gATCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCat             >  1:1109363/1‑151 (MQ=255)
                                                                                                                                 aTCATGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCata            <  4:137161/151‑1 (MQ=255)
                                                                                                                                    aTGAGCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGt         <  4:34787/151‑1 (MQ=255)
                                                                                                                                        gCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTatat     >  4:309814/1‑151 (MQ=255)
                                                                                                                                        gCATGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTatat     <  3:22230/151‑1 (MQ=255)
                                                                                                                                         cATGATTCTCGTGCTGCTGGTGACCGGGGTGGTTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATg    <  3:131762/151‑1 (MQ=255)
                                                                                                                                          aTGATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATgg   <  2:226002/151‑1 (MQ=255)
                                                                                                                                            gATTCTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGc  <  2:166383/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CATCACAACATCCCGGATAAGAAAGATATTCCGTGGCTGTTGAACATTGTCGAAGTATTGAAAGGCAATGAGCATAAAGTGGCGGATGTCGGTAAGTACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGCTACAGCCTGCTGATCCACGCGGCTGCGGGTATCATCCTGATCCACGCCATCCTGATCCATATGTATATGGC  >  NC_000913/1551566‑1551855

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