Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,676,245 C→T 5.6% R194H (CGC→CAC)  yphD ← putative sugar ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,676,2450CT 5.6% 225.1 / 2.3 71R194H (CGC→CAC) yphDputative sugar ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (36/31);  new base T (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

CCGCCAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGTCATAAACAGCCCCATTCC  >  NC_000913/2676099‑2676390
                                                                                                                                                  |                                                                                                                                                 
ccgccAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCgg                                                                                                                                                 <  2:142711/149‑1 (MQ=255)
   ccAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGTACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCt                                                                                                                                               >  5:76701/1‑148 (MQ=255)
   ccAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCt                                                                                                                                               >  5:138478/1‑148 (MQ=255)
     aacaaAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACGCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGACCGCCCGAAGGCGGTTTTGCGGCtgat                                                                                                                                            >  7:144771/1‑149 (MQ=255)
     aacaaAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCtgat                                                                                                                                            >  8:230732/1‑149 (MQ=255)
     aacaaAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCtgat                                                                                                                                            <  2:240392/149‑1 (MQ=255)
       caaAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGCCGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCtgatg                                                                                                                                           >  2:349912/1‑148 (MQ=255)
            tGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGc                                                                                                                                                   <  5:228158/135‑1 (MQ=255)
            tGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGc                                                                                                                                                   >  6:228158/1‑135 (MQ=255)
              ccGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACaaa                                                                                                                                   <  5:292841/149‑1 (MQ=255)
              ccGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACaaa                                                                                                                                   >  6:292841/1‑149 (MQ=255)
                 gTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGTGGCTGATGAACACAAACAg                                                                                                                                >  6:6421/1‑149 (MQ=255)
                 gTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAg                                                                                                                                >  7:253635/1‑149 (MQ=255)
                  tCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAgc                                                                                                                               <  1:91975/149‑1 (MQ=255)
                     ccgccgCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAgcgc                                                                                                                             >  1:53105/1‑148 (MQ=255)
                          gcTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAaacaa                                                                                                                       <  8:60846/149‑1 (MQ=255)
                          gcTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAaacaa                                                                                                                       <  8:93185/149‑1 (MQ=255)
                           cTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAaacaac                                                                                                                      <  1:70216/149‑1 (MQ=255)
                            taataaTCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAacaaca                                                                                                                     >  2:86273/1‑149 (MQ=255)
                                   cccGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGTGGCTGATGAACACAAACAGCGCAaaca                                                                                                                        >  7:172637/1‑139 (MQ=255)
                                   cccGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGTGGCTGATGAACACAAACAGCGCAaaca                                                                                                                        <  8:172637/139‑1 (MQ=255)
                                    ccGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCAt                                                                                                              <  8:253635/148‑1 (MQ=255)
                                    ccGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCAt                                                                                                              >  5:5162/1‑148 (MQ=255)
                                      gAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGAt                                                                                                           <  6:5162/149‑1 (MQ=255)
                                       aaaGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGTGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATc                                                                                                          <  5:6421/149‑1 (MQ=255)
                                        aaGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCa                                                                                                         <  7:230732/149‑1 (MQ=255)
                                           ggTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATc                                                                                                          >  6:299390/1‑145 (MQ=255)
                                           ggTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATc                                                                                                          <  5:299390/145‑1 (MQ=255)
                                            gTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAgc                                                                                                                               <  3:254853/123‑1 (MQ=255)
                                            gTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAgc                                                                                                                               >  4:254853/1‑123 (MQ=255)
                                            gTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAgcgc                                                                                                     <  4:57859/149‑1 (MQ=255)
                                             tAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGcgcg                                                                                                    <  4:244793/149‑1 (MQ=255)
                                                   tAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACCAACAGCGCAAACAACACAATCATGATCAGCGCGGATAcc                                                                                              <  5:335656/149‑1 (MQ=255)
                                                     aGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCgg                                                                                            <  5:171097/149‑1 (MQ=255)
                                                     aGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCgg                                                                                            <  6:354694/149‑1 (MQ=255)
                                                      gAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGaa                                                                                                                                         >  8:77489/1‑103 (MQ=255)
                                                      gAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGaa                                                                                                                                         <  7:77489/103‑1 (MQ=255)
                                                              aCCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGcaca                                                                                        <  5:276362/144‑1 (MQ=255)
                                                              aCCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGcaca                                                                                        >  6:276362/1‑144 (MQ=255)
                                                              aCCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCgag                                                                                   <  8:102648/149‑1 (MQ=255)
                                                                   acgaACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAAtt                                                                              >  8:198779/1‑149 (MQ=255)
                                                                    cgaACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTg                                                                             >  8:196180/1‑149 (MQ=255)
                                                                     gaACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACAGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCCCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGt                                                                            >  5:130570/1‑149 (MQ=255)
                                                                                   cacaACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTc                                                              <  3:135890/149‑1 (MQ=255)
                                                                                        cTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCa                                                         >  3:381614/1‑149 (MQ=255)
                                                                                                ccgTAGCATTACCGCCACCGGCAAAACCCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAACCAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCCGCCCCGCGttt                                                 >  7:335630/1‑149 (MQ=255)
                                                                                                ccgTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGttt                                                 >  1:357438/1‑149 (MQ=255)
                                                                                                  gTAGCATTACCGCCAACGGCAAAACCCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACACACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTtc                                               <  1:349912/149‑1 (MQ=255)
                                                                                                   tAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATc                                                                                                                >  3:56765/1‑83 (MQ=255)
                                                                                                   tAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATc                                                                                                                <  4:56765/83‑1 (MQ=255)
                                                                                                   tAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTtca                                              >  6:103587/1‑149 (MQ=255)
                                                                                                   tAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTtca                                              <  2:53105/149‑1 (MQ=255)
                                                                                                       aTTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCaa                                          >  7:88616/1‑149 (MQ=255)
                                                                                                        ttACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACAAAATCATGATCCGCGCGGATACCGGCCCACCGAGAAATTGTCCGCCCAGCCAGTCGAGCACCTCGTTTTCATCAAt                                         >  5:16968/1‑149 (MQ=255)
                                                                                                           ccGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCgg                                      >  2:162685/1‑149 (MQ=255)
                                                                                                           ccGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCgg                                      >  5:182365/1‑149 (MQ=255)
                                                                                                            cGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAATCAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTcc                                                             <  6:142708/125‑1 (MQ=255)
                                                                                                            cGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCCTGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCgg                                      >  5:378259/1‑148 (MQ=255)
                                                                                                              ccAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCGCACCGAGAAATTGTCCGCCCAGCAAGTCCAGCACCTCGTTTTCATCAATCGGCAc                                   >  7:207748/1‑149 (MQ=255)
                                                                                                              ccAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCAc                                   >  6:86707/1‑149 (MQ=255)
                                                                                                                     cAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGt                            <  1:86273/149‑1 (MQ=255)
                                                                                                                       aaaaCCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGc                          >  2:210512/1‑149 (MQ=255)
                                                                                                                            cGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCg                                       >  7:218719/1‑131 (MQ=255)
                                                                                                                            cGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCg                                       <  8:218719/131‑1 (MQ=255)
                                                                                                                             gAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGt                    <  7:196180/149‑1 (MQ=255)
                                                                                                                              aGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGt                    >  7:412145/1‑148 (MQ=255)
                                                                                                                              aGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGTc                   <  7:198779/149‑1 (MQ=255)
                                                                                                                                cgcCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGTCAt                 <  5:103587/149‑1 (MQ=255)
                                                                                                                                   ccGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGCTCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGGTTTCAGCAATCGGCACTGGCCGTGCGTTCGTCATaaa              >  5:401730/1‑149 (MQ=255)
                                                                                                                                         gCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGTCATAAACAGccc        >  3:73223/1‑149 (MQ=255)
                                                                                                                                             ttttGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGTCATAAACAGCCCCAtt    >  5:274410/1‑149 (MQ=255)
                                                                                                                                                 gCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGTCATAAACAGCCCCATTcc  <  2:351668/147‑1 (MQ=255)
                                                                                                                                                 gCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGTCATAAACAGCCCCATTcc  >  1:351668/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCGCCAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGATGCCGCACAACTGCGCCGCCGTAGCATTACCGCCAACGGCAAAAACCGAGCGCCCGAAGGCGGTTTTGCGGCTGATGAACACAAACAGCGCAAACAACACAATCATGATCAGCGCGGATACCGGCACACCGAGAAATTGTCCGCCCAGCCAGTCCAGCACCTCGTTTTCATCAATCGGCACTGGCAGTGCGTTCGTCATAAACAGCCCCATTCC  >  NC_000913/2676099‑2676390

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: