Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 32.7% intergenic (+266/+376) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT32.7% 28.1 / 41.7 49intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (21/12);  new base T (6/10);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-01

ATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGC  >  NC_000913/4295840‑4296273
                                                                                                                                                                                                                             |                                                                                                                                                                                                                      
aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAACGTG‑ATTTGGTAGTTCGGATAAGGCGTTTACGCCGCATCCGCACTCAATGCCTgt                                                                                                                                                                                                                  >  1:573198/1‑226 (MQ=255)
           gCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCACCATCAATGCCCGAGCGGACGCTTGc                                                                                                                                                                                                      >  2:646170/1‑228 (MQ=255)
                     aTCCGGCCTGCGCAAGCCAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                            <  1:496851/228‑1 (MQ=255)
                     aTCAGGACTACGCAGGACAGCGGAATAGCCTGATTTAGCGTGATTTTGGCGGTCGGATAAGGGATTTATGCCGCATCCGACATCAACGCCTGATGCGCAGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                            <  1:484310/228‑1 (MQ=11)
                          gccaccGCCAGGCAGGGCAATAGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                       <  1:164497/223‑1 (MQ=255)
                           ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGGATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                       <  1:346087/227‑1 (MQ=18)
                               cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATACGCAGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                                                                                  >  1:399114/1‑226 (MQ=12)
                               cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGAATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTTTCAGGCCTACg                                                                                                                                                                                  >  1:342401/1‑228 (MQ=21)
                                    agatgcgcAATTGCCTGATTTTGGGGGATTTTGGAGGTCGGATAAGGCTGTTATGCCGCATCCGACACCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                             <  1:68117/224‑1 (MQ=11)
                                       ggcgcAATAGCCTGAGTTAGCGTGATTTTGGAGGTCTGAAAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTAGGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                                                                                                                                                                           <  1:214451/226‑1 (MQ=255)
                                                 ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                                                                                                   >  1:230552/1‑193 (MQ=34)
                                                  cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                                                                                                                                                                                  <  2:230552/193‑1 (MQ=34)
                                                               aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACAGCGCAAAAGCCTGATTTTGCGTGAt                                                                                                                                                  >  1:403015/1‑228 (MQ=21)
                                                                          cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCCTCGCAATAACTTTTTTTAACGGCATTTTGtgtttcg                                                                                                                                       >  2:596102/1‑222 (MQ=11)
                                                                          cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                                                                                        >  2:125538/1‑179 (MQ=21)
                                                                          cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                                                                                                                                                                                  >  1:366905/1‑169 (MQ=35)
                                                                          cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCACGACAACGAAATAGCCTGATTTAGCGTGATTTTGTTGGtgc                                                                                                                                       >  1:619717/1‑226 (MQ=255)
                                                                           ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                       <  1:125538/179‑1 (MQ=21)
                                                                           ggATAAGGCGTTTATGCCGCATCCGACATCAACGCATGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                 <  2:366905/169‑1 (MQ=25)
                                                                             aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACGTCAATGCCTGATGCGACGCTTGCCGCGCCTTACCAGGCCTACCCCAGCAAGGGCAAACGCCTCATTTAGCTTAATTTGGtcgaccgga                                                                                                                                     >  2:522489/1‑219 (MQ=11)
                                                                                                tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                 <  2:142947/148‑1 (MQ=25)
                                                                                                tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                 >  1:142947/1‑148 (MQ=37)
                                                                                                  cggcATCAACGCCGGATGCGACGTTTAACGCGTCTTAGCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTT‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCg                                                                                                               <  1:32237/225‑1 (MQ=255)
                                                                                                               tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                                                                                                   <  1:26733/147‑2 (MQ=21)
                                                                                                               tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                                                                                                   >  2:26733/1‑146 (MQ=37)
                                                                                                                 aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGACTCCGCCATCAACACCTGa                                                                                                 >  1:241933/1‑227 (MQ=255)
                                                                                                                                 gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                      <  2:61694/126‑1 (MQ=25)
                                                                                                                                 gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                      >  1:61694/1‑126 (MQ=38)
                                                                                                                                  tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAc                                                                                                                                                                                                            >  2:281479/1‑103 (MQ=32)
                                                                                                                                  gcTTATCAGGCCTAGACCAGGCAGCGTAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑Gt                                                                               <  1:30101/226‑1 (MQ=19)
                                                                                                                                   cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACg                                                                                                                                                                                                           <  1:281479/103‑1 (MQ=32)
                                                                                                                                           gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATCGCCTGATTTCGCGTGATTTTGTGTGTCGGATAAGGCGTTTACGCCGCATCCGTCCTCAACACCTGCTGCGCCCCTTCACGC‑GTCTTTTCAgtc                                                                     >  2:519142/1‑226 (MQ=11)
                                                                                                                                           ctcttgacaaGACAGCCCAATAGCCTTATTTAGCGTG‑AGTTTGTAGGTCGGATAAGGCGGTTACGCCGGATCCGACATCAACGCCTGATGCGGCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                     <  1:161651/219‑1 (MQ=11)
                                                                                                                                            gttaagtcggACAGCGCAATAGCTTGATTT‑GGGTGTATTGTGTGGGACGGATAAGGCGTTTAAGCCGGGTACGACATCAACGCCGGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                     <  2:288570/218‑1 (MQ=11)
                                                                                                                                                cggcAGGCAGAGCAATAGCCTGATTTAGCGTG‑ATTTTGGAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACg                                                                <  1:193671/225‑1 (MQ=255)
                                                                                                                                                cGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGGCGGATAAGCCGTTTATGCCGCATCCGCCCTCAACGCCGAGTCCGACGCTTACCGC‑ATCTTACCAGCACTACg                                                                >  1:462237/1‑228 (MQ=11)
                                                                                                                                                cGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGt                                                                                                                                         <  1:74792/156‑1 (MQ=35)
                                                                                                                                                      aCAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGAATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTCGCGCGTTCTTGTCGGTCGGCTAAGGCGCTTATGCCGCATCCGACTTCAACGCCCAATGAGACGCTTTACCA‑CTCTATTCACGCCTACACCacgc                                                          >  2:161651/1‑225 (MQ=11)
                                                                                                                                                          cgcAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                                                                                                      <  2:307593/85‑1 (MQ=37)
                                                                                                                                                          cgcAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTACGCCAGCCATCGCAATATCCTGATTTCGCGTCATTTTGTAGGTCGGATAAGGCGTTTATGCCCCATCCGACATCAACGTCTGATGCGACGCTT‑ACGCTGCCTTACCGGGCCAACGCCAGAcccgc                                                      >  2:297567/1‑224 (MQ=11)
                                                                                                                                                          cgcAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                                                                                      >  1:307593/1‑85 (MQ=37)
                                                                                                                                                                      aTTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCAGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGAATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GCCTACTCAGGCCTACGCCAGCCAGCGCAATAGCCTGAt                                          >  1:275609/1‑228 (MQ=11)
                                                                                                                                                                                 aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                                                                                                      <  1:629770/63‑1 (MQ=255)
                                                                                                                                                                                 aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                                                                                      >  2:629770/1‑63 (MQ=39)
                                                                                                                                                                                            cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcg                                                                                                                                                                                                    >  2:414895/1‑54 (MQ=255)
                                                                                                                                                                                             ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                                                                                                   <  1:414895/54‑1 (MQ=255)
                                                                                                                                                                                                       tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTACCGGTATTTTTCAAGTCGGATAAgggatt          >  2:390828/1‑224 (MQ=11)
                                                                                                                                                                                                         tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAATGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATACCCTGATTTAGCGGGATTTTGTAGGTCGACAAAGGCGATTa        >  2:514447/1‑228 (MQ=11)
                                                                                                                                                                                                               ttaTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGACGCTTACCGC‑GTCTTATCAGGCCTACTCCAGACAGCGCAAAAGCTCGTATTAGCGTGATTTTGCTGGTCGGATCAGGCCTTTaacccca  >  2:242908/3‑222 (MQ=11)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                      
ATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTG‑ATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGC  >  NC_000913/4295840‑4296273

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