Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,050,717 A→T 5.1% intergenic (‑50/+398) gcvT ← / ← ubiI aminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex/2‑octaprenylphenol hydroxylase, FAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,050,7170AT 5.1% 390.8 / 6.2 129intergenic (‑50/+398)gcvT/ubiIaminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex/2‑octaprenylphenol hydroxylase, FAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (61/60);  new base T (5/2);  total (66/63)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

TTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGCGAC  >  NC_000913/3050568‑3050865
                                                                                                                                                     |                                                                                                                                                     
            aaTGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGATAATTTCACCATGaaaa                                                                                                                                           <  1:1732868/150‑2 (MQ=255)
               gCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCACAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTACTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATGACAACGATAcccc                                                                                                                                        >  1:94163/1‑149 (MQ=255)
               gCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAACCCCa                                                                                                                                        >  2:1952880/1‑150 (MQ=255)
                cAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAg                                                                                                                                       >  1:1156602/1‑150 (MQ=255)
                 aGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCCCGCAACCCAGc                                                                                                                                      >  2:276516/1‑150 (MQ=255)
                 aGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGGGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCGTGGTGATAATTTCACCATGaaaaagttg                                                                                                                                      <  1:1577997/150‑7 (MQ=255)
                     gCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCTGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCg                                                                                                                                  <  1:1096919/150‑1 (MQ=255)
                     gCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCCGCATCg                                                                                                                                  >  2:1695991/1‑150 (MQ=255)
                        caagccATCCGTGGATATACACCATGCGAGCGCCGCAAATCGGGGGTTGTGCGTACAAAGGAGTCTGTTGTGACATATTGTCCTCATTGAATAAGCGTGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAAACCAGCATCGGAg                                                                                                                               <  2:825472/146‑1 (MQ=255)
                          atcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGcc                                                                                                                             <  1:212026/150‑1 (MQ=255)
                          atcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGcc                                                                                                                             <  1:910223/150‑1 (MQ=255)
                          atcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGATGTGCCATCTTGTCCTCATTGAATAAGCGCGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGcc                                                                                                                             >  2:1677819/1‑150 (MQ=255)
                            ccAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACTAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCAc                                                                                                                           <  2:701137/150‑1 (MQ=255)
                             cAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACTCAGCATCGGAGCCACt                                                                                                                          >  2:1188582/1‑150 (MQ=255)
                                ccGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCCCTccc                                                                                                                       >  1:672410/1‑150 (MQ=255)
                                ccGTGGAAATCCACCATGCGAGCGCCGCAAAGCGAGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTCTCCTCATAGAATAAGCGGGGCTGACAACTTTTTCATGGTCAAATTATCACCCCGAATCCCAGCATCGGAGCCACTccc                                                                                                                       >  1:33975/1‑150 (MQ=255)
                                 cGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCCGCATCGGAGCCACTCCCg                                                                                                                      >  2:639050/1‑150 (MQ=255)
                                  gTGGAACTCCACCATGCGAGCGCCGCAAAGCGGGTGTTGGTCGTACAAAGGAGTCTGTTGGGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCgg                                                                                                                     <  2:1231228/150‑1 (MQ=255)
                                  gTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCgg                                                                                                                     >  2:906334/1‑150 (MQ=255)
                                   tGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGCAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATAATCACCACGAAACCCAGCATCGGAGCCACACCCGGt                                                                                                                    >  1:1551991/1‑150 (MQ=255)
                                    ggAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTc                                                                                                                   <  1:1931510/150‑1 (MQ=255)
                                     gAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTcc                                                                                                                  <  1:1388190/150‑1 (MQ=255)
                                      aaaTCCACCATGCGAGCGCCGCAAAGCGGGGGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTccc                                                                                                                 <  2:977602/150‑1 (MQ=255)
                                      aaaTCCACCATGCGAGCGCCGCAAAACGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTccc                                                                                                                 <  2:139099/150‑1 (MQ=255)
                                       aaTCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCCCTCCCGGTcccc                                                                                                                >  1:633489/1‑150 (MQ=255)
                                        aTCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCa                                                                                                               <  1:1642005/150‑1 (MQ=255)
                                        aTCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCa                                                                                                               >  1:529432/1‑150 (MQ=255)
                                          ccaccaTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGTTGCTGACATCTTTTTCATGGTGAAATTATCACCACGAAACGCAGCATCGGAGCCACTCCCGGTCCCCAAc                                                                                                             <  2:1162172/150‑1 (MQ=255)
                                          ccaccaTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAAc                                                                                                             <  2:1185103/150‑1 (MQ=255)
                                           caccaTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                                                            >  1:747406/1‑150 (MQ=255)
                                               aTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTAGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAAt                                                                                                        <  2:1102954/150‑1 (MQ=255)
                                                tGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATc                                                                                                       >  2:1153374/1‑150 (MQ=255)
                                                 gCGACCGCCGCAAAGCGGGTGTTGTTCGTACAAAGGAGGCTGTTGTGCCATCTTGGCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCg                                                                                                      <  1:1142692/150‑1 (MQ=255)
                                                    aGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTtc                                                                                                   <  1:1230653/150‑1 (MQ=255)
                                                      ctccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTTTTGTGCCATCTTGTCCTCATTGAATAAGCGGGTCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTtctc                                                                                                 <  1:1239168/148‑1 (MQ=255)
                                                       gccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTctct                                                                                                >  2:103251/1‑150 (MQ=255)
                                                        ccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtt                                                                                               >  1:1814598/1‑150 (MQ=255)
                                                        ccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtt                                                                                               <  1:1804346/150‑1 (MQ=255)
                                                         cgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCACCGGAATCGTTCTCttt                                                                                              >  2:157066/1‑150 (MQ=255)
                                                         cgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCACCGCAATCGTTCTCttt                                                                                              >  2:878446/1‑150 (MQ=255)
                                                          gcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGAGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtttt                                                                                             >  1:597038/1‑150 (MQ=255)
                                                           ctaattgtgtgtTTTTCGTACAATGGAGTCTGTTGTGCCATGTTTTCCTCATTGAATAAGAGGGGCTGTCTACTTTTTCATGGTGAAATTATCACCAGGAATCCCAGCATCGGATCCACTCCAGGTCCCCAACGCAATCGTTCTCTTTTg                                                                                            <  1:2009130/144‑1 (MQ=255)
                                                             aaGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGcc                                                                                          <  2:1146763/150‑1 (MQ=255)
                                                             aaGCGTGTGTTGGTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGcc                                                                                          <  1:978394/150‑1 (MQ=255)
                                                              aGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCt                                                                                         <  2:1274221/150‑1 (MQ=255)
                                                                 gtgtgtTGTTCGTACTTAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAATTTATCACCACGAAACCCAGGATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGaa                                                                                      <  2:1619499/150‑1 (MQ=255)
                                                                  tgtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGTGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGTGCCTCTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAAc                                                                                     <  1:1040964/150‑1 (MQ=255)
                                                                  tgtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAAc                                                                                     <  1:85193/150‑1 (MQ=255)
                                                                  tgtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAAc                                                                                     <  2:230400/150‑1 (MQ=255)
                                                                  tgtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAAc                                                                                     <  1:1972165/150‑1 (MQ=255)
                                                                     gttgttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTa                                                                                  >  2:1586562/1‑150 (MQ=255)
                                                                        gttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCCCCACGAAACACAGCATCGGCGCCACTCACGGTCCCCAACTCACTCGTTCTCTTTTGCATGAACTTcccc                                                                               >  1:1223161/1‑149 (MQ=255)
                                                                         ttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTaccac                                                                              <  2:758613/150‑1 (MQ=255)
                                                                         ttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTaccac                                                                              >  1:1149890/1‑150 (MQ=255)
                                                                          tCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCTTCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTATGCCTGAACTTaccacc                                                                             >  2:488916/1‑150 (MQ=255)
                                                                            gTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGa                                                                           >  2:794543/1‑150 (MQ=255)
                                                                              aCAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCATGAACTTACCACCGaaa                                                                         >  2:1481540/1‑150 (MQ=255)
                                                                                 aaGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAg                                                                      >  1:1953454/1‑150 (MQ=255)
                                                                                   ggAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGAc                                                                    <  1:1313270/150‑1 (MQ=255)
                                                                                      gTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGt                                                                 >  1:1670880/1‑150 (MQ=255)
                                                                                          gttgtGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAc                                                             <  1:775626/150‑1 (MQ=255)
                                                                                           ttgtGCCATCTAGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAAATATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAcc                                                            <  2:1850488/150‑1 (MQ=255)
                                                                                            tgtgCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCa                                                           <  1:1631438/150‑1 (MQ=255)
                                                                                                   tCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTa                                                    >  2:1287598/1‑150 (MQ=255)
                                                                                                        tCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAAtt                                               <  2:1814598/150‑1 (MQ=255)
                                                                                                        tCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAAtt                                               >  2:570099/1‑150 (MQ=255)
                                                                                                          cTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGa                                             >  2:1484788/1‑150 (MQ=255)
                                                                                                          cTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGa                                             >  2:1246682/1‑150 (MQ=255)
                                                                                                          cTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATAACCACGAAACCCAGCATCGGAGCCACTCACGGTCCCCAACGCAATCGTTCTATTTTGCCTGAACTTACCACCGAAACAGACAGTTAACCATAAGGTAAAATTGa                                             >  1:1861408/1‑150 (MQ=255)
                                                                                                           acATTGAATAAGCTGGGCGGACTACTTTTTAATGGTGAAATTCTCACCACGAAACCCAGCATCGGATCCACTCCCGGGACCCACCGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGat                                            <  2:355610/149‑1 (MQ=255)
                                                                                                            cATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatc                                           <  2:1473679/150‑1 (MQ=255)
                                                                                                             aTTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATAACCACGAAACCCAGCATCGGAGCCACTCACGGTCCCCAACGCAAACGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTATACCATAAGGTAAAATTGatca                                          >  1:1226868/1‑150 (MQ=255)
                                                                                                                gAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACCGACTGTTAACCATAAGGTAAAATTGAtcatca                                      <  2:1037509/150‑1 (MQ=255)
                                                                                                                gAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGAtcatca                                      >  1:429691/1‑150 (MQ=255)
                                                                                                                  aTAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATcaca                                    >  1:957245/1‑150 (MQ=255)
                                                                                                                   tAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACAt                                   >  2:1056242/1‑150 (MQ=255)
                                                                                                                   tAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGGTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACAt                                   >  1:964726/1‑150 (MQ=255)
                                                                                                                     aGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTa                                 >  2:1544481/1‑150 (MQ=255)
                                                                                                                      gCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAg                                <  1:909198/150‑1 (MQ=255)
                                                                                                                       cGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGc                               >  1:458728/1‑150 (MQ=255)
                                                                                                                        ggggCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGACACTCCCGGTCCCCAACGCAATCGTTCTCATTTGCCTGAACTAACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATAACATTAGCt                              >  2:1427396/1‑150 (MQ=255)
                                                                                                                            cTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATAGTTCTCTTTTGCCTGAACTAACCACCCAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTAtt                          >  2:1554178/1‑149 (MQ=255)
                                                                                                                             tGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATgg                         <  2:1136895/150‑1 (MQ=255)
                                                                                                                             tGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCAAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATgg                         <  2:1746335/150‑1 (MQ=255)
                                                                                                                               tgctgggtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGtt                       >  2:491025/8‑150 (MQ=255)
                                                                                                                                cAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCTGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTa                      <  1:169920/150‑1 (MQ=255)
                                                                                                                                  tgggtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaa                    >  1:588921/5‑150 (MQ=255)
                                                                                                                                  aCTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaa                    >  2:1912529/1‑150 (MQ=255)
                                                                                                                                   cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCTTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGTTTaaaa                   >  2:280340/1‑150 (MQ=255)
                                                                                                                                   cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaa                   >  2:1764092/1‑150 (MQ=255)
                                                                                                                                   cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaa                   >  1:1367180/1‑150 (MQ=255)
                                                                                                                                    tttttCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGCTAATTTGATGATCAATTTTACATTATGGTTaacag                  <  1:1449273/150‑4 (MQ=255)
                                                                                                                                     gtttCGTGGTGATAATTTCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACGACCGAAACAGACGGTTAACAATAAGGTAAAATTGATCATCACTTTAGCTTATGGTTaaaaaa                 >  1:1024830/2‑150 (MQ=255)
                                                                                                                                       ttCATTGTGAAAATATCAAAACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTTAACTTACCACCGAAACCGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGGTAAAAAATg               <  2:1718548/150‑1 (MQ=255)
                                                                                                                                         cATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCa             <  2:529432/150‑1 (MQ=255)
                                                                                                                                           tGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACTCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaa           <  1:1433592/150‑1 (MQ=255)
                                                                                                                                            ggTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaa          >  2:1912831/1‑150 (MQ=255)
                                                                                                                                               gAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATc       <  1:1056242/150‑1 (MQ=255)
                                                                                                                                                  aTTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATcgcg    <  1:1695991/150‑1 (MQ=255)
                                                                                                                                                    tATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGCGAc  <  1:1246682/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGCGAC  >  NC_000913/3050568‑3050865

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: