Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,083,441 C→T 5.8% P794P (CCG→CCA)  fdoG ← formate dehydrogenase‑O, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,083,4410CT 5.8% 228.6 / 2.1 69P794P (CCG→CCA) fdoGformate dehydrogenase‑O, large subunit
Reads supporting (aligned to +/- strand):  ref base C (35/30);  new base T (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGCCATCTGGTTGCCTTCC  >  NC_000913/4083303‑4083585
                                                                                                                                          |                                                                                                                                                
aTCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATcc                                                                                                                                            <  7:53653/145‑1 (MQ=255)
 tCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACg                                                                                                                                        <  8:114381/148‑1 (MQ=255)
    cgggCCAACCGACGCCCCAGGCCTTCCTGCTGCATGATAAACGGCCCGACGCAGCTACCCGGAGGCGCTGCGCGGTAGTCCGGCATCTCCCAGCCGGGCCACTTAGTGCCGGCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGt                                                                                                                                     <  3:222930/147‑1 (MQ=255)
           aCAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTg                                                                                                                              >  1:71145/1‑148 (MQ=255)
               aCGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGa                                                                                                                          <  2:71145/148‑1 (MQ=255)
                cGCCCCATGCCTTCCTGCTGAATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGAt                                                                                                                         <  2:280799/148‑1 (MQ=255)
                   cccATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTg                                                                                                                      <  1:7078/148‑1 (MQ=255)
                   cccATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTg                                                                                                                      <  7:18253/148‑1 (MQ=255)
                       tGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCCGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGa                                                                                                                  <  5:164652/148‑1 (MQ=255)
                        gCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGa                                                                                                                          >  6:228407/1‑139 (MQ=255)
                        gCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGa                                                                                                                          <  5:228407/139‑1 (MQ=255)
                         ccTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAgg                                                                                                                >  7:32874/1‑148 (MQ=255)
                              ctgctgCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGcgcg                                                                                                            <  2:324233/147‑1 (MQ=255)
                               tgctgcATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGt                                                                                                          >  7:12724/1‑148 (MQ=255)
                                gctgcATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGt                                                                                                          >  7:116496/1‑147 (MQ=255)
                                 ctgcATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGtt                                                                                                         <  8:12724/147‑1 (MQ=255)
                                    cATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTtat                                                                                                       <  5:159500/146‑1 (MQ=255)
                                    cATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTtat                                                                                                       <  6:148392/146‑1 (MQ=255)
                                      tGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTTGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTtata                                                                                                      >  5:200712/1‑145 (MQ=255)
                                      tGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTTGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTtata                                                                                                      >  3:159166/1‑145 (MQ=255)
                                         tAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAAt                                                                                                >  7:194912/1‑148 (MQ=255)
                                             cGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATgcgg                                                                                            >  3:61291/1‑148 (MQ=255)
                                              ggCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTTGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATgcggc                                                                                           <  4:159166/148‑1 (MQ=255)
                                              ggCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTTGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATgcggc                                                                                           <  6:200712/148‑1 (MQ=255)
                                                  cGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGt                                                                                                          <  7:27013/129‑1 (MQ=255)
                                                  cGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGt                                                                                                          >  8:27013/1‑129 (MQ=255)
                                                   gACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACGTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACTGAATGCGGCGGtt                                                                                       >  4:258332/1‑147 (MQ=255)
                                                        cgcTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCgg                                                                                 <  8:292472/148‑1 (MQ=255)
                                                          cTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGcc                                                                               <  7:371828/148‑1 (MQ=255)
                                                          cTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGcc                                                                               <  1:340783/148‑1 (MQ=255)
                                                             cacGGAGGCGCTGCTCTTTAGTACGGAATATCCCAGCAGGGCCACTTAGTGCCGTCCCATTTCATCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGccatg                                                                            <  2:311723/146‑1 (MQ=255)
                                                               cGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGccatgcc                                                                          >  8:115695/1‑148 (MQ=255)
                                                                ggAGGCGCTTCGAGGTATTCCTGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGTCGCTTCGGATCCCACGTGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGATTGCGGCGGTTAAGCGGccatgccc                                                                         <  6:259084/148‑1 (MQ=255)
                                                                ggAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGccatgccc                                                                         <  6:182347/148‑1 (MQ=255)
                                                                ggAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGccatgcc                                                                          >  6:1851/1‑147 (MQ=255)
                                                                       ctgcgctgTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCAtgcccatgccca                                                                  <  7:59970/148‑1 (MQ=255)
                                                                               tAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCccagcccagc                                                           <  4:327874/147‑1 (MQ=255)
                                                                                 gTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCgtg                                                        <  4:61291/148‑1 (MQ=255)
                                                                                  tCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCgtgt                                                       >  2:158708/1‑148 (MQ=255)
                                                                                   ccGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCCGCgtgt                                                       >  1:269416/1‑147 (MQ=255)
                                                                                    cGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTa                                                     <  7:115695/148‑1 (MQ=255)
                                                                                     ggAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCCCGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTAc                                                    >  1:367860/1‑148 (MQ=255)
                                                                                      gAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGccatgcccatg                                                                      <  2:100057/129‑1 (MQ=255)
                                                                                      gAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGccatgcccatg                                                                      >  1:100057/1‑129 (MQ=255)
                                                                                       aaTATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTAcc                                                   >  1:122476/1‑147 (MQ=255)
                                                                                              cAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGc                                              >  8:60352/1‑145 (MQ=255)
                                                                                              cAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAAATGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGGGTGTTACCGAGGCCAg                                           >  2:290323/1‑148 (MQ=255)
                                                                                                     tCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATc                                    >  2:251133/1‑148 (MQ=255)
                                                                                                     tCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATc                                    >  3:352107/1‑148 (MQ=255)
                                                                                                      ccACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCg                                   >  5:302081/1‑148 (MQ=255)
                                                                                                        aCTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGc                                 >  6:82049/1‑148 (MQ=255)
                                                                                                         cTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCCGCGTGTTACCGAGGCCAGACGGATCGGCg                                >  1:302392/1‑148 (MQ=255)
                                                                                                         cTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCg                                <  8:116496/148‑1 (MQ=255)
                                                                                                            aGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGt                               <  2:122476/146‑1 (MQ=255)
                                                                                                               gCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCa                          >  5:147252/1‑148 (MQ=255)
                                                                                                                ccGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCAc                         >  3:212422/1‑148 (MQ=255)
                                                                                                                        tttCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGcc                                                                               <  2:8649/86‑1 (MQ=255)
                                                                                                                        tttCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGcc                                                                               >  1:8649/1‑86 (MQ=255)
                                                                                                                        tttCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGcc                 <  5:1851/148‑1 (MQ=255)
                                                                                                                        tttCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGcc                 >  8:141248/1‑148 (MQ=255)
                                                                                                                        tttCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGc                  >  3:339755/1‑147 (MQ=255)
                                                                                                                        tttCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGc                  >  8:363906/1‑147 (MQ=255)
                                                                                                                         ttCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGCCa                >  3:300851/1‑148 (MQ=255)
                                                                                                                             gCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTACGCGGCCATGCCCATGCCCAGCCCGGAGTGTTACCGCGGCCAGACGGATCGGCGTTATCCGGGCGTGcctgctg            >  6:264399/1‑143 (MQ=255)
                                                                                                                               aaCTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGCCATCTGGt          >  5:374178/1‑148 (MQ=255)
                                                                                                                                aCTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGCCATCTGGtt         >  4:299071/1‑148 (MQ=255)
                                                                                                                                  tGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGATAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGTCATCTTGTTGc       >  6:203167/1‑148 (MQ=255)
                                                                                                                                   gccgattCGGATACCACGGGTTACCCTGCGGATCTGCGTAGTCGCGTTTATACAGAATGCGGCGGTTAAGAGGCCATGCCCAAGCCCCGCCCCGCGTCTTACCGAGGCCAGACGGATCGGCGGTAAACCGGCGTGCAATCTGGTTGcg      >  2:202157/6‑147 (MQ=255)
                                                                                                                                       cTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGCCATCTGGTTGCCTTcc  <  6:147252/148‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                
ATCGAGGGCAAACAGACGCCCCATGCCTTCCTGCTGCATGATAAACGGCCCGACGCCGCTACCCGGAGGCGCTGCGCTGTAGTCCGGAATATCCCAGCCGGTCCACTTAGTGCCGTCCCATTTCAGCAACTGACGCTTCGGATCCCACGGGTTACCCTGCGGATCTGCGGAGGCGCGGTTATACAGAATGCGGCGGTTAAGCGGCCATGCCCATGCCCAGCCCAGCGTGTTACCGAGGCCAGACGGATCGGCGTTATCACGGCGTGCCATCTGGTTGCCTTCC  >  NC_000913/4083303‑4083585

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