Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,616 C→A 6.6% G205C (GGC→TGC)  fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,6160CA 6.6% 320.1 / 11.2 91G205C (GGC→TGC) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base C (40/45);  new base A (3/3);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

TTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTA  >  NC_000913/4080469‑4080764
                                                                                                                                                   |                                                                                                                                                    
tttCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCa                                                                                                                                                     >  3:199999/1‑149 (MQ=255)
 ttCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCaa                                                                                                                                                    >  2:245796/1‑149 (MQ=255)
   cGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTCACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAAtt                                                                                                                                                  <  4:258100/149‑1 (MQ=255)
     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                                                                                                                                          <  1:581469/120‑2 (MQ=255)
                                        gTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCac                                                                                                                                          >  2:581469/1‑119 (MQ=255)
                                         tAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                            <  2:350113/149‑1 (MQ=255)
                                         tAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                            <  2:463801/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                                                               >  4:340691/1‑149 (MQ=255)
                                                                                      aCTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCggg                                                               <  3:472850/149‑1 (MQ=255)
                                                                                      aCTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCggg                                                               >  1:203511/1‑149 (MQ=255)
                                                                                      aCTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCggg                                                               <  3:499676/149‑1 (MQ=255)
                                                                                        tACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGa                                                             <  3:619937/149‑1 (MQ=255)
                                                                                         aCGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGAt                                                            <  2:457422/149‑1 (MQ=255)
                                                                                           gTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCa                                                          <  2:266082/149‑1 (MQ=255)
                                                                                           gTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCa                                                          >  3:485679/1‑149 (MQ=255)
                                                                                             gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGa                                                        <  2:60394/149‑1 (MQ=255)
                                                                                             gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGa                                                        <  4:438184/149‑1 (MQ=255)
                                                                                               cATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAtt                                                      <  2:27082/149‑1 (MQ=255)
                                                                                               cATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAtt                                                      <  4:263139/149‑1 (MQ=255)
                                                                                               cATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAtt                                                      >  1:580406/1‑149 (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                                                                                                                               <  3:490956/59‑1 (MQ=255)
                                                                                                                gTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCg                                                                                                                               >  4:490956/1‑59 (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           <  2:177055/149‑1 (MQ=255)
                                                                                                                                          gAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAAc           <  4:112936/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)
                                                                                                                                                   |                                                                                                                                                    
TTTCGGCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTA  >  NC_000913/4080469‑4080764

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: