Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,236,658 C→A 7.0% intergenic (+158/+85) preA → / ← mglC dihydropyrimidine dehydrogenase, NADH‑dependent, subunit C/methyl‑galactoside transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,236,6580CA 7.0% 234.5 / 6.4 71intergenic (+158/+85)preA/mglCdihydropyrimidine dehydrogenase, NADH‑dependent, subunit C/methyl‑galactoside transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base A (3/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-01

GCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTG  >  NC_000913/2236511‑2236806
                                                                                                                                                   |                                                                                                                                                    
gCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAAGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATgcg                                                                                                                                                     >  4:179020/1‑149 (MQ=255)
  aTCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATgcggc                                                                                                                                                   <  2:356084/149‑1 (MQ=255)
   tCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATgcggcg                                                                                                                                                  >  4:232825/1‑149 (MQ=255)
   tCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATgcggcg                                                                                                                                                  >  6:413711/1‑149 (MQ=255)
     cGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCCGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTg                                                                                                                                                <  7:417807/149‑1 (MQ=255)
              gTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCATCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCATTa                                                                                                                                       >  8:163851/1‑149 (MQ=255)
                 ccTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                     <  8:273412/148‑1 (MQ=255)
                  cTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCTCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    <  4:84802/148‑1 (MQ=255)
                     tCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                <  3:409918/149‑1 (MQ=255)
                     tCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                >  1:246628/1‑149 (MQ=255)
                       gCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                              <  3:232825/149‑1 (MQ=255)
                       gCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                              >  1:353670/1‑149 (MQ=255)
                            aTGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                                                         <  2:353670/149‑1 (MQ=255)
                                 gCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACcaca                                                                                                                    <  4:363686/149‑1 (MQ=255)
                                         gCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCat                                                                                                            >  7:343883/1‑149 (MQ=255)
                                         gCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCat                                                                                                            >  5:260300/1‑149 (MQ=255)
                                         gCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCat                                                                                                            >  5:396521/1‑149 (MQ=255)
                                           cTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCatga                                                                                                          >  7:1259/1‑149 (MQ=255)
                                              aTCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCatgat                                                                                                         <  3:872/147‑1 (MQ=255)
                                               tCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATAAGGCGTTCACGCCGCATCCGGcaaca                                                                                                                             <  4:163444/126‑5 (MQ=38)
                                               tCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATAAGGCGTTCACGCCGCATCCGGcaaca                                                                                                                             >  3:163444/1‑122 (MQ=38)
                                                 cGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTAt                                                                                                    >  6:372850/1‑149 (MQ=255)
                                                 cGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTAt                                                                                                    <  4:131358/149‑1 (MQ=255)
                                                  ggCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTAtg                                                                                                   >  4:179365/1‑149 (MQ=255)
                                                   gCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTAtgt                                                                                                  >  8:364628/1‑149 (MQ=255)
                                                   gCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTAtgt                                                                                                  >  4:155333/1‑149 (MQ=255)
                                                      tACAAAACCATACAAATCAGTATGTTGAGCTTATTGTATGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTACTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCg                                                                                               >  3:208470/1‑149 (MQ=255)
                                                         aaaaCCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGttt                                                                                            <  5:413711/149‑1 (MQ=255)
                                                         aaaaCCAAACAAAACAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTTCATCAGGCATTTTGATCCTGTTGACGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATAAAGATTTTATGTGCGtg                                                                                             >  1:274030/1‑147 (MQ=37)
                                                               aTACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCatgatg                                                                                                        <  1:24104/131‑1 (MQ=255)
                                                               aTACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCatgatg                                                                                                        >  2:24104/1‑131 (MQ=255)
                                                                  cAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGaa                                                                                   <  7:127253/149‑1 (MQ=255)
                                                                   aaaTCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAAt                                                                                  <  3:362619/149‑1 (MQ=255)
                                                                   aaaTCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAAt                                                                                  >  1:100050/1‑149 (MQ=255)
                                                                     aTCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGc                                                                                >  3:386694/1‑149 (MQ=255)
                                                                        aGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTTGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGa                                                                             <  6:260300/149‑1 (MQ=255)
                                                                        aGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGa                                                                             <  6:396521/149‑1 (MQ=255)
                                                                              ttGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACcaca                                                                        <  4:285740/148‑1 (MQ=255)
                                                                                gAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCTGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACAtt                                                                      <  2:388356/148‑1 (MQ=255)
                                                                                  gCTTATTGTAGGCAGGATAAAACGCTGCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTTTCCGGCCTACAAAACCACACAGATCTTGATTTTATGTGCGTTTATGACCGATTGCGGACCACATTcac                                                                   <  2:274030/149‑1 (MQ=38)
                                                                                   cTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTcac                                                                   <  6:402966/148‑1 (MQ=255)
                                                                                         gTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCAttt                                                            >  2:136233/1‑149 (MQ=255)
                                                                                             gCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTa                                                        >  4:42440/1‑149 (MQ=255)
                                                                                             gCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTa                                                        >  6:112899/1‑149 (MQ=255)
                                                                                                  aTAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTAcgcgc                                                   <  3:179020/149‑1 (MQ=255)
                                                                                                   tAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTca                                                                    <  2:365130/131‑1 (MQ=255)
                                                                                                   tAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTca                                                                    >  1:365130/1‑131 (MQ=255)
                                                                                                     aaaCGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTa                                                >  5:253871/1‑149 (MQ=255)
                                                                                                      aaCGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTAt                                               >  6:71939/1‑149 (MQ=255)
                                                                                                      aaCGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTAt                                               >  8:372771/1‑149 (MQ=255)
                                                                                                       aCGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTAtt                                              <  3:155333/149‑1 (MQ=255)
                                                                                                        cgcgCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTAttt                                             <  2:246628/149‑1 (MQ=255)
                                                                                                           gcCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAg                                          >  7:188032/1‑149 (MQ=255)
                                                                                                             cAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGt                                         >  2:226073/1‑148 (MQ=255)
                                                                                                             cAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGt                                         >  3:55602/1‑148 (MQ=255)
                                                                                                              aGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCgg                                                                              <  7:99492/110‑1 (MQ=255)
                                                                                                              aGCGTTGCATCCGGCATTTGGCTACTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACACATCATGATGTTATGTGCGTTTATGACATAATGCgg                                                                              >  8:99492/1‑110 (MQ=25)
                                                                                                                 gTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGCACCACATTCACATCATTTCTAACGCGCGAATTTCAGTGAATc                                    >  2:348568/1‑149 (MQ=255)
                                                                                                                     cATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAgc                                <  8:451232/149‑1 (MQ=255)
                                                                                                                     cATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGACTACAAAACCAAAAATATCATGATGTTATTTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAgc                                >  4:59105/1‑149 (MQ=39)
                                                                                                                       tCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAgcgc                              >  5:253587/1‑149 (MQ=255)
                                                                                                                           gCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTAtgt                                                                                                  <  4:28585/77‑1 (MQ=255)
                                                                                                                           gCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTAtgt                                                                                                  >  3:28585/1‑77 (MQ=255)
                                                                                                                           gCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTAATCCGGCCTACAAAACCACACAGATCATGAAGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACg                          >  6:33160/1‑149 (MQ=255)
                                                                                                                             aTTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACATATCATGATGTTATGTTCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGc                        <  1:381667/149‑1 (MQ=255)
                                                                                                                               ttGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCg                       >  3:291777/1‑148 (MQ=255)
                                                                                                                                 ggCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGaa                     <  3:42440/148‑1 (MQ=255)
                                                                                                                                 ggCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGaa                     <  7:364628/148‑1 (MQ=255)
                                                                                                                                  gCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGa                   <  7:372771/149‑1 (MQ=255)
                                                                                                                                    tCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGat                  >  1:291339/1‑148 (MQ=255)
                                                                                                                                      cTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGataa                >  1:74436/1‑148 (MQ=255)
                                                                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGata                 <  1:456912/142‑1 (MQ=255)
                                                                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGata                 >  2:456912/4‑145 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATc           >  2:172228/2‑149 (MQ=255)
                                                                                                                                           gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATc           <  4:279955/148‑1 (MQ=255)
                                                                                                                                                   cggcgTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTtg  <  8:1259/149‑1 (MQ=255)
                                                                                                                                                   cggcgTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTtg  <  2:100050/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATCGCCCCTTTG  >  NC_000913/2236511‑2236806

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