Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,191 C→A 16.5% intergenic (+200/+207) 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,1910CA16.5% 177.5 / 9.3 89intergenic (+200/+207)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (20/51);  new base A (12/5);  total (33/56)
Fisher's exact test for biased strand distribution p-value = 1.80e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGC  >  NC_000913/1552042‑1552338
                                                                                                                                                     |                                                                                                                                                   
atatTGTAGACAACGTAGTCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                     <  1:1210118/150‑1 (MQ=14)
atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTc                                                                                                                                                     <  1:1789637/150‑1 (MQ=31)
  atTGTAGACAACGTAGGCTTTGTTCATGACGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTTCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaa                                                                                                                                                   <  1:408895/150‑1 (MQ=17)
   ttGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaaa                                                                                                                                                  <  1:1245072/150‑1 (MQ=33)
        gACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcg                                                                                                                                             <  1:1292104/150‑1 (MQ=34)
           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGACCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                          <  1:1794159/150‑1 (MQ=25)
           aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                                                                                                          <  2:567805/150‑1 (MQ=35)
             cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTa                                                                                                                                        <  2:517971/150‑1 (MQ=37)
                  gCTTTGTTCATGCCGGATGCTGCGTGAACGCCTTATCCGGCATGAAAACCGTTCATATCCAATAGATTGCAGTGAACGTGTAGTTCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCTTCAGTTTCAAATGGCGCTGTAAAAgg                                                                                                                                   <  1:847910/150‑1 (MQ=14)
                  gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                                                                                                                                   <  2:1466872/150‑1 (MQ=255)
                      tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTc                                                                                                                               <  2:1236219/150‑1 (MQ=255)
                      tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAggagac                                                                                                                               >  2:2008799/1‑146 (MQ=255)
                        ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAt                                                                                                                             <  2:655497/150‑1 (MQ=255)
                        ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAATGAGTCAt                                                                                                                             >  2:876096/1‑150 (MQ=11)
                          cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTGAAATGGCGCTGTAAAAGGCGTGAttt                                                                                                                           >  1:1161790/1‑150 (MQ=255)
                          cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAttt                                                                                                                           <  2:1602809/150‑1 (MQ=255)
                           aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGGCAtttt                                                                                                                          >  1:549193/1‑150 (MQ=255)
                            tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                         <  1:1298363/150‑1 (MQ=255)
                            tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                         <  2:2174843/150‑1 (MQ=255)
                            tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                         <  2:1057495/150‑1 (MQ=255)
                             gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa                                                                                                                        <  1:1480387/150‑1 (MQ=255)
                             gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa                                                                                                                        <  1:468030/150‑1 (MQ=255)
                             gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCa                                                                                                                        <  1:1187481/150‑1 (MQ=255)
                              ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                                                                                                       <  1:563663/150‑1 (MQ=255)
                               cGGATGCGGCGTGAACGCCTTATCCGGTATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCata                                                                                                                      <  1:812618/150‑1 (MQ=255)
                                   ttaggcgTGAACGCCTTATCCTTCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                  <  1:603575/147‑1 (MQ=255)
                                    gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTTTAGGCATGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa                                                                                                                 <  1:38150/150‑1 (MQ=255)
                                    gcggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa                                                                                                                 <  1:1443829/150‑1 (MQ=255)
                                    gcggcgTGAACGCCTTATCCGGCATGAAAACCATTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTa                                                                                                                 <  1:341464/150‑1 (MQ=255)
                                     cggcgTGAACGCGTTATCCTGCATGAAAACCCTTCAAATCCAATATATTGCAGTGAACGTGTAGGACTCATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                <  2:1968702/150‑1 (MQ=17)
                                     cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCTAATCAGGAGATCTTTAGATAGTCATCAGTTTAAAAAGTCGATTTAAAAGGAGTAATTTTAATATTGTat                                                                                                                >  2:1373759/1‑149 (MQ=14)
                                     cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                <  1:35136/150‑1 (MQ=255)
                                      ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGa                                                                                                               <  1:654704/150‑1 (MQ=255)
                                       gcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGAc                                                                                                              <  1:1751709/150‑1 (MQ=255)
                                        cgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACa                                                                                                             <  2:914776/150‑1 (MQ=255)
                                         gTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACaa                                                                                                            <  1:386715/150‑1 (MQ=255)
                                          tGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATATATTGCAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc                                                                                                           <  2:1283879/150‑1 (MQ=255)
                                           gAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTATAAGGAGACATTTTAATATTGAAGACaaag                                                                                                          >  2:1512691/1‑148 (MQ=11)
                                               gCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAgg                                                                                                      <  2:2157796/150‑1 (MQ=255)
                                                acTTATCCTGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGTCCTGATAAGAGTATCGCTTCATTAAATGTTGCTTTTTTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                                                                                                     <  2:1403/149‑1 (MQ=18)
                                                   tATCCTGAATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAAGTGTAGGCCTGATAAGCGTAGCGCATCAGTCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCttt                                                                                                  <  2:2121737/150‑1 (MQ=255)
                                                     tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGt                                                                                                <  1:655866/150‑1 (MQ=255)
                                                      ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                               >  1:1575914/1‑150 (MQ=255)
                                                      ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                               <  2:836428/150‑1 (MQ=255)
                                                         gCATGAAAACCCTTCAAATCAAATAGATTTCCGTGAACGTGTAGGCTTGTTTATGTTATCTATTCAGTTAATGTTGTTTTTTTCATAAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCAt                                                                                            <  2:1476175/150‑1 (MQ=12)
                                                             gAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAATGAGTCATTTTAATATTGTAGAAAACGTAGGCATTGTTCATGCCg                                                                                        >  2:1654772/1‑150 (MQ=18)
                                                                   ccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcg                                                                                  >  2:463954/1‑150 (MQ=255)
                                                                    cTTCAATTTCATTAGATTGCATTGTACGTGTAGTAATTATATGCGTAGCGAATCAGGCAATGTTGAGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcgg                                                                                 <  2:1959072/150‑1 (MQ=21)
                                                                      tCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                               <  2:727366/150‑1 (MQ=255)
                                                                       cAAATCCAATAGATAGCAGTGAACGTGTAGGCCTGATAAGCGTAGCTCATCAGGCAATGTTGAGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGt                                                                              <  2:802040/150‑1 (MQ=255)
                                                                        aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                             <  2:678651/150‑1 (MQ=255)
                                                                           tCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                          <  1:1859320/150‑1 (MQ=255)
                                                                            ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGACATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                         <  1:1390066/150‑1 (MQ=255)
                                                                             cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGACATTTTCATATTGAAGAAAACGAAGGCTTTGTTCATCCCCGATGCGGATTGAAAGc                                                                        >  2:5695/1‑150 (MQ=18)
                                                                             cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGAGTAGCGCATCAGTCAATGTTGCTATTGTCATAAGTTTCAAATGGCGCTGTAAAAGGAGACAATTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACgg                                                                        >  1:1119982/1‑149 (MQ=25)
                                                                              aaTAGATTGCAGTGAACTTGTATGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTTTTCATGCCGGATGCGGCGTGAACGcc                                                                       <  2:2087032/150‑1 (MQ=25)
                                                                              aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                       <  2:1037377/150‑1 (MQ=255)
                                                                                       ctgttaaCGTGTAGGCCTGATAAGCGTATCGCATCAGGCAATGTTGAGTTTGTCTTCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                              <  2:858539/145‑1 (MQ=21)
                                                                                       ctgttaaCGTGTAGGCATGATATTCGTAGCTCTTCAGGAAATGTTGCTTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                              <  2:850854/145‑1 (MQ=18)
                                                                                       cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCATTAAAAGGAGACATATTAATATAGTAGACAACGTAAACTTTGTTAATACCTGAAGCTCAGTGAACTCCatataagaa                                                              >  2:1379161/1‑145 (MQ=9)
                                                                                        aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                             <  2:1923693/150‑1 (MQ=255)
                                                                                          tGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGTCATATTCATATAGTAGACAACGTAGGCTTTGTTAATTCCGGATGCGGAGTGAACGCCTTATAAGGCAta                                                           >  2:470993/1‑149 (MQ=21)
                                                                                            aaCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaa                                                         <  1:1763913/150‑1 (MQ=255)
                                                                                              cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGATGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaaa                                                       >  2:2051798/1‑150 (MQ=255)
                                                                                              cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaaa                                                       >  2:8530/1‑150 (MQ=255)
                                                                                               gtgtAGGCCTGATAAGCGGAGGGCATGAGGCAAGGTTGCGTTTGTCATCAGATTCAAATGGCGCTGTAAAAGGCGTCATTTTCATAATGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAc                                                      <  2:1110067/150‑1 (MQ=18)
                                                                                                tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAAAGACTTATACGGAATGAAAAcc                                                     >  1:2182851/1‑150 (MQ=25)
                                                                                                tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                     >  2:753192/1‑150 (MQ=255)
                                                                                                tgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAcc                                                     >  1:575127/1‑150 (MQ=255)
                                                                                                   aGGCCTGATAAGCGTAGCGCATAAGGCAATGTTGAGTTTGTCATCAGTTTCAAATGGCGCTGAAAAAGTCGTCATTTTCATATTGTAGACAAAGTAGGCTTTGTTCATGCCGGATGCGGCGTGACCGCCGTATCCGGCATGCAAGcccgt                                                  >  2:119077/1‑148 (MQ=21)
                                                                                                      ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCaa                                               <  2:1575914/150‑1 (MQ=255)
                                                                                                           tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGAAAACGTAGGCTTTGTTCATTCCGGAAGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTCAAATCCa                                          >  2:712139/1‑150 (MQ=255)
                                                                                                              gCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                       >  1:1009524/1‑150 (MQ=255)
                                                                                                                    cgcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                 <  2:575127/150‑1 (MQ=255)
                                                                                                                    cgaaTCAGGCAATGTTGCGTTTGGCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                 <  2:1416964/147‑1 (MQ=255)
                                                                                                                     gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTATAAAACGTAGGCTTTGTACATGCCGGAAGCGGCGTGAACGACTTAACCGGAAAGAAAACCCTTCAAAACCAAAAGATTgaa                                >  1:874179/1‑148 (MQ=37)
                                                                                                                          aGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTAAAAGGAGACATTTTAATATTGTATACAACGAAGGCTTTGTACATGCCCGATGCGGCGTGAAACACAAAAACGCCATGACAACACATCAAATACAATAGAATGCATTGaa                           >  2:1172479/1‑150 (MQ=16)
                                                                                                                           ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATTCCCGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATATATTGCAGTGAAc                          >  1:1613395/1‑150 (MQ=255)
                                                                                                                           ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTAAAATACAATATATTGCAGTTAAc                          >  1:1157324/1‑150 (MQ=17)
                                                                                                                            ttaaTGTTGCGTTTTTCATCAGGTGAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                         <  2:1203718/148‑1 (MQ=255)
                                                                                                                            gCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCCGAAGCGGCGTGAACGACTGATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACg                         >  2:2069021/1‑150 (MQ=17)
                                                                                                                                 gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCATAACCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                    >  2:2122554/1‑150 (MQ=255)
                                                                                                                                  ttGCGTTTGTCATCAGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAgg                   <  2:660884/150‑1 (MQ=255)
                                                                                                                                   tGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTAAAATCCAATAGATTGCAGTGAACGTGTAGGc                  >  2:2026068/1‑150 (MQ=255)
                                                                                                                                    gCGTTTGTCATCAGTTCCAAATGGCGCTGAAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGATCATGCCGTACGCGGCTTTAACGCCTTATCCGTCATGAATACCCTTCTAATCCAATAGATTGCAGTGAACGTGTAGGcc                 <  2:1119982/150‑1 (MQ=14)
                                                                                                                                     cGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGCTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                >  2:665821/1‑150 (MQ=255)
                                                                                                                                          gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGACCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa           >  2:1414906/1‑150 (MQ=255)
                                                                                                                                            catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc         >  1:300113/1‑150 (MQ=40)
                                                                                                                                            catcaGTTTCAAATGGAGCTGTAAATGGCGGCATTTTCATATTGTAGTCAACGTAGGCTTTGTGCATGCCGGATGCGTCGTGAACGCCTTATCCGGCATGAAAACCCTTAAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGTTAAGc         <  1:119077/150‑1 (MQ=11)
                                                                                                                                                  tttCAAATGGCGCTGTAAAAGGCGTCATTTTCATTTTGTAGACAACGTAGGCTTTGTTCATGCCGTATGCGGCGTGAACGCCTTATCCGGCATGAAAACACTTCAAATCCAAATGATTGCAGGAAACGTGTAGGCCTCATAACAGTAgag   >  1:406462/1‑148 (MQ=11)
                                                                                                                                                   ttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATACAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcgc  >  2:1395580/1‑150 (MQ=21)
                                                                                                                                                     |                                                                                                                                                   
ATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGC  >  NC_000913/1552042‑1552338

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: