Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,325,980 A→T 16.0% Y412F (TAC→TTC)  yciQ → enhancer of membrane protein expression; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,325,9800AT16.0% 125.1 / 2.2 58Y412F (TAC→TTC) yciQenhancer of membrane protein expression; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (10/38);  new base T (8/2);  total (18/40)
Fisher's exact test for biased strand distribution p-value = 6.92e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTAT  >  NC_000913/1325754‑1326201
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
ggggcgtcagaaaaTAACTCGCTATTCCGGCCCGAAGAGCAGTTATGACGGGAAGGGTTTCTGCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTGACGCAGGGCATTTTTACGGATGCAGAAATTTTTTCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGTGGTGGACTGGCAATATTGGCGGTCATTCATTACg                                                                                                                                                                                                                               <  2:563797/218‑1 (MQ=255)
  cgcgACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGt                                                                                                                                                                                                                             <  2:102499/229‑1 (MQ=255)
       ccAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTAGTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACtt                                                                                                                                                                                                                        <  1:377261/229‑1 (MQ=255)
                 tGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACTGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATCCTTTTTTGCACTGGGGTGGACTTGCAACATTGGCGGTCCTTTTTTTCTGTTACTTTTGTccgcta                                                                                                                                                                                                              >  2:678818/1‑225 (MQ=255)
                 tGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCCTTCTTTATGGTAACTTGAGTACTGTa                                                                                                                                                                                                              >  1:431849/1‑229 (MQ=255)
                   gggTATTAGCTCCCACAGAGCAGTTATTTCGGGGAACGTTTTTTCGTAAAGGCGATAAGGTCGTTATGACGGGGAGACGCATCAGCGTCGTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTTCCGGGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAgg                                                                                                                                                                                                            <  2:488943/228‑1 (MQ=255)
                   cGGTATTACTGCCCGCAGAGCAGTTAGTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAgg                                                                                                                                                                                                            <  1:458476/229‑1 (MQ=255)
                            tGCCCGCAGGGGTGTTAATACCGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCGTAAGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAgg                                                                                                                                                                                                  <  1:333427/229‑1 (MQ=255)
                              cccGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGAGGACTGGCAATATTCGCGGTTGTTCTTTACTGTCAATTTAGTCGCGTTGG‑ATggggatggg                                                                                                                                                                                                >  2:513856/1‑223 (MQ=255)
                               ccGCAGAGCAGTTAATACTGGAAGGGTTATTTCGTGAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTGGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGaa                                                                                                                                                                                                <  1:294422/228‑1 (MQ=255)
                               cagCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGaa                                                                                                                                                                                                <  1:392196/226‑1 (MQ=255)
                                cGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTTTTTACTGTAACTTGAGTTCTGTTGG‑CTTGGCAGGagag                                                                                                                                                                                              >  1:691082/1‑226 (MQ=255)
                                              taCTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGGAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTTTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCg                                                                                                                                                                                 <  1:383777/228‑1 (MQ=255)
                                               tcTGGCAGCGTTATTTCGGAAGGGTGATTAGGGCGTTCTTACGGGGAGACGCAACGAGGTTGTACGAAGGGGATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATTGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGa                                                                                                                                                                                <  2:34773/227‑1 (MQ=255)
                                               aCTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGc                                                                                                                                                                                                     <  1:499758/207‑1 (MQ=255)
                                                cTGGAAGGGTTTGTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAATTTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACGGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGat                                                                                                                                                                               <  1:638718/228‑1 (MQ=255)
                                                cTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGat                                                                                                                                                                               <  1:199583/228‑1 (MQ=255)
                                                    acgctttATTCTGTAAAGGCGATCAGGTCGGTCTTACGGGGAGACGCAACCGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGACGGGTAAAAAGTTTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGt                                                                                                                                                                           <  2:406366/223‑1 (MQ=255)
                                                    aaGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTggg                                                                                                                                                                                                      <  1:220851/201‑1 (MQ=255)
                                                         tttttttGTAAAGGGCATAAGTGTGTTTTGTCGGGGGGACGAAACCGTGTTTTTCCCAGGGCATTTTGACGGCTGCAAAAATTTTATCTGCCGCGTAAAAAGTCTTCATTTGACCGGTCTGATAC‑‑TGTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGTATGATTCTGGTGGGCGATATGTTTATTATg                                                                                                                                                                   <  2:255530/226‑1 (MQ=255)
                                                          ttttggttAAAGG‑GATTAGGGCTTTTTTACGGGTGGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCATAAATTTTATCGGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGTTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATg                                                                                                                                                                   <  2:250114/222‑1 (MQ=255)
                                                           ttttCCTAAAGGCGAAAAGGTCGTTCTTGCGGGGGAACTCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATg                                                                                                                                                                   <  2:431849/228‑1 (MQ=255)
                                                              tattAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGAt                                                                                                                                                                 <  1:726699/225‑1 (MQ=255)
                                                               cGTAAAGGCGATAAGGTCGTTCTTACGGGGCGACGCAACAGAGTCTTACGGAGGGGATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTC‑TTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATctg                                                                                                                                                              <  1:480942/229‑1 (MQ=255)
                                                                                         ggggAGACGCAACAGAGTTTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGGGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTg                                                                                                                                     <  2:685931/228‑1 (MQ=255)
                                                                                             agaCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCCTTCTTTACGGTAACTTTAGTCCTGTAGG‑CTGGGCTGGGATGGGGCTGGTGGGGGGTATGTTTATTTTGGTCTTCTGGGTTATTTCCTTTTTTTTTTTTTTcc                                                                                                                                 >  1:431217/1‑229 (MQ=255)
                                                                                             agaCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACCGTAAATTTAGTCCCGTGAG‑CTGGGCAGGAATGAGTTTCGGGGGCGCTATGGTTATTATTGTCTGGTGGATTT‑TTCCTTTTTTATTTTTTTgcc                                                                                                                                >  1:722660/1‑229 (MQ=255)
                                                                                             agaCGCAACAGAAGCTTACGCAGGGCATTTTGACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTGTGCAATGGGGTGGACTGGCAATATTGGCGGGCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTccc                                                                                                                                <  1:625808/229‑1 (MQ=255)
                                                                                               aCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGGCATCCTTTTCGGTGAATTGAGGCCTGTAGGTCTGGG‑AGGGCTGAGTTTGGGGGGCGCGATGGTTATTTTGAGTTGGTGTGTTT‑TTCATTTTTTTTTTTTTTTCCtt                                                                                                                              >  2:522511/1‑229 (MQ=255)
                                                                                                             tACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTTAGTCCTGTAGG‑CTGGGCAGGAATTAGTCTGGTTGGTGATATGTTTATTATGTTCTGTCGGATTT‑TTCCTTTTTTTTTTTTTTCCCTTTTGCCttttttttt                                                                                                                >  1:215317/1‑227 (MQ=255)
                                                                                                               cGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTTCGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGGGATATGTTTATTTTGATCTTCTGGGTTA‑TTCCTTTTTTTTTTTTTTCCCTTGTGCtttttgtttccc                                                                                                              >  1:221959/1‑221 (MQ=255)
                                                                                                                  aGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGGCATTTTTTTTCGGAAATTTAGTTCTGTAGG‑CTGGGCAGGGATGAGTCTTGTGGGGGATTTGTTTTTTTTTTTTGTTTGGGTTA‑TTTATTTTTTTTTTTTTTTCCTTTTGGtttttttttttcccg                                                                                                           >  2:664335/1‑222 (MQ=255)
                                                                                                                    ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAACGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTTATTTTTTATTGTTACTTGAGTCCTTTGGG‑CTGGGCAGGGATGAGTTCTGTGGGGGATTTTTTTATTTTGGTCTTCCTGGTTA‑TTTCTTTTTTGTTTTGTTTCCCTGAGATTTTTTTTTCCcgtgtt                                                                                                         >  1:359071/1‑225 (MQ=255)
                                                                                                                    ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCCATCTTTTCTGTGACTTGAGGCCTGTAGG‑CTGGGCAGGGATGGGGCTGGTTGGCGATATGTTTTTTTTGTGTTTGTGGGGTA‑TTCCTTTTTTTTTTTTTTCCCTTTTGTTtttttttgtcgcgc                                                                                                           >  2:115767/1‑221 (MQ=255)
                                                                                                                    ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTTAGTCCTGTAGG‑CTGGGGAGGAATGAGTCTGGTGGGCGGTATTTTTATTTTGTTCTTCTGGGTTT‑TTCCTTTTTTTTTTTTTTTCCTTTGGCTTTTTTTTGCCcccgg                                                                                                          >  1:161042/1‑227 (MQ=255)
                                                                                                                    ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTTAGGCCTGTTGG‑GTGGGCAGGAAATAATCTGGGGGGCGGTATGTTTATTATTGTTTTGTGGGTTATTTCTTTTTTTATTTTTTTCCTTTGGCTTTTTGTTGGccccgg                                                                                                           >  1:697396/1‑225 (MQ=255)
                                                                                                                         ttttTCGGATGCAGAATTTTTTTCTGCCGGTTAAAAAGTCTTCGTTTTACCGGTCTGATTCTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGAtgatga                                                                                                    <  2:221959/229‑1 (MQ=255)
                                                                                                                          ttttCGGGTGCAAAAATTTTATCTGCCGGGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGAAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGAc                                                                                                   <  2:522977/229‑1 (MQ=255)
                                                                                                                          cttACGGATGCATAAATTTTAGCTGTCGCTTAAAAAGTTTTCGTTTTACCGGTCTTATACTTTTTTGCAATGGGGTGGACTTGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGA‑CTTGGCAGGAATGAGTCTGGTGGGCTATATGTTTATTATTATCTGCTGTATTA‑TTCCTTTTTTATTTTGTTCCCTTGTGCTTTTGTTTGCCCGCGatgatgtc                                                                                                   <  1:667230/228‑3 (MQ=255)
                                                                                                                                 aTGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTAGCGGTCTGATTCTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCt                                                                                             <  1:743399/228‑1 (MQ=255)
                                                                                                                                            tttttcggcgtgTAAAAAGTTTTTGTTTTACGGGTCTGATGCTTTTTGGCAATGGGGTGGACTGGCAATAGTGGTGGTCATTCTTTACGGTAGTTTGAGTCCTGTAGT‑CTGCGGCGGAATTAGTCTTGTAGGTGTTATGTTTATTATGTTCTGATTGATTA‑TTCCTTTCTTATTTTGTTCCCTTGACCTTTTGTTTGCCTTCGATGATGACAAGCCTTGCGTTAATCgt                                                                                 <  2:163285/218‑1 (MQ=255)
                                                                                                                                                    cgcgTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTTCTGGAACTTGAGTCCTGTAGG‑CGGGGCGGGGATGAGGTTGGTGGGCGATATGTTTTTTTTTGTCTTCGTGGTTT‑TTCCTTTTTTATTTTTTTCCCTTTAGGTTTTGTTTGTCCCCGGTGATTAAACACCTTTCGTTTATTGTTTatttat                                                                         >  2:13818/1‑224 (MQ=255)
                                                                                                                                                               ttttCGTTTTACCGGTCTGATTCTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGtttt                                                              <  1:664336/227‑1 (MQ=255)
                                                                                                                                                                 ttCGTTTTACCGGTCTGATTCTTTTTTGCAAGGGGGTGGACTTGAAATAGTGGCGGTTATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGttttt                                                             <  1:678819/228‑1 (MQ=255)
                                                                                                                                                                 ttCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTTTGGGAAATTGGGTCCTGGGTG‑GTTGGGAGGAATGAGTATGGTGGGCGATTTGTTTTTTTTGATCTTGTGGGTTA‑TACCCTTTTTTTTTTGTTTCCCTGTGCTTTTGTTTGTCCGCGTTGATTGCAAGTGTTTTGCTTAACGTGTAATACTACCTTTGTTTTTa                                                            >  2:339540/1‑229 (MQ=255)
                                                                                                                                                                          cctgttttatatTTTTTTGGAATGGGGGGGACTGGCAATTTGGGGGGTCATTCTTTACGGTAACTTGAGTCCGGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGAt                                                    <  2:215317/220‑1 (MQ=255)
                                                                                                                                                                                      ttttttGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGGAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTg                                       <  2:159604/229‑1 (MQ=255)
                                                                                                                                                                                        ttttGCAATGGGGTGGACTGGCAATATTGGCGGTTAATCTTTTCGGTTACTTTGGTCCTGTAGG‑GCGGGGAGGGATAGGTCTGGGTGGCGGTATGTTTTTTTTTATTTGCTGGGTTA‑TTTCTTTTTTTTTTTTTTTCCTTTAGGTTTTTTTTTCCCGCGAGTATGGAAAGCCTTTCTTTTATCTATGAATTATTCCTTTGTTTTTTTCACCGCTTTTTTCatggggg                                       >  2:724709/1‑221 (MQ=255)
                                                                                                                                                                                        ttttGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAAATTGAGTTCTGTAGG‑CTGGGCAGGAATGAGGTTGGTGGGCGATATGTTTATTATGATCTGCTTGATTG‑TTCCTTTTTTATTTTGTTCCCTTGTGCTTTTGTTTGCCCGCGTTGATGTCAAGCCTTTCGTTTATCGTGTAATCCTCACTTTGTTTTTTCACCTGATTTGTTCAGTCGTGGc                                     >  1:746826/1‑229 (MQ=255)
                                                                                                                                                                                           tGCAATGGGGTGGACTGGGAATATTGGTGGTTATTCATTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGTAATGAGTATGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCTTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCtt                                  <  1:320508/229‑1 (MQ=255)
                                                                                                                                                                                                tGGGGTGGACTGGCAATTTTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTAtt                             <  1:522511/229‑1 (MQ=255)
                                                                                                                                                                                                  gggTGGTCTGGCAATTTTGGCGGTTGTTTTTTTGGGTAACTTGAGTCCTGTTGG‑CTGGGCAGGAATTGGTCTTGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATtct                           <  2:359071/229‑1 (MQ=255)
                                                                                                                                                                                                  gggTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑TTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATtct                           <  1:534385/229‑1 (MQ=255)
                                                                                                                                                                                                          tGGAAATATTGGGGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTTTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACAt                   <  1:11208/229‑1 (MQ=255)
                                                                                                                                                                                                                   tGGCGGTCATTCTTTACGGTAACTTGGTTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGTGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCg           <  1:513856/228‑1 (MQ=255)
                                                                                                                                                                                                                   tGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTTGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTGAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCg           <  1:149521/228‑1 (MQ=255)
                                                                                                                                                                                                                          taTTCTTTACGGTAATTTGAGTCCTTGAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTa   <  2:161042/228‑1 (MQ=255)
                                                                                                                                                                                                                           aTTCTTGTCGGTACCTTTAGTTCGGTAGG‑CTTGGCAGGAAGGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑GTCCTTTTTTATTTTGTTCCCTTGAGCTTTTTTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCCCTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTAt  <  1:596927/229‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
GGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGG‑CTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTA‑TTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTAT  >  NC_000913/1325754‑1326201

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