Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,942,635 (T)7→6 intergenic (‑68/+283) gcvA ← / ← ygdI glycine cleavage system transcriptional activator; autorepressor/DUF903 family verified lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,942,6290T.98.6% 329.3 / ‑1.3 70intergenic (‑62/+289)gcvA/ygdIglycine cleavage system transcriptional activator; autorepressor/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (38/31);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CTCTTCTGCTGCGCGAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTG  >  NC_000913/2942481‑2942775
                                                                                                                                                    |                                                                                                                                                  
ctctTCTGCTGCGCGAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGttt                                                                                                                                                  >  2:273565/1‑151 (MQ=255)
ctctTCTGCTGCGCGAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGttt                                                                                                                                                  >  2:30559/1‑151 (MQ=255)
 tcttctGCTGCGCGAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGtttt                                                                                                                                                 >  1:631152/1‑151 (MQ=255)
 tcttctGCTGCGCGAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGtttt                                                                                                                                                 >  2:473340/1‑151 (MQ=255)
    tctGCTGCGCGAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGtttttt                                                                                                                                               <  2:1216167/150‑1 (MQ=255)
         tgcGCGAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTg                                                                                                                                        >  2:1154320/1‑151 (MQ=255)
          gcgcGAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGa                                                                                                                                       <  2:818100/151‑1 (MQ=255)
           cgcgAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTg                                                                                                                                        >  1:226109/1‑149 (MQ=255)
                gTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCAtt                                                                                                                                 <  2:764208/151‑1 (MQ=255)
                   aaaCTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATtata                                                                                                                              <  1:32569/151‑1 (MQ=255)
                       ttAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAAtt                                                                                                                          >  1:587805/1‑151 (MQ=255)
                        tAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAAt                                                                                                                           <  1:1258437/149‑1 (MQ=255)
                        tAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAAt                                                                                                                           <  1:584922/149‑1 (MQ=255)
                        tAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTg                                                                                                                         >  2:1053305/1‑151 (MQ=255)
                               cGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTg                                                                                                                  <  2:1185978/151‑1 (MQ=255)
                               cGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTg                                                                                                                  >  1:46137/1‑151 (MQ=255)
                                 tgctgcGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTg                                                                                                                  >  2:10297/1‑149 (MQ=255)
                                  gctgcGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAg                                                                                                                <  2:878672/150‑1 (MQ=255)
                                  gctgcGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGc                                                                                                               <  2:215846/151‑1 (MQ=255)
                                       ggCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCt                                                                                                           >  1:673616/1‑150 (MQ=255)
                                            cAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAg                                                                                                      <  2:807547/150‑1 (MQ=255)
                                               aaaCTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCaaa                                                                                                  <  2:1337658/151‑1 (MQ=255)
                                                  cTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATAc                                                                                               >  1:846831/1‑151 (MQ=255)
                                                   tCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATAcc                                                                                              >  2:41318/1‑151 (MQ=255)
                                                   tCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATAcc                                                                                              >  1:1271142/1‑151 (MQ=255)
                                                      tAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCtat                                                                                           >  2:323424/1‑151 (MQ=255)
                                                       aaGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCtata                                                                                          <  2:1310222/151‑1 (MQ=255)
                                                         ggCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGt                                                                                        <  2:295035/151‑1 (MQ=255)
                                                          gCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTg                                                                                       <  2:949711/151‑1 (MQ=255)
                                                            aTTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTggc                                                                                     >  1:416913/1‑151 (MQ=255)
                                                              ttAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTggcgg                                                                                   >  2:372843/1‑151 (MQ=255)
                                                                  cGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACt                                                                               <  1:830424/151‑1 (MQ=255)
                                                                         aaTCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTg                                                                          <  1:14607/150‑1 (MQ=255)
                                                                         aaTCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTg                                                                          >  1:697620/1‑149 (MQ=255)
                                                                         aaTCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAg                                                                        <  2:617093/151‑1 (MQ=255)
                                                                            cGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCg                                                                     >  2:1127626/1‑151 (MQ=255)
                                                                             gTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCgg                                                                    >  2:349565/1‑151 (MQ=255)
                                                                             gTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCgg                                                                    <  2:1244551/151‑1 (MQ=255)
                                                                               ttAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGa                                                                   <  1:880664/150‑1 (MQ=255)
                                                                               ttAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGa                                                                   >  2:1053027/1‑150 (MQ=255)
                                                                               ttAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGa                                                                   <  2:1285153/150‑1 (MQ=255)
                                                                                tAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGaa                                                                  <  1:1294006/150‑1 (MQ=255)
                                                                                 aGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAAc                                                                 <  2:1027400/150‑1 (MQ=255)
                                                                                 aGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACg                                                                >  2:872726/1‑151 (MQ=255)
                                                                                   aCATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGa                                                               >  2:1329532/1‑150 (MQ=255)
                                                                                     aTGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGaaa                                                             >  1:805202/1‑150 (MQ=255)
                                                                                       ggCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGt                                                          >  2:509850/1‑151 (MQ=255)
                                                                                              aaaCTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCg                                                   <  1:473357/151‑1 (MQ=255)
                                                                                               aaCTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCgg                                                  >  2:758718/1‑151 (MQ=255)
                                                                                                aCTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGa                                                 <  2:237638/151‑1 (MQ=255)
                                                                                                   ttGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATg                                              >  2:299901/1‑151 (MQ=255)
                                                                                                    tGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATg                                              <  1:739542/150‑1 (MQ=255)
                                                                                                               gAATTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGt                                   <  2:587781/150‑1 (MQ=255)
                                                                                                                aaTTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTg                                  >  2:803915/1‑150 (MQ=255)
                                                                                                                aaTTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTg                                  >  2:93037/1‑150 (MQ=255)
                                                                                                                 aTTTAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGa                                 >  2:1223222/1‑150 (MQ=255)
                                                                                                                    tAACAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACtt                             >  1:1272807/1‑151 (MQ=255)
                                                                                                                     aaCAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACttt                            <  1:663374/151‑1 (MQ=255)
                                                                                                                      aCAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACtttt                           <  1:1127119/151‑1 (MQ=255)
                                                                                                                       cAATTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTg                          <  2:1313108/151‑1 (MQ=255)
                                                                                                                        aaTTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTgg                         >  2:1327462/1‑151 (MQ=255)
                                                                                                                        aaTTAGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTgg                         >  1:750997/1‑151 (MQ=255)
                                                                                                                            aGATCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCtt                      >  1:1136659/1‑150 (MQ=255)
                                                                                                                              aTCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACg                   >  1:646581/1‑151 (MQ=255)
                                                                                                                              aTCACACTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACg                   >  1:879642/1‑151 (MQ=255)
                                                                                                                                cacaCTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGt                 >  1:285424/1‑151 (MQ=255)
                                                                                                                                  cacTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTtg               >  1:1314805/1‑151 (MQ=255)
                                                                                                                                   acTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTtgt              <  2:414184/151‑1 (MQ=255)
                                                                                                                                   acTATGTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTtgt              >  1:394107/1‑151 (MQ=255)
                                                                                                                                        gTAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgt         <  1:758136/151‑1 (MQ=255)
                                                                                                                                         tAACCTATTAG‑TTTTTTAATCTGAGCCATTATAGATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgtt        >  2:281135/1‑151 (MQ=255)
                                                                                                                                         tAACCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgtt        >  1:245770/1‑151 (MQ=255)
                                                                                                                                          aaCCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttg       >  2:545610/1‑151 (MQ=255)
                                                                                                                                           aCCTATTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgt      <  1:1281599/151‑1 (MQ=255)
                                                                                                                                               aTTAG‑TTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGAtgttgtgttg  <  2:918726/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTCTTCTGCTGCGCGAGTGAAACTTAAATGGCGTGCTGCGGCATCAAAAACTCGTAAGGCATTTAGCGGTGGTAATCGTTTAGACATGGCTATTAAACTTTGATGTTAAATGAATTTAACAATTAGATCACACTATGTAACCTATTAGTTTTTTTAATCTGAGCCATTATAAATTGTCCGTTGAGCTTCTACCAGCAAATACCTATAGTGGCGGCACTTCCTGAGCCGGAACGAAAAGTTTTATCGGAATGCGTGTTCTGGTGAACTTTTGGCTTACGGTTGTGATGTTGTGTTG  >  NC_000913/2942481‑2942775

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑