Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,479 A→G 5.2% I250I (ATT→ATC fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,4790AG 5.2% 241.3 / 3.2 77I250I (ATT→ATCfdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base G (2/2);  total (36/41)
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.82e-01

ATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCAC  >  NC_000913/4080333‑4080625
                                                                                                                                                  |                                                                                                                                                  
aTGGAACTGCGGGCATAGCCTTCTTGCTCAATTTGTGCGTCCAGTACCAGAGAGTCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTaa                                                                                                                                                   <  4:44805/148‑1 (MQ=255)
  ggAACTGAGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATcg                                                                                                                                                >  4:159350/1‑149 (MQ=255)
   gAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATcgc                                                                                                                                               <  8:208111/149‑1 (MQ=255)
    aaCTGCGGGCATAGCCTTCTTGCTCCCTGCGTGCGACCAGTACCCGCGAGGCGCGGGCAGCAGCTACGGCTGCGATTTTTGGAGCTTTTTCCTGATAGAGAATTATCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATcgcg                                                                                                                                              <  2:289350/149‑1 (MQ=255)
    aaCTGCGGGCATAGCCTTCTTGCTCCATTCGTGCTTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATcgcg                                                                                                                                              <  2:330608/149‑1 (MQ=255)
     aCTGCGGTCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGATAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATcgcg                                                                                                                                              <  1:168155/148‑1 (MQ=255)
      cTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGc                                                                                                                                            <  2:149460/149‑1 (MQ=255)
            gCATAGCCTTCTTGCTGCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGAGCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTc                                                                                                                                      <  5:309167/149‑1 (MQ=255)
            gCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTc                                                                                                                                      <  6:161660/149‑1 (MQ=255)
            gCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCATGTTc                                                                                                                                      >  1:90516/1‑149 (MQ=255)
               tAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGt                                                                                                                                    >  4:59067/1‑148 (MQ=255)
                 gCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCAt                                                                                                                                 <  1:219339/149‑1 (MQ=255)
                          cTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGATTTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGAc                                                                                                                         >  6:124021/1‑148 (MQ=255)
                          cTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGAc                                                                                                                         <  4:30563/148‑1 (MQ=255)
                           tCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGc                                                                                                                                            >  8:216183/1‑128 (MQ=255)
                           tCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGc                                                                                                                                            <  7:216183/128‑1 (MQ=255)
                           tCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCGTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGAc                                                                                                                         >  2:149828/1‑147 (MQ=255)
                           tCCATTCGTGCGTCCAGGACCAGAGAGGCGATGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCTCGGCCTGTTc                                                                                                                                      <  7:129897/134‑1 (MQ=255)
                           tCCAATCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCATCGATTAATGGATCTTTTTACTGATAGAGAATTTTCAGGTAAGTCTCGCAGTCATCGCAGCTATCGTCTTTAATCGCGGCCTGTTc                                                                                                                                      >  8:129897/1‑134 (MQ=255)
                            ccATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAg                                                                                                                      >  1:226069/1‑149 (MQ=255)
                              aTTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGt                                                                                                                     >  2:231285/1‑148 (MQ=255)
                               ttCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTaa                                                                                                                   <  6:262342/149‑1 (MQ=255)
                               ttCGTGCGTCCAGTACCAGAGAGGCGAGGTCAGCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTAATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTaa                                                                                                                   <  1:45536/149‑1 (MQ=255)
                                       tCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGttt                                                                                                           >  5:43404/1‑149 (MQ=255)
                                         cAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGc                                                                                                         <  7:23125/149‑1 (MQ=255)
                                          aGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGcc                                                                                                        >  2:343689/1‑149 (MQ=255)
                                                   gagGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTc                                                                                               <  3:159350/149‑1 (MQ=255)
                                                    tggCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAAGTTTCAGGTATGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCACCGACCAGTAATGCAGTTTGCCGCTCTGCTCg                                                                                              >  1:275038/2‑149 (MQ=255)
                                                      gCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCg                                                                                              >  1:345157/1‑147 (MQ=255)
                                                            tCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTtgct                                                                                      >  4:229215/1‑149 (MQ=255)
                                                            tCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTtgct                                                                                      <  2:226069/149‑1 (MQ=255)
                                                            tCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTtgct                                                                                      <  1:343689/149‑1 (MQ=255)
                                                            tCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTtgct                                                                                      <  1:231285/149‑1 (MQ=255)
                                                             cATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGGCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTtgctg                                                                                     <  5:215473/149‑1 (MQ=255)
                                                               tCTGCTACGGCTTCGATTTTTGGATAATTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTACTCGCAGCTTTCGGATTTAATCGCGGACTGTTCGTCATCCAGAGACCAGTAATTCAGTTTGCATCTCTGATCGGAGTTGCTGCa                                                                                   >  5:317319/1‑149 (MQ=255)
                                                                  gCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCAtt                                                                                 >  5:256243/1‑148 (MQ=255)
                                                                   cTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTGATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTtg                                                                                        <  7:102528/140‑1 (MQ=255)
                                                                   cTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTGATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTtg                                                                                        >  8:102528/1‑140 (MQ=255)
                                                                   cTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCAttt                                                                                >  7:233543/1‑148 (MQ=255)
                                                                   cTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCAttt                                                                                >  8:339716/1‑148 (MQ=255)
                                                                   cTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCATGTTCGTCCTCCAGCGCCCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCGACAttt                                                                                >  6:268706/1‑148 (MQ=255)
                                                                    tACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTa                                                                              >  7:32862/1‑149 (MQ=255)
                                                                           tCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCa                                                                        >  4:101151/1‑148 (MQ=255)
                                                                                 tttGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGt                                                                  <  7:339716/148‑1 (MQ=255)
                                                                                   tGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGc                                                                <  6:256243/148‑1 (MQ=255)
                                                                                      aTCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCAtt                                                            >  5:110577/1‑149 (MQ=255)
                                                                                                 gATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTcaca                                                  >  8:274361/1‑148 (MQ=255)
                                                                                                  aTAGAGTATTGTCAGGAACGTGTCGCAGGCATCGCAGCTTTCGGCTTTAATGGCGGCCTGTTCGGCATCCAGCGACCAGTAATGCAGTATGCCGCACTGCTCGCAGTTGCTGCATTATACGCGCACTACGTGCCATTCAGTTTcaccc                                                 <  5:65325/148‑3 (MQ=255)
                                                                                                  aTAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTcacaca                                                <  4:163644/149‑1 (MQ=255)
                                                                                                  aTAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTGCGGCTTTACTCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTcacaca                                                <  2:275038/149‑1 (MQ=255)
                                                                                                   tAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTGATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAg                                               >  8:43646/1‑149 (MQ=255)
                                                                                                   tAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGAGACCAGTAATGCAGTTTGCAGCTCTGCTCGCAGTTGCTGCATTTTACGAGCACTGCGTGCCATTCAGTTTCAGACAg                                               >  4:26956/1‑149 (MQ=255)
                                                                                                         aTTTTCAGGTTAGTGTCGACGTCATCGCAGCTTTCGGCTTTACTCGCGGCCTGTTCGTCATCCCGCGACCAGTACTGCAGTGTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTtgcc                                         <  6:305168/149‑2 (MQ=255)
                                                                                                         aTTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTGCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGAAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTtgca                                         <  3:242758/149‑1 (MQ=255)
                                                                                                         aTTTTCAGGTAAGGGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCACTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTtgca                                         <  4:310026/149‑1 (MQ=255)
                                                                                                                ggTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAgg                                  >  7:179604/1‑149 (MQ=255)
                                                                                                                 gTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTTCAgg                                  >  2:73763/1‑148 (MQ=255)
                                                                                                                  tAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGt                                 >  6:302341/1‑148 (MQ=255)
                                                                                                                       gtCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACg                             >  6:173762/1‑147 (MQ=255)
                                                                                                                        tCGCAGGCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAg                          <  8:187760/149‑1 (MQ=255)
                                                                                                                            aGTCATCGCAGCTTGCGGCTTTAATCGCGGCCTGGTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGcccc                      <  7:175540/149‑1 (MQ=255)
                                                                                                                             gTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCt                     >  8:81154/1‑149 (MQ=255)
                                                                                                                              tCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCt                     >  6:17498/1‑148 (MQ=255)
                                                                                                                                aTCGCAGCTTTCGGCTTTACTCGCGGCCTGTTCGTCATCCCCCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGGTGCTGCATTTTACGCGCACTACGTGCCTTTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGa                   <  8:272648/148‑1 (MQ=255)
                                                                                                                                aTCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGa                   <  3:59067/148‑1 (MQ=255)
                                                                                                                                   gCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGATCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGAAGTTAACGTAGCCCCTGAgtgg               >  8:106243/1‑149 (MQ=255)
                                                                                                                                    cAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAgtggt              <  3:101151/149‑1 (MQ=255)
                                                                                                                                    cAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAgtggt              <  2:148633/149‑1 (MQ=255)
                                                                                                                                          tCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCaa         <  2:345157/148‑1 (MQ=255)
                                                                                                                                          tCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCaa         <  5:124021/148‑1 (MQ=255)
                                                                                                                                           cGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAAt        <  7:274361/148‑1 (MQ=255)
                                                                                                                                           cGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAAt        <  6:53308/148‑1 (MQ=255)
                                                                                                                                              cTTTGATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGc    <  7:43646/149‑1 (MQ=255)
                                                                                                                                              cTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCAGTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTCGCCCCTGAGTGGTGCCAATTTg     >  5:216731/1‑148 (MQ=255)
                                                                                                                                               tttAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa   >  4:277737/1‑149 (MQ=255)
                                                                                                                                               tttAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGAAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa   >  6:294115/1‑149 (MQ=255)
                                                                                                                                                atAATCGCGGCCTGGGCCTCATCCAGCGCCCAGTACGGCAGGTTGCCGCTCTGCTCGCAGGGGCGGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCAc  <  2:224990/148‑1 (MQ=255)
                                                                                                                                                 tAATCGCGGCCGGTGCGCCATCCAGCGACCAGGAACGCAGTGTGCCGCTCTGCTCGCAGTTGCTCCATTTAACGCGCACAACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCAc  <  6:138499/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATGGAACTGCGGGCATAGCCTTCTTGCTCCATTCGTGCGTCCAGTACCAGAGAGGCGAGGTCATCTGCTACGGCTTCGATTTTTGGATCTTTTTCCTGATAGAGAATTTTCAGGTAAGTGTCGCAGTCATCGCAGCTTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCAC  >  NC_000913/4080333‑4080625

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