Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,991,679 C→T 100% S176F (TCC→TTC)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,991,6790CT99.9% 3009.4 / ‑5.4 806S176F (TCC→TTC) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (399/406);  total (399/407)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 806 total aligned reads displayed.

CGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTG  >  NC_000913/3991533‑3991824
                                                                                                                                                  |                                                                                                                                                 
      gcgcCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCGGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGc                                                                                                                                           >  4:2226866/1‑149 (MQ=255)
       cgcCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCa                                                                                                                                          <  3:4012566/149‑1 (MQ=255)
         ccAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCata                                                                                                                                        <  2:1383093/149‑1 (MQ=255)
          cAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCatat                                                                                                                                       >  4:1464279/1‑149 (MQ=255)
          cAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCatat                                                                                                                                       >  1:3994021/1‑149 (MQ=255)
             ttGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATAct                                                                                                                                    >  2:1594447/1‑149 (MQ=255)
              tGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATActg                                                                                                                                   <  4:2538694/149‑1 (MQ=255)
                cTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATActgct                                                                                                                                 <  3:964669/149‑1 (MQ=255)
                 tACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATActgctg                                                                                                                                >  4:4640878/1‑149 (MQ=255)
                  aCAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACtgctgc                                                                                                                               >  3:3536698/1‑149 (MQ=255)
                   cAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTgctgct                                                                                                                              >  3:474288/1‑149 (MQ=255)
                   cAACGTAAATGGAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTCTGGAAGTCTGCTTCTTCCTGATTGATGAAAACCGCTTCTGTCATAATGAAAGCGGCAGCCTGGGGGGCGAATATTGTGGCTTCACCCAGCATATACTgctgct                                                                                                                              <  1:4728791/149‑1 (MQ=255)
                      cGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGa                                                                                                                           <  3:4640878/149‑1 (MQ=255)
                       gTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGAc                                                                                                                          >  4:5886594/1‑149 (MQ=255)
                       gTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGAc                                                                                                                          >  1:4353800/1‑149 (MQ=255)
                       gTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGAc                                                                                                                          >  2:5349598/1‑149 (MQ=255)
                       gTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGAc                                                                                                                          <  3:2528485/149‑1 (MQ=255)
                            tGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAAtt                                                                                                                     <  2:2045261/149‑1 (MQ=255)
                              tAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAAtttt                                                                                                                   <  2:4932067/149‑1 (MQ=255)
                               aGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTa                                                                                                                  <  4:3246576/149‑1 (MQ=255)
                                gCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTAt                                                                                                                 <  2:6867207/149‑1 (MQ=255)
                                 ccTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGc                                                                                                                                           <  1:5605894/122‑1 (MQ=255)
                                     ctgGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTAc                                                                                                            >  1:3068303/1‑149 (MQ=255)
                                       ggAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTAccg                                                                                                          >  3:4344295/1‑149 (MQ=255)
                                        gAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTAccgc                                                                                                         >  2:6509899/1‑149 (MQ=255)
                                         aaaaCTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGATTTTTATCGTAccgcc                                                                                                        <  4:5268664/149‑1 (MQ=255)
                                           agcTGGGCCGCCTCGCTGAGTGTGGAAGTCAGCTTCTTCGTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACACAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGt                                                                                                      <  3:743122/147‑1 (MQ=255)
                                             cTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGACTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGc                                                                                                    <  4:1134068/149‑1 (MQ=255)
                                                ggCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGtc                                                                                                 <  2:4487742/149‑1 (MQ=255)
                                                 gccgccTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGtct                                                                                                >  1:5838603/1‑149 (MQ=255)
                                                    gccTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGAGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGc                                                                                                    >  2:2519916/1‑142 (MQ=255)
                                                     ccTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGc                                                                                             >  2:2390414/1‑148 (MQ=255)
                                                      cTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCg                                                                                           <  3:208736/149‑1 (MQ=255)
                                                       tCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCgg                                                                                          >  2:810216/1‑149 (MQ=255)
                                                       tCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCgg                                                                                          <  1:3201544/149‑1 (MQ=255)
                                                            gggTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGc                                                                                     <  3:6383856/149‑1 (MQ=255)
                                                              gtgtGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGCCAGCCGGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAAGTTTATCGTACCGCCGTGCGTCTCGCGGGTAAGCGt                                                                                   <  2:6296544/149‑1 (MQ=255)
                                                               tgtgGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTa                                                                                  >  3:680313/1‑149 (MQ=255)
                                                                 tgGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTAtt                                                                                >  4:2751815/1‑149 (MQ=255)
                                                                 tgGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTAtt                                                                                >  4:4074950/1‑149 (MQ=255)
                                                                  ggAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTc                                                                               >  2:411465/1‑149 (MQ=255)
                                                                   gAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCt                                                                              <  2:2732993/149‑1 (MQ=255)
                                                                     aGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCtgt                                                                            <  3:3113898/149‑1 (MQ=255)
                                                                     aGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCtgt                                                                            <  3:6118285/149‑1 (MQ=255)
                                                                     aGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCtgt                                                                            <  3:5957683/149‑1 (MQ=255)
                                                                      gTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCTGTAAGCGTATTCtgtg                                                                           <  3:2276326/149‑1 (MQ=255)
                                                                      gTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCtgtg                                                                           <  3:3767696/149‑1 (MQ=255)
                                                                      gTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCtgtg                                                                           <  3:6973199/149‑1 (MQ=255)
                                                                         aGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAGGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGaa                                                                        >  1:1110164/1‑149 (MQ=255)
                                                                           cttcttCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAata                                                                      >  3:5264648/1‑149 (MQ=255)
                                                                            ttcttcCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAatat                                                                     >  3:1184515/1‑149 (MQ=255)
                                                                             tcttcCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGTCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATg                                                                    <  1:7193669/149‑1 (MQ=255)
                                                                             tcttcCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATg                                                                    <  1:2525883/149‑1 (MQ=255)
                                                                             tcttcCTGATGGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATg                                                                    <  2:3467060/149‑1 (MQ=255)
                                                                              cttcCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATgg                                                                   <  3:5716680/149‑1 (MQ=255)
                                                                              cttcCTGATTGACGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATATTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATgg                                                                   >  4:6026851/1‑149 (MQ=255)
                                                                                      ttgatGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGc                                                           <  4:4418712/149‑1 (MQ=255)
                                                                                      ttgatGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGc                                                           <  1:1806610/149‑1 (MQ=255)
                                                                                      ttgatGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGc                                                           <  3:5076825/149‑1 (MQ=255)
                                                                                        gatgaAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGcga                                                         >  1:619185/1‑149 (MQ=255)
                                                                                         atgaAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGcgac                                                        <  3:3345851/149‑1 (MQ=255)
                                                                                          tgaAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGc                                                                                             <  1:2078027/111‑1 (MQ=255)
                                                                                          tgaAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGcgacg                                                       >  4:4871898/1‑149 (MQ=255)
                                                                                               aCCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACgaagaa                                                  <  1:6494860/149‑1 (MQ=255)
                                                                                               aCCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACgaagaa                                                  <  2:6146586/149‑1 (MQ=255)
                                                                                               aCCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACgaagaa                                                  <  1:7026211/149‑1 (MQ=255)
                                                                                               aCCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACgaagaa                                                  <  4:6309192/149‑1 (MQ=255)
                                                                                                ccGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATGCTGTGGAATATGGTGCCGTGCGACGaagaag                                                 >  4:68925/1‑149 (MQ=255)
                                                                                                 cGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAagaaga                                                <  1:4178642/149‑1 (MQ=255)
                                                                                                 cGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAagaaga                                                <  4:2036125/149‑1 (MQ=255)
                                                                                                  gCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAagag                                               >  1:988084/1‑149 (MQ=255)
                                                                                                    ttCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCa                                             <  1:6919188/149‑1 (MQ=255)
                                                                                                        gTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATGTTACCGCCTGGCGTCTCGCCGGTAAGCTTATTTTGTGGATTATGGTGCGGTGGGACGAAGAAGAGCATTac                                         >  4:6644562/1‑149 (MQ=255)
                                                                                                          cATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTacga                                       <  4:2896449/149‑1 (MQ=255)
                                                                                                           aTAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTacgac                                      <  1:3488108/149‑1 (MQ=255)
                                                                                                             aaTGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTAcgacga                                    <  2:3068303/149‑1 (MQ=255)
                                                                                                               tGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACt                                  >  3:1872914/1‑149 (MQ=255)
                                                                                                               tGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACt                                  <  1:2847557/149‑1 (MQ=255)
                                                                                                               tGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACt                                  >  2:4333314/1‑149 (MQ=255)
                                                                                                                gAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTa                                 >  4:174991/1‑149 (MQ=255)
                                                                                                                 aaaGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTAt                                >  1:1563272/1‑149 (MQ=255)
                                                                                                                 aaaGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTAt                                <  1:2344555/149‑1 (MQ=255)
                                                                                                                  aaGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTAtg                               >  3:1517296/1‑149 (MQ=255)
                                                                                                                  aaGCGGCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTAtg                               >  1:528132/1‑149 (MQ=255)
                                                                                                                      ggCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGtgat                           <  1:329966/149‑1 (MQ=255)
                                                                                                                       gCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGtgatg                          >  2:3688062/1‑149 (MQ=255)
                                                                                                                       gCAGCCTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGtgatg                          >  1:6615882/1‑149 (MQ=255)
                                                                                                                            cTGGGGGGCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCt                     <  4:5625453/149‑1 (MQ=255)
                                                                                                                              ggggggCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTc                                                                               >  3:2933009/1‑89 (MQ=255)
                                                                                                                              ggggggCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCttt                   <  3:2690233/149‑1 (MQ=255)
                                                                                                                                ggggCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTAc                 <  4:510931/149‑1 (MQ=255)
                                                                                                                                ggggCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTAc                 <  2:32979/149‑1 (MQ=255)
                                                                                                                                 gggCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCtgtg                                                                           <  4:3713383/90‑1 (MQ=255)
                                                                                                                                   gCGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTAcgcg              >  1:5056291/1‑149 (MQ=255)
                                                                                                                                    cGAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACgcgc             >  2:5610811/1‑149 (MQ=255)
                                                                                                                                     gAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTacgac                                      <  2:705563/123‑1 (MQ=255)
                                                                                                                                     gAAGATTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCa            <  2:5119174/149‑1 (MQ=255)
                                                                                                                                         aTTGTGGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGc        <  2:5867419/149‑1 (MQ=255)
                                                                                                                                            gtgGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTg     <  1:3688062/149‑1 (MQ=255)
                                                                                                                                            gtgGCTTCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTg     <  3:1041816/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCTTCTTCCTGATTGATGAAAACCGCTTCCGTCATAATGAAAGCGGCAGCCTGGGGGGCGAAGATTGTGGCTCCACCCAGCATATACTGCTGCTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTG  >  NC_000913/3991533‑3991824

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