Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,157,563 T→C 5.3% intergenic (+113/‑87) dkgA → / → yqhG 2,5‑diketo‑D‑gluconate reductase A/DUF3828 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,157,5630TC 5.3% 244.2 / 3.6 76intergenic (+113/‑87)dkgA/yqhG2,5‑diketo‑D‑gluconate reductase A/DUF3828 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base C (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

TCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCAC  >  NC_000913/3157417‑3157704
                                                                                                                                                  |                                                                                                                                             
tctcGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAAtata                                                                                                                                              <  1:38886/148‑1 (MQ=255)
    ggTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACCTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAg                                                                                                                                            <  8:279709/146‑1 (MQ=255)
      tCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAca                                                                                                                                        >  5:115558/1‑148 (MQ=255)
       cccGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAcat                                                                                                                                       >  5:45433/1‑148 (MQ=255)
         cGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGATCAGGAATGCCAAATATAAGGAcatca                                                                                                                                     >  5:268391/1‑148 (MQ=255)
          gATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAcatcat                                                                                                                                    >  8:320467/1‑148 (MQ=255)
           aTCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGATCAGGAATGCCAAATATAAGGAcatcatc                                                                                                                                   <  6:268391/148‑1 (MQ=255)
                  ccAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAcatcat                                                                                                                                    >  4:5578/1‑140 (MQ=255)
                  ccAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAcatcat                                                                                                                                    <  3:5578/140‑1 (MQ=255)
                      ttCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCg                                                                                                                        >  6:233875/1‑148 (MQ=255)
                      ttCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCg                                                                                                                        <  3:158162/148‑1 (MQ=255)
                                 aCATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTaa                                                                                                             <  6:45433/148‑1 (MQ=255)
                                 aCATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTa                                                                                                              <  1:153398/147‑1 (MQ=255)
                                  cATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAcatcatc                                                                                                                                   >  3:58804/1‑125 (MQ=255)
                                  cATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAcatcatc                                                                                                                                   <  4:58804/125‑1 (MQ=255)
                                   aTGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTaaaa                                                                                                           <  7:320467/148‑1 (MQ=255)
                                    tGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAg                                                                                                          <  5:233875/148‑1 (MQ=255)
                                    tGCAAATTCTCCCGGTGGCGGTAAAGTTCCGCGACCGGACGGTTCAGAAATCATTGATTCCCCTCTCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGACAAATATAATGACATCATCATGCAGAACCGGAAGCTCTTaaac                                                                                                           >  5:43074/1‑146 (MQ=255)
                                     gCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAg                                                                                                          >  7:235125/1‑147 (MQ=255)
                                      cAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTTCTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGaa                                                                                                        >  2:49425/1‑148 (MQ=255)
                                      cAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGa                                                                                                         >  2:249494/1‑147 (MQ=255)
                                      cAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGa                                                                                                         <  1:249494/147‑1 (MQ=255)
                                        aaTTCTCCCGGTGGCGGTACTGTTCCGCTACCGGACTTTTCAGATATCATTTCTTCCCCTCGCGTCCCGCCCGTTGTTACGCTTCCTGGTTCAGGAATGCCAAATATAAGGACCTCATCATGCAGAGCCGGAAGCTCTTAAAAGaaca                                                                                                      <  6:43074/148‑1 (MQ=255)
                                        aaTTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGaaca                                                                                                      <  6:185085/148‑1 (MQ=255)
                                               ccGGTGGCGGTAATGATGCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGTACAACTCATc                                                                                               <  8:111378/148‑1 (MQ=255)
                                                cGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGACGCTCTTAAAAGAACAACTCATc                                                                                               >  1:54057/1‑147 (MQ=255)
                                                cGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGGTCAGGAATGCCAAATATAAGGACATCATCATGCAGCGCCGGAAGCTCTTACACGaacaacggatc                                                                                               >  8:124674/1‑142 (MQ=255)
                                                          aaTGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGaac                                                                                                       >  1:178369/1‑129 (MQ=255)
                                                          aaTGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGaac                                                                                                       <  2:178369/129‑1 (MQ=255)
                                                            tGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATa                                                                                  <  6:115558/148‑1 (MQ=255)
                                                             gTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATa                                                                                  >  8:4633/1‑147 (MQ=255)
                                                               tCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAAc                                                                               >  4:3322/1‑148 (MQ=255)
                                                                     aCCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACg                                                                         >  1:147928/1‑148 (MQ=255)
                                                                        ggACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGa                                                                                                                                          <  1:277121/80‑1 (MQ=255)
                                                                        ggACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGa                                                                                                                                          >  2:277121/1‑80 (MQ=255)
                                                                             tttCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGt                                                                  >  1:240040/1‑147 (MQ=255)
                                                                             tttCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGt                                                                  >  2:9718/1‑147 (MQ=255)
                                                                             tttCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGt                                                                  >  4:181627/1‑147 (MQ=255)
                                                                              ttCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTg                                                                 >  7:40825/1‑147 (MQ=255)
                                                                               tCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATc                                                                                         >  6:70335/1‑122 (MQ=255)
                                                                               tCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATCCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATc                                                                                         <  5:70335/122‑1 (MQ=255)
                                                                                     aTCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATCAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAg                                                         <  4:6157/148‑1 (MQ=255)
                                                                                     aTCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAg                                                         <  7:4633/148‑1 (MQ=255)
                                                                                       cATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGAt                                                       >  2:249406/1‑148 (MQ=255)
                                                                                        aTTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGAt                                                       <  3:212201/147‑1 (MQ=255)
                                                                                         tttATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATg                                                      >  4:63848/1‑147 (MQ=255)
                                                                                            actCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGaaa                                                   <  2:27266/145‑1 (MQ=255)
                                                                                             ttCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAACAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAaat                                                 <  6:283534/148‑1 (MQ=255)
                                                                                              tCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAaat                                                 >  7:54505/1‑147 (MQ=255)
                                                                                               ccccTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAaata                                                <  8:235125/147‑1 (MQ=255)
                                                                                                cccTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAaataat                                              <  3:3322/148‑1 (MQ=255)
                                                                                                cccTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAaataat                                              <  2:147928/148‑1 (MQ=255)
                                                                                                cccTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAaataat                                              <  8:40825/148‑1 (MQ=255)
                                                                                                    cgcgTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATaa                                                                                 >  3:119190/1‑109 (MQ=255)
                                                                                                    cgcgTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATaa                                                                                 <  4:119190/109‑1 (MQ=255)
                                                                                                    cgcgTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACtt                                          >  7:268139/1‑148 (MQ=255)
                                                                                                    cgcgTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACtt                                          >  2:276617/1‑148 (MQ=255)
                                                                                                     gcgTCCCGCCCGTTGTTTCTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTc                                         <  1:49425/148‑1 (MQ=255)
                                                                                                     gcgTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACtt                                          <  1:276617/147‑1 (MQ=255)
                                                                                                          ccgcccgTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGtt                                     <  1:249406/147‑1 (MQ=255)
                                                                                                            gcccgTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGttttt                                  <  3:181627/148‑1 (MQ=255)
                                                                                                            gcccgTTGTTACTCTTCCTTGTTCAGGAATGCCAAATACAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGttttt                                  >  1:86413/1‑148 (MQ=255)
                                                                                                            gcccgTTGTTACTCTTCCTTGTTCAGGAATGCCAAATACAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGtttt                                   >  1:90827/1‑147 (MQ=255)
                                                                                                               cgTTGTTACTCTTCCTTGTTCAGGAATGCCAAATACAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTagc                               <  2:90827/148‑1 (MQ=255)
                                                                                                               cgTTGTTACTCTTCCTTGTTCAGGAATGCCAAATACAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTagc                               <  2:86413/148‑1 (MQ=255)
                                                                                                                  tgttACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTagcagc                            >  4:103402/1‑148 (MQ=255)
                                                                                                                   gttACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTagcagc                            <  3:191831/147‑1 (MQ=255)
                                                                                                                         ctTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTGATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAACCAGATGAAAATAATAGTTCTGTTTTTAGCAGCCCTGGc                      >  2:213081/1‑147 (MQ=255)
                                                                                                                         ctTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTa                                 >  4:208278/1‑136 (MQ=255)
                                                                                                                         ctTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTa                                 <  3:208278/136‑1 (MQ=255)
                                                                                                                           tCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAATCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAg                   >  8:251054/1‑148 (MQ=255)
                                                                                                                             cTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGt                  <  2:240040/147‑1 (MQ=255)
                                                                                                                                     ggAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTa         <  8:54505/148‑1 (MQ=255)
                                                                                                                                          gCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTagc                               >  7:312495/1‑121 (MQ=255)
                                                                                                                                          gCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTagc                               <  8:312495/121‑1 (MQ=255)
                                                                                                                                            cAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGca   >  8:294901/1‑147 (MQ=255)
                                                                                                                                             aaaTATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGcac  >  6:37263/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCAC  >  NC_000913/3157417‑3157704

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