Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,621 T→G 6.5% Q203P (CAA→CCA)  fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,6210TG 6.5% 323.5 / 12.7 92Q203P (CAA→CCA) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base T (40/46);  new base G (3/3);  total (43/49)
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.28e-01

GCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCG  >  NC_000913/4080474‑4080767
                                                                                                                                                   |                                                                                                                                                  
gCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTg                                                                                                                                                   >  2:614276/1‑149 (MQ=255)
 cTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGc                                                                                                                                                  <  4:112957/149‑1 (MQ=255)
  tttAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa                                                                                                                                                 >  1:507807/1‑149 (MQ=255)
  tttAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa                                                                                                                                                 >  3:530928/1‑149 (MQ=255)
        cgcgGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGc                                                                                                                                           <  4:62158/149‑1 (MQ=255)
          cgGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCactcat                                                                                                                                         <  4:310490/149‑6 (MQ=255)
          cgGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCactcag                                                                                                                                         <  4:74796/149‑6 (MQ=255)
             ccTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGac                                                                                                                                      >  1:449474/1‑149 (MQ=255)
              cTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGaca                                                                                                                                     >  2:194655/1‑149 (MQ=255)
                gTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACAcc                                                                                                                                   >  4:472850/1‑149 (MQ=255)
                 ttCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCg                                                                                                                                  <  2:74922/149‑1 (MQ=255)
                     tCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCAt                                                                                                                              <  2:387402/149‑1 (MQ=255)
                      cATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATa                                                                                                                             <  4:228021/149‑1 (MQ=255)
                         ccAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGaa                                                                                                                          <  1:504359/149‑1 (MQ=255)
                          cAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAAc                                                                                                                         >  1:463801/1‑149 (MQ=255)
                          cAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAAc                                                                                                                         >  3:550232/1‑149 (MQ=255)
                          cAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAAc                                                                                                                         <  4:15746/149‑1 (MQ=255)
                            gCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCTTTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACca                                                                                                                       <  3:516432/149‑1 (MQ=255)
                            gCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACca                                                                                                                       >  1:341664/1‑149 (MQ=255)
                                   gTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCac                                                                                                                                             >  2:581469/1‑119 (MQ=255)
                                   gTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCac                                                                                                                                             <  1:581469/120‑2 (MQ=255)
                                    tAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                               <  2:463801/149‑1 (MQ=255)
                                    tAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                               <  2:350113/149‑1 (MQ=255)
                                        gCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCaa                                                                                                           >  2:61557/1‑149 (MQ=255)
                                            tttGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATAtt                                                                                                       <  2:294218/149‑1 (MQ=255)
                                              tGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCGCAAACCGGGCAATATTGa                                                                                                     <  2:507807/149‑1 (MQ=255)
                                               gCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGAc                                                                                                    <  1:183978/149‑1 (MQ=255)
                                                 cGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGt                                                                                                  >  3:438184/1‑149 (MQ=255)
                                                 cGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGt                                                                                                  >  2:539242/1‑149 (MQ=255)
                                                   ctctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCAGCATGCTGGAAACCGGCATAGAACCACAAACCGGGCAATATTGACGTTg                                                                                                >  4:47376/1‑149 (MQ=255)
                                                   ctctGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTg                                                                                                <  4:530928/149‑1 (MQ=255)
                                                      tGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCTGCATAGAACCACAACCCGAGCAATATTGACGTTGCTc                                                                                             <  4:79922/149‑1 (MQ=255)
                                                         tCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGcc                                                                                          >  2:487837/1‑149 (MQ=255)
                                                                tgctgcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGGTGCTCGCCGTATTCa                                                                                   >  1:543542/1‑149 (MQ=255)
                                                                 gctgcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAg                                                                                  <  1:22325/149‑1 (MQ=255)
                                                                   tgcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAgcg                                                                                >  2:362212/1‑149 (MQ=255)
                                                                      aTTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCgg                                                                             <  1:128493/149‑1 (MQ=255)
                                                                       ttttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCggg                                                                            >  2:41656/1‑149 (MQ=255)
                                                                            cgcgCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTg                                                                       <  2:327969/149‑1 (MQ=255)
                                                                              cgcACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGcc                                                                     <  3:438797/149‑1 (MQ=255)
                                                                                cACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgg                                                                   <  2:560282/149‑1 (MQ=255)
                                                                                 aCTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCggg                                                                  >  1:203511/1‑149 (MQ=255)
                                                                                 aCTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCggg                                                                  <  3:499676/149‑1 (MQ=255)
                                                                                 aCTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCggg                                                                  <  3:472850/149‑1 (MQ=255)
                                                                                 aCTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCggg                                                                  >  4:340691/1‑149 (MQ=255)
                                                                                   tACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGa                                                                <  3:619937/149‑1 (MQ=255)
                                                                                    aCGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGAt                                                               <  2:457422/149‑1 (MQ=255)
                                                                                      gTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCa                                                             >  3:485679/1‑149 (MQ=255)
                                                                                      gTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCa                                                             <  2:266082/149‑1 (MQ=255)
                                                                                        gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGa                                                           <  4:438184/149‑1 (MQ=255)
                                                                                        gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGa                                                           <  2:60394/149‑1 (MQ=255)
                                                                                          cATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAtt                                                         <  4:263139/149‑1 (MQ=255)
                                                                                          cATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAtt                                                         >  1:580406/1‑149 (MQ=255)
                                                                                          cATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAtt                                                         <  2:27082/149‑1 (MQ=255)
                                                                                           aTTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTg                                                        <  1:354675/149‑1 (MQ=255)
                                                                                            ttCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGtt                                                                                                 >  4:618726/1‑107 (MQ=255)
                                                                                            ttCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGtt                                                                                                 <  3:618726/107‑1 (MQ=255)
                                                                                            ttCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTgg                                                       >  2:105656/1‑149 (MQ=255)
                                                                                                 tttCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATc                                                  >  3:112936/1‑149 (MQ=255)
                                                                                                     acacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTgcg                                              >  3:6172/1‑149 (MQ=255)
                                                                                                       acaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGcc                                            >  4:163216/1‑149 (MQ=255)
                                                                                                       acaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGcc                                            >  1:468990/1‑149 (MQ=255)
                                                                                                         aGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCa                                          <  1:211135/149‑1 (MQ=255)
                                                                                                           gTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCg                                                                                                                                  >  4:490956/1‑59 (MQ=255)
                                                                                                           gTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCg                                                                                                                                  <  3:490956/59‑1 (MQ=255)
                                                                                                           gTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGt                                        >  2:9287/1‑149 (MQ=255)
                                                                                                              gcagtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTaga                                     <  1:61557/149‑1 (MQ=255)
                                                                                                              gcagtgcagGTAACGTAGCCCCTGAGTGGTACCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTaga                                     <  2:62847/149‑1 (MQ=255)
                                                                                                                agtgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTGTTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGc                                   >  4:18585/1‑149 (MQ=255)
                                                                                                                  tgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGa                                 <  1:365322/149‑1 (MQ=255)
                                                                                                                  tgcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGa                                 <  4:199999/149‑1 (MQ=255)
                                                                                                                   gcagGTAACGTAGCCCCTGTGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGAc                                >  3:493905/1‑149 (MQ=255)
                                                                                                                   gcagGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGAc                                <  3:254186/149‑1 (MQ=255)
                                                                                                                     agGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAg                              >  3:65040/1‑149 (MQ=255)
                                                                                                                       gTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTg                            >  4:242891/1‑149 (MQ=255)
                                                                                                                            gTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGcc                       <  4:440643/149‑1 (MQ=255)
                                                                                                                                ccccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGa                   >  1:87714/1‑149 (MQ=255)
                                                                                                                                  ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATa                 >  1:554085/1‑149 (MQ=255)
                                                                                                                                  ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATa                 >  2:146265/1‑149 (MQ=255)
                                                                                                                                   cagcatggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATaa                >  1:68650/6‑149 (MQ=255)
                                                                                                                                   cTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATaa                >  1:139783/1‑149 (MQ=255)
                                                                                                                                     gAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAAc              <  4:112936/149‑1 (MQ=255)
                                                                                                                                     gAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAAc              <  2:177055/149‑1 (MQ=255)
                                                                                                                                      catggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACg             >  3:520606/3‑149 (MQ=255)
                                                                                                                                       gtggtgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACgg            >  1:400517/1‑149 (MQ=255)
                                                                                                                                           tgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCt        <  1:487837/149‑1 (MQ=255)
                                                                                                                                           tgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCt        <  1:194655/149‑1 (MQ=255)
                                                                                                                                            gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCtt       <  1:41656/149‑1 (MQ=255)
                                                                                                                                              cAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTa     >  1:609191/1‑149 (MQ=255)
                                                                                                                                               aaTTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGAGCTTTAt    <  2:341664/149‑1 (MQ=255)
                                                                                                                                                aTTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATc   <  2:468990/149‑1 (MQ=255)
                                                                                                                                                aTTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGAAGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATc   >  4:81441/1‑149 (MQ=255)
                                                                                                                                                 tttGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATGTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCg  >  3:512331/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTATCG  >  NC_000913/4080474‑4080767

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