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

GGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTAT  >  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)
                            tGCCCGCAGGGGTGTTAATACCGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCGTAAGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAgg                                                                                                                                                                                                  <  1:333427/229‑1 (MQ=255)
                              cccGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGAGGACTGGCAATATTCGCGGTTGTTCTTTACTGTCAATTTAGTCGCGTTGGATggggatggg                                                                                                                                                                                                >  2:513856/1‑223 (MQ=255)
                               ccGCAGAGCAGTTAATACTGGAAGGGTTATTTCGTGAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTGGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGaa                                                                                                                                                                                                <  1:294422/228‑1 (MQ=255)
                               cagCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGaa                                                                                                                                                                                                <  1:392196/226‑1 (MQ=255)
                                cGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTTTTTACTGTAACTTGAGTTCTGTTGGCTTGGCAGGagag                                                                                                                                                                                              >  1:691082/1‑226 (MQ=255)
                                              taCTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGGAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTTTAGGCTGGGCAGGAATGAGTCTGGTGGGCg                                                                                                                                                                                 <  1:383777/228‑1 (MQ=255)
                                               tcTGGCAGCGTTATTTCGGAAGGGTGATTAGGGCGTTCTTACGGGGAGACGCAACGAGGTTGTACGAAGGGGATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATTGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGa                                                                                                                                                                                <  2:34773/227‑1 (MQ=255)
                                               aCTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGc                                                                                                                                                                                                     <  1:499758/207‑1 (MQ=255)
                                                cTGGAAGGGTTTGTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAATTTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACGGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGat                                                                                                                                                                               <  1:638718/228‑1 (MQ=255)
                                                cTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGat                                                                                                                                                                               <  1:199583/228‑1 (MQ=255)
                                                    acgctttATTCTGTAAAGGCGATCAGGTCGGTCTTACGGGGAGACGCAACCGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGACGGGTAAAAAGTTTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGt                                                                                                                                                                           <  2:406366/223‑1 (MQ=255)
                                                    aaGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTggg                                                                                                                                                                                                      <  1:220851/201‑1 (MQ=255)
                                                         tttttttGTAAAGGGCATAAGTGTGTTTTGTCGGGGGGACGAAACCGTGTTTTTCCCAGGGCATTTTGACGGCTGCAAAAATTTTATCTGCCGCGTAAAAAGTCTTCATTTGACCGGTCTGATAC‑‑TGTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGTATGATTCTGGTGGGCGATATGTTTATTATg                                                                                                                                                                   <  2:255530/226‑1 (MQ=255)
                                                          ttttggttAAAGGGATTAGGGCTTTTTTACGGGTGGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCATAAATTTTATCGGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGTTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATg                                                                                                                                                                   <  2:250114/222‑1 (MQ=255)
                                                           ttttCCTAAAGGCGAAAAGGTCGTTCTTGCGGGGGAACTCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATg                                                                                                                                                                   <  2:431849/228‑1 (MQ=255)
                                                              tattAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGAt                                                                                                                                                                 <  1:726699/225‑1 (MQ=255)
                                                               cGTAAAGGCGATAAGGTCGTTCTTACGGGGCGACGCAACAGAGTCTTACGGAGGGGATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATctg                                                                                                                                                              <  1:480942/229‑1 (MQ=255)
                                                                                         ggggAGACGCAACAGAGTTTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGGGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTg                                                                                                                                     <  2:685931/228‑1 (MQ=255)
                                                                                             agaCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCCTTCTTTACGGTAACTTTAGTCCTGTAGGCTGGGCTGGGATGGGGCTGGTGGGGGGTATGTTTATTTTGGTCTTCTGGGTTATTTCCTTTTTTTTTTTTTTcc                                                                                                                                 >  1:431217/1‑229 (MQ=255)
                                                                                             agaCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACCGTAAATTTAGTCCCGTGAGCTGGGCAGGAATGAGTTTCGGGGGCGCTATGGTTATTATTGTCTGGTGGATTTTTCCTTTTTTATTTTTTTgcc                                                                                                                                >  1:722660/1‑229 (MQ=255)
                                                                                             agaCGCAACAGAAGCTTACGCAGGGCATTTTGACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTGTGCAATGGGGTGGACTGGCAATATTGGCGGGCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTccc                                                                                                                                <  1:625808/229‑1 (MQ=255)
                                                                                               aCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGGCATCCTTTTCGGTGAATTGAGGCCTGTAGGTCTGGGAGGGCTGAGTTTGGGGGGCGCGATGGTTATTTTGAGTTGGTGTGTTTTTCATTTTTTTTTTTTTTTCCtt                                                                                                                              >  2:522511/1‑229 (MQ=255)
                                                                                                             tACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTTAGTCCTGTAGGCTGGGCAGGAATTAGTCTGGTTGGTGATATGTTTATTATGTTCTGTCGGATTTTTCCTTTTTTTTTTTTTTCCCTTTTGCCttttttttt                                                                                                                >  1:215317/1‑227 (MQ=255)
                                                                                                               cGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTTCGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGGGATATGTTTATTTTGATCTTCTGGGTTATTCCTTTTTTTTTTTTTTCCCTTGTGCtttttgtttccc                                                                                                              >  1:221959/1‑221 (MQ=255)
                                                                                                                  aGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGGCATTTTTTTTCGGAAATTTAGTTCTGTAGGCTGGGCAGGGATGAGTCTTGTGGGGGATTTGTTTTTTTTTTTTGTTTGGGTTATTTATTTTTTTTTTTTTTTCCTTTTGGtttttttttttcccg                                                                                                           >  2:664335/1‑222 (MQ=255)
                                                                                                                    ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAACGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTTATTTTTTATTGTTACTTGAGTCCTTTGGGCTGGGCAGGGATGAGTTCTGTGGGGGATTTTTTTATTTTGGTCTTCCTGGTTATTTCTTTTTTGTTTTGTTTCCCTGAGATTTTTTTTTCCcgtgtt                                                                                                         >  1:359071/1‑225 (MQ=255)
                                                                                                                    ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCCATCTTTTCTGTGACTTGAGGCCTGTAGGCTGGGCAGGGATGGGGCTGGTTGGCGATATGTTTTTTTTGTGTTTGTGGGGTATTCCTTTTTTTTTTTTTTCCCTTTTGTTtttttttgtcgcgc                                                                                                           >  2:115767/1‑221 (MQ=255)
                                                                                                                    ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTTAGTCCTGTAGGCTGGGGAGGAATGAGTCTGGTGGGCGGTATTTTTATTTTGTTCTTCTGGGTTTTTCCTTTTTTTTTTTTTTTCCTTTGGCTTTTTTTTGCCcccgg                                                                                                          >  1:161042/1‑227 (MQ=255)
                                                                                                                    ggCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTTAGGCCTGTTGGGTGGGCAGGAAATAATCTGGGGGGCGGTATGTTTATTATTGTTTTGTGGGTTATTTCTTTTTTTATTTTTTTCCTTTGGCTTTTTGTTGGccccgg                                                                                                           >  1:697396/1‑225 (MQ=255)
                                                                                                                         ttttTCGGATGCAGAATTTTTTTCTGCCGGTTAAAAAGTCTTCGTTTTACCGGTCTGATTCTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGAtgatga                                                                                                    <  2:221959/229‑1 (MQ=255)
                                                                                                                          ttttCGGGTGCAAAAATTTTATCTGCCGGGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGAAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGAc                                                                                                   <  2:522977/229‑1 (MQ=255)
                                                                                                                          cttACGGATGCATAAATTTTAGCTGTCGCTTAAAAAGTTTTCGTTTTACCGGTCTTATACTTTTTTGCAATGGGGTGGACTTGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGACTTGGCAGGAATGAGTCTGGTGGGCTATATGTTTATTATTATCTGCTGTATTATTCCTTTTTTATTTTGTTCCCTTGTGCTTTTGTTTGCCCGCGatgatgtc                                                                                                   <  1:667230/228‑3 (MQ=255)
                                                                                                                                 aTGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTAGCGGTCTGATTCTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCt                                                                                             <  1:743399/228‑1 (MQ=255)
                                                                                                                                            tttttcggcgtgTAAAAAGTTTTTGTTTTACGGGTCTGATGCTTTTTGGCAATGGGGTGGACTGGCAATAGTGGTGGTCATTCTTTACGGTAGTTTGAGTCCTGTAGTCTGCGGCGGAATTAGTCTTGTAGGTGTTATGTTTATTATGTTCTGATTGATTATTCCTTTCTTATTTTGTTCCCTTGACCTTTTGTTTGCCTTCGATGATGACAAGCCTTGCGTTAATCgt                                                                                 <  2:163285/218‑1 (MQ=255)
                                                                                                                                                    cgcgTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTTCTGGAACTTGAGTCCTGTAGGCGGGGCGGGGATGAGGTTGGTGGGCGATATGTTTTTTTTTGTCTTCGTGGTTTTTCCTTTTTTATTTTTTTCCCTTTAGGTTTTGTTTGTCCCCGGTGATTAAACACCTTTCGTTTATTGTTTatttat                                                                         >  2:13818/1‑224 (MQ=255)
                                                                                                                                                               ttttCGTTTTACCGGTCTGATTCTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGtttt                                                              <  1:664336/227‑1 (MQ=255)
                                                                                                                                                                 ttCGTTTTACCGGTCTGATTCTTTTTTGCAAGGGGGTGGACTTGAAATAGTGGCGGTTATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGttttt                                                             <  1:678819/228‑1 (MQ=255)
                                                                                                                                                                 ttCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTTTGGGAAATTGGGTCCTGGGTGGTTGGGAGGAATGAGTATGGTGGGCGATTTGTTTTTTTTGATCTTGTGGGTTATACCCTTTTTTTTTTGTTTCCCTGTGCTTTTGTTTGTCCGCGTTGATTGCAAGTGTTTTGCTTAACGTGTAATACTACCTTTGTTTTTa                                                            >  2:339540/1‑229 (MQ=255)
                                                                                                                                                                          cctgttttatatTTTTTTGGAATGGGGGGGACTGGCAATTTGGGGGGTCATTCTTTACGGTAACTTGAGTCCGGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGAt                                                    <  2:215317/220‑1 (MQ=255)
                                                                                                                                                                                      ttttttGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGGAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTg                                       <  2:159604/229‑1 (MQ=255)
                                                                                                                                                                                        ttttGCAATGGGGTGGACTGGCAATATTGGCGGTTAATCTTTTCGGTTACTTTGGTCCTGTAGGGCGGGGAGGGATAGGTCTGGGTGGCGGTATGTTTTTTTTTATTTGCTGGGTTATTTCTTTTTTTTTTTTTTTCCTTTAGGTTTTTTTTTCCCGCGAGTATGGAAAGCCTTTCTTTTATCTATGAATTATTCCTTTGTTTTTTTCACCGCTTTTTTCatggggg                                       >  2:724709/1‑221 (MQ=255)
                                                                                                                                                                                        ttttGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAAATTGAGTTCTGTAGGCTGGGCAGGAATGAGGTTGGTGGGCGATATGTTTATTATGATCTGCTTGATTGTTCCTTTTTTATTTTGTTCCCTTGTGCTTTTGTTTGCCCGCGTTGATGTCAAGCCTTTCGTTTATCGTGTAATCCTCACTTTGTTTTTTCACCTGATTTGTTCAGTCGTGGc                                     >  1:746826/1‑229 (MQ=255)
                                                                                                                                                                                           tGCAATGGGGTGGACTGGGAATATTGGTGGTTATTCATTACGGTAACTTGAGTCCTGTAGGCTGGGCAGTAATGAGTATGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCTTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCtt                                  <  1:320508/229‑1 (MQ=255)
                                                                                                                                                                                                tGGGGTGGACTGGCAATTTTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTAtt                             <  1:522511/229‑1 (MQ=255)
                                                                                                                                                                                                  gggTGGTCTGGCAATTTTGGCGGTTGTTTTTTTGGGTAACTTGAGTCCTGTTGGCTGGGCAGGAATTGGTCTTGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATtct                           <  2:359071/229‑1 (MQ=255)
                                                                                                                                                                                                  gggTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGTTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATtct                           <  1:534385/229‑1 (MQ=255)
                                                                                                                                                                                                          tGGAAATATTGGGGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTTTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACAt                   <  1:11208/229‑1 (MQ=255)
                                                                                                                                                                                                                   tGGCGGTCATTCTTTACGGTAACTTGGTTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGTGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCg           <  1:513856/228‑1 (MQ=255)
                                                                                                                                                                                                                   tGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTTGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTGAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCg           <  1:149521/228‑1 (MQ=255)
                                                                                                                                                                                                                          taTTCTTTACGGTAATTTGAGTCCTTGAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTa   <  2:161042/228‑1 (MQ=255)
                                                                                                                                                                                                                           aTTCTTGTCGGTACCTTTAGTTCGGTAGGCTTGGCAGGAAGGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTAGTCCTTTTTTATTTTGTTCCCTTGAGCTTTTTTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCCCTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTAt  <  1:596927/229‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
GGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGGGTGGACTGGCAATATTGGCGGTCATTCTTTACGGTAACTTGAGTCCTGTAGGCTGGGCAGGAATGAGTCTGGTGGGCGATATGTTTATTATGATCTGCTGGATTATTCCTTTTTTATTTTGTTCCCTTGAGCTTTTGTTTGCCCGCGATGATGACAAGCCTTGCGTTAATCGTGTAATCATCACTTTGTTTTTACCACTGATTTGTTCAGGCGTGGCCTTCTATTCTCTCTACATCAATGTCGGAGATGTAT  >  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: