Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,236,657 G→A 7.0% intergenic (+157/+86) 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,6570GA 7.0% 215.5 / 4.9 71intergenic (+157/+86)preA/mglCdihydropyrimidine dehydrogenase, NADH‑dependent, subunit C/methyl‑galactoside transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (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.26e-01

GCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATC  >  NC_000913/2236511‑2236797
                                                                                                                                                  |                                                                                                                                            
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            <  7:364628/148‑1 (MQ=255)
                                                                                                                                 ggCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGaa            <  3:42440/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)
                                                                                                                                                  |                                                                                                                                            
GCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCATACAAATCAGTATGTTGAGCTTATTGTAGGCCGGATAAAACGCGCCAGCGTTGCATCCGGCATTTGGCTCCTGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCACACAGATCATGATGTTATGTGCGTTTATGACCGAATGCGGACCACATTCACATCATTTCTTACGCGCGTATTTCAGTGAATCCAGCGCTACGGCGAAGATAATAATC  >  NC_000913/2236511‑2236797

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: