Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,221,839 T→C 5.1% pseudogene (2239/2574 nt) ycgH → pseudogene; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,221,8390TC 5.1% 247.6 / 2.7 78pseudogene (2239/2574 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base C (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTCCGCACCTTTATTCGT  >  NC_000913/1221693‑1221982
                                                                                                                                                  |                                                                                                                                               
ccTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAata                                                                                                                                                >  4:179960/1‑148 (MQ=255)
 cTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAata                                                                                                                                                >  5:130150/1‑147 (MQ=255)
  tGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAAtata                                                                                                                                              >  8:300697/1‑148 (MQ=255)
    ggCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACg                                                                                                                                            <  1:34408/148‑1 (MQ=255)
    ggCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTAAGGTTTTGCCCTGAATAAATCCCTTGGTCTGGAAGCGGCGAAATATACg                                                                                                                                            <  6:109989/148‑1 (MQ=255)
     gcgcAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTATGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGt                                                                                                                                           <  6:425188/148‑1 (MQ=255)
     gcgcAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGt                                                                                                                                           >  4:179933/1‑148 (MQ=255)
      cgcAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTg                                                                                                                                          <  1:281197/148‑1 (MQ=255)
                cTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAg                                                                                                                                <  4:405059/148‑1 (MQ=255)
                cTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGAGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAg                                                                                                                                >  5:191633/1‑148 (MQ=255)
                 tGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGc                                                                                                                               <  7:53006/148‑1 (MQ=255)
                         aaaCTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCcacaggca                                                                                                                       >  8:152633/1‑148 (MQ=255)
                           aCTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCcacaggcaca                                                                                                                     <  7:300697/148‑1 (MQ=255)
                                 tGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTAt                                                                                                               <  3:179933/148‑1 (MQ=255)
                                   aaaGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCcac                                                                                                                            >  4:389761/1‑133 (MQ=255)
                                   aaaGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCcac                                                                                                                            <  3:389761/133‑1 (MQ=255)
                                      ggggATGAATTAGCCTATGAATCCTATAGCGCGTTAGGTGCAAGCGTCTCGCTGGAAGCGGGTTTCGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCATCCACAGGCACAGGCTATCTGGa                                                                                                          <  6:191633/148‑1 (MQ=255)
                                             aaTTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGa                                                                                                                                        <  2:359734/111‑1 (MQ=255)
                                             aaTTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGa                                                                                                                                        >  1:359734/1‑111 (MQ=255)
                                                  gCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATc                                                                                              <  6:130150/148‑1 (MQ=255)
                                                     tATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGATGGAAGCGGGTTACGGTTTTGCACTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCAt                                                                                            >  3:252545/1‑147 (MQ=255)
                                                      aTGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTgg                                                                                                           >  7:446610/1‑131 (MQ=255)
                                                      aTGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTgg                                                                                                           <  8:446610/131‑1 (MQ=255)
                                                      aTGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATaa                                                                                          >  5:377635/1‑148 (MQ=255)
                                                      aTGAATCCTATAGCGCGAAAGGTGCAACCGTCGCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTgg                                                                                                           >  7:349844/1‑131 (MQ=255)
                                                      aTGAATCCTATAGCGCGAAAGGTGCAACCGTCGCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTgg                                                                                                           <  8:349844/131‑1 (MQ=255)
                                                       tGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATa                                                                                           <  6:377635/146‑1 (MQ=255)
                                                       tGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATa                                                                                           <  6:78179/146‑1 (MQ=255)
                                                       tGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATa                                                                                           <  7:152633/146‑1 (MQ=255)
                                                         aaTCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCtt                                                                                                                                    >  3:39392/1‑103 (MQ=255)
                                                         aaTCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCtt                                                                                                                                    <  4:39392/103‑1 (MQ=255)
                                                             cTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATg                                                                                                        >  8:247492/1‑127 (MQ=255)
                                                             cTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATg                                                                                                        <  7:247492/127‑1 (MQ=255)
                                                               ataGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTGTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCacac                                                                                 >  2:90429/1‑148 (MQ=255)
                                                                taGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGcaca                                                                                  <  6:274907/146‑1 (MQ=255)
                                                                taGCGCGAAAGGTGCAACCGTCGCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATCCGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGGCGATCATAATGCGcaca                                                                                  <  4:356102/146‑1 (MQ=255)
                                                                    gcgAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAACATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTTCATGGGCGTCGATCATAATGCGCACACGGAAg                                                                            >  6:251203/1‑148 (MQ=255)
                                                                    gcgAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAACATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAg                                                                            >  4:335319/1‑148 (MQ=255)
                                                                       aaaGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCa                                                                         <  7:56027/148‑1 (MQ=255)
                                                                         aGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAACCCAAt                                                                       >  7:237070/1‑148 (MQ=255)
                                                                             gCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTa                                                                                                                <  8:339963/103‑1 (MQ=255)
                                                                             gCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTa                                                                                                                >  7:339963/1‑103 (MQ=255)
                                                                             gCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCt                                                                   >  8:38655/1‑148 (MQ=255)
                                                                              cAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCAt                                                                                            <  5:299132/122‑1 (MQ=255)
                                                                              cAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCAt                                                                                            >  6:299132/1‑122 (MQ=255)
                                                                               aaCCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAACATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCa                                                                 <  3:335319/148‑1 (MQ=255)
                                                                                aCCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATa                                                                                           >  1:409895/1‑121 (MQ=255)
                                                                                aCCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATa                                                                                           <  2:409895/121‑1 (MQ=255)
                                                                                 ccGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAACATACGTGGATCTTCCAGCCACAGGCACAGGCTAGCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCa                                                                 <  5:251203/146‑1 (MQ=255)
                                                                                       cGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCGCACGGAAGCCAATGGCTCACGTAtt                                                           >  8:75449/1‑146 (MQ=255)
                                                                                             aaGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGa                                                   <  8:237070/148‑1 (MQ=255)
                                                                                                          ttttGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATaacaa                                       >  2:243923/1‑147 (MQ=255)
                                                                                                             tGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACAt                                    <  4:228696/147‑1 (MQ=255)
                                                                                                              gCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATgg                                                                     >  1:346509/1‑113 (MQ=255)
                                                                                                              gCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATgg                                                                     <  2:346509/113‑1 (MQ=255)
                                                                                                                 cTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCa                                                                         <  4:71684/106‑1 (MQ=255)
                                                                                                                 cTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCa                                                                         >  3:71684/1‑106 (MQ=255)
                                                                                                                  tGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGAc                              >  7:166161/1‑148 (MQ=255)
                                                                                                                   gAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGAcc                             >  1:274377/1‑148 (MQ=255)
                                                                                                                       aaaTCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATaa                                                                                          <  3:305215/83‑1 (MQ=255)
                                                                                                                       aaaTCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATaa                                                                                          >  4:305215/1‑83 (MQ=255)
                                                                                                                       aaaTCCTTTGGTCTGGAAGCGGCGAAATATACGTGGAGCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACCCGGAAGCCAATGGCTCACGTATTGAGAGTGACGCAAATCACAACATCCAGCCCCgcc                         >  4:57646/1‑146 (MQ=255)
                                                                                                                          tCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCa                                                                 >  1:187045/1‑105 (MQ=255)
                                                                                                                          tCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCa                                                                 <  2:187045/105‑1 (MQ=255)
                                                                                                                              ttGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTgag                                                        <  6:229158/110‑1 (MQ=255)
                                                                                                                              ttGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTgag                                                        >  5:229158/1‑110 (MQ=255)
                                                                                                                               tGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCtt                  <  7:38655/147‑1 (MQ=255)
                                                                                                                                ggTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAGTGACGCAAATAACAACATCCAGACCCGACTCGGCTTc                 <  3:57646/147‑1 (MQ=255)
                                                                                                                                ggTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTc                 <  8:166161/147‑1 (MQ=255)
                                                                                                                                  tCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGCCGCAAATAACAACATCCCGCCCCGACTCGGCTTCCGc              >  5:5507/1‑148 (MQ=255)
                                                                                                                                     ggAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCa                                                                         >  6:102902/1‑86 (MQ=255)
                                                                                                                                     ggAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCa                                                                         <  5:102902/86‑1 (MQ=255)
                                                                                                                                     ggAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTCCGCAcc           >  2:294514/1‑148 (MQ=255)
                                                                                                                                       aaGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGACGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTCCGCACCtt         >  6:168561/1‑148 (MQ=255)
                                                                                                                                       aaGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTCCGCACCtt         >  5:320808/1‑148 (MQ=255)
                                                                                                                                       aaGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGCATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTCCGCACAtt         >  5:358912/1‑148 (MQ=255)
                                                                                                                                        aGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTCCGCACCtt         >  8:94658/1‑147 (MQ=255)
                                                                                                                                          cggcgAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTCCGCACCTTTAt      <  1:90429/148‑1 (MQ=255)
                                                                                                                                             cgAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGCCGCAAATAACAACATCCAGCCCCGACTCGGCTTCCGCACCTTTATTCg   >  5:265406/1‑148 (MQ=255)
                                                                                                                                              gAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTCCGCACCTTTATTCGt  <  5:168561/148‑1 (MQ=255)
                                                                                                                                               aaaTATACGTGGATCTTCCAGACACAGGCACAGGCTATCTGGATGGGCTTCGATCATAATGCGAACACGGAAGACAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGCCTCGGATTCAGCACCTTTAttatt  >  1:156839/1‑144 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGATGGGCGTCGATCATAATGCGCACACGGAAGCCAATGGCTCACGTATTGAGAATGACGCAAATAACAACATCCAGACCCGACTCGGCTTCCGCACCTTTATTCGT  >  NC_000913/1221693‑1221982

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: