Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,806,935 A→T 5.2% K189I (AAA→ATA)  pfkB → 6‑phosphofructokinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,806,9350AT 5.2% 263.2 / 2.8 83K189I (AAA→ATA) pfkB6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base A (32/45);  new base T (2/2);  total (35/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 = 9.39e-01

CCGCCAGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGT  >  NC_000913/1806790‑1807083
                                                                                                                                                 |                                                                                                                                                    
ccgccAGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCaaaaa                                                                                                                                                     <  8:386281/147‑1 (MQ=255)
 cgccAGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGaa                                                                                                                                                  <  3:301317/149‑1 (MQ=255)
    cAGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCaaaaa                                                                                                                                                     >  2:387904/1‑143 (MQ=255)
    cAGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCaaaaa                                                                                                                                                     <  1:387904/143‑1 (MQ=255)
    cAGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACt                                                                                                                                                >  3:196236/1‑148 (MQ=255)
     aGGTGTGAAGCTGTAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGATGCATCGTCGACAGTTCTGGCGAAGCGGTAAGTGCAGCACTGGCAATTGGTAACATCTAGTTGGTTAAGCCTAACCAAAAAGAACTCa                                                                                                                                              <  8:156447/149‑1 (MQ=255)
     aGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACt                                                                                                                                                >  5:198114/1‑147 (MQ=255)
     aGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCa                                                                                                                                              <  2:48502/149‑1 (MQ=255)
     aGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCa                                                                                                                                              <  8:301570/149‑1 (MQ=255)
     aGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCa                                                                                                                                              <  1:95218/149‑1 (MQ=255)
       gtgtGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGt                                                                                                                                            <  1:251981/149‑1 (MQ=255)
           gAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTgcgc                                                                                                                                        >  7:356556/1‑149 (MQ=255)
            aaGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTgcgc                                                                                                                                        <  2:200738/148‑1 (MQ=255)
            aaGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTgcgc                                                                                                                                        <  5:101811/148‑1 (MQ=255)
            aaGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCAGCTGCATCGTCGACACTTCTGGCGAAGCGTTAACTGCAGCACTGGCAATTGGTAACAACGAGTTGGTTAAGCCTAACCAAAAAGAACTAAGTGCGCt                                                                                                                                       >  2:162764/1‑149 (MQ=255)
             aGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTg                                                                                                                                      <  6:186943/149‑1 (MQ=255)
             aGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTg                                                                                                                                      <  4:80222/149‑1 (MQ=255)
             aGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTg                                                                                                                                      >  3:296258/1‑149 (MQ=255)
              gCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAATCAAGTGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTATGTGCAGCACTGGCTATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTgg                                                                                                                                     <  1:4080/149‑1 (MQ=255)
               cTGGAAAAATTCACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCAGCGTCGACAGTTCAGGCGAAGCGTTAAGTGCATCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGt                                                                                                                                    <  5:282681/149‑1 (MQ=255)
                  gataaaTTAAGCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCTATTGGTAACATGGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGt                                                                                                                                    <  1:162764/143‑1 (MQ=255)
                  gAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGa                                                                                                                                  <  5:370279/148‑1 (MQ=255)
                  gAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCATCACCATCGACGACACTTCTGGCGAAGCGTAAAGTGCAGCACTGTCAATAGTTAACATCGAGTTTTTTAAGCCTAACAAAATACAACTCAGTGAGCTgcaga                                                                                                                                  >  2:35432/1‑144 (MQ=255)
                   aaaaaTTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAAt                                                                                                                                <  5:220704/149‑1 (MQ=255)
                      aaTTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAATAGAACTCAGTGCGCTGGTGAATcgc                                                                                                                             >  7:287411/1‑149 (MQ=255)
                            cccAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGGTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTc                                                                                                                       <  1:133238/149‑1 (MQ=255)
                               aaCTGATTTCCGCTGCGCAAATACAAGGGATCCGCTGCATGGTCGTCAGTTCTGGCTTAGTGTTATGTGGAGCACTTGCAATTGGTATCATCGAGTTTGTTAAGGCTAACCAAAAAGAACTCAGTGCGCTGGTGATTGGCGAACTCTcc                                                                                                                    <  6:223288/149‑1 (MQ=255)
                               aaCTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTACGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCAcc                                                                                                                    >  1:381464/1‑149 (MQ=255)
                               aaCTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTGAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGCACTCAcc                                                                                                                    >  2:296238/1‑149 (MQ=255)
                                 cTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCa                                                                                                                  <  6:145757/149‑1 (MQ=255)
                                      ttCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAGAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCg                                                                                                              >  4:13007/1‑148 (MQ=255)
                                         cGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGAc                                                                                                           <  7:331450/148‑1 (MQ=255)
                                              cgcAAACACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGAta                                                                                                       <  3:78620/147‑2 (MQ=255)
                                              cgcAAAAACAAGGGATCCGCTGCATCGTCGACAGTGCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGAta                                                                                                       <  5:127352/147‑2 (MQ=255)
                                              agaaaaaaCAAGGGATCCTCTGCATCGTCGACAGTTCTCGCGAAGCGTTCAGAGCAGAACTGGCAATTGGTACCATCGAGTTGGTTAAGCCTAACCCAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGt                                                                                                      <  8:274803/145‑1 (MQ=255)
                                               gcAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATcc                                                                                                    >  4:188352/1‑149 (MQ=255)
                                               gcAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCAGGACGATATcc                                                                                                    >  4:246005/1‑149 (MQ=255)
                                                 aaaaaCAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGc                                                                                                  <  4:296258/149‑1 (MQ=255)
                                                   aaaCAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCa                                                                                                 >  8:253964/1‑148 (MQ=255)
                                                    aaCAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCa                                                                                                 >  4:267899/1‑147 (MQ=255)
                                                     aCAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAg                                                                                              >  3:126877/1‑149 (MQ=255)
                                                      cAAGGGATCCGCAGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGc                                                                                             >  4:353060/1‑149 (MQ=255)
                                                       aaGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGcc                                                                                            >  7:143605/1‑149 (MQ=255)
                                                           gATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCgcgc                                                                                        <  1:234624/149‑1 (MQ=255)
                                                                  tGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGa                                                                                                         <  7:161258/125‑1 (MQ=255)
                                                                  tGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGa                                                                                                         >  8:161258/1‑125 (MQ=255)
                                                                    cATCGTCGACAGTTCTCGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGCAAATCg                                                                               >  7:2842/1‑149 (MQ=255)
                                                                     aTCGTCGCCAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGt                                                                              <  8:381589/149‑1 (MQ=255)
                                                                      tcgtcgACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAGAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGtt                                                                             <  3:13007/149‑1 (MQ=255)
                                                                      tcgtcgACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGtt                                                                             <  1:202538/149‑1 (MQ=255)
                                                                            aCAGTTCTGGCTAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGc                                                                       <  7:28599/149‑1 (MQ=255)
                                                                             cAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCg                                                                      <  4:196236/149‑1 (MQ=255)
                                                                                ttCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGc                                                                    >  3:364742/1‑148 (MQ=255)
                                                                                  cTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAg                                                                 >  2:52510/1‑149 (MQ=255)
                                                                                     gCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGc                                                               >  1:205565/1‑148 (MQ=255)
                                                                                     gCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGc                                                               <  8:356556/148‑1 (MQ=255)
                                                                                         aGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAATAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAAc                                                          <  1:272900/149‑1 (MQ=255)
                                                                                         aGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAAc                                                          <  5:239908/149‑1 (MQ=255)
                                                                                          gCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACg                                                         >  7:320444/1‑149 (MQ=255)
                                                                                            gTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAATAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACggg                                                       <  5:95377/149‑1 (MQ=255)
                                                                                             ttAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGt                                                      >  8:13867/1‑149 (MQ=255)
                                                                                                 gTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAATTCGTTAATAGCGGCAAGGCCAAACGGGTTGTc                                                  <  8:2842/149‑1 (MQ=255)
                                                                                                 gTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTc                                                  <  5:293690/149‑1 (MQ=255)
                                                                                                     agcaCTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGtt                                                                             <  1:173009/118‑1 (MQ=255)
                                                                                                     agcaCTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGtt                                                                             >  2:173009/1‑118 (MQ=255)
                                                                                                      gcaCTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAATCCGCGCAGGACATCGTTAATAGCGGCAAGGCCAAACGGGGAGTCGTTTc                                             >  5:307820/1‑149 (MQ=255)
                                                                                                         cTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCt                                          <  1:143122/149‑1 (MQ=255)
                                                                                                          tGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTg                                         >  8:138821/1‑149 (MQ=255)
                                                                                                              aaTTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTc                                     <  4:376838/149‑1 (MQ=255)
                                                                                                              aaTTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTc                                     <  1:52510/149‑1 (MQ=255)
                                                                                                               aTTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTcc                                    <  3:353060/149‑1 (MQ=255)
                                                                                                                  ggTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCcac                                  <  5:232653/148‑1 (MQ=255)
                                                                                                                   gTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATTTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCcac                                  >  3:275530/1‑147 (MQ=255)
                                                                                                                    tAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAg                               <  7:138821/149‑1 (MQ=255)
                                                                                                                    tAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAg                               >  2:342585/1‑149 (MQ=255)
                                                                                                                         tCGAGTTGGTTAAGCCTAACCAAAAAGACCTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAACGCCAAACGGGTTCTCGTTTCCCTGGCTCCACAAGGAgcg                          >  5:293582/1‑149 (MQ=255)
                                                                                                                           gAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCt                        <  3:267899/149‑1 (MQ=255)
                                                                                                                                 gTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGgtgt                  <  4:126877/149‑1 (MQ=255)
                                                                                                                                  ttAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGtt                 >  2:81634/1‑149 (MQ=255)
                                                                                                                                     aGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGa               <  3:188352/148‑1 (MQ=255)
                                                                                                                                      gCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATa             <  3:246005/149‑1 (MQ=255)
                                                                                                                                         tAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCAAAAGGAGCGCTGGGTGTTGATAGt           >  5:193840/1‑148 (MQ=255)
                                                                                                                                         tAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGGCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGt           >  3:316667/1‑148 (MQ=255)
                                                                                                                                             cAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCCCCCAGCCGGACGAAATCCGCAAAGCCGCGCAGGAAATCGATAATAGCGGCAAGGCCAAACGCATTGTCGTTTCCCAGCGTCCACAAGGAGCGCTCGGTGTTGATAGTGaaa       >  3:260376/1‑148 (MQ=255)
                                                                                                                                             cAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGaaaa      >  7:194976/1‑149 (MQ=255)
                                                                                                                                             cAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCAAGCCGGAAGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAATCGGGTTGTCGTTTCCCTGGGTCCACAAGGTGCGCTGGGTGTTGATAGTGaaaa      <  6:293582/149‑1 (MQ=255)
                                                                                                                                              aaaaaGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAAt                                                                          <  4:50485/80‑1 (MQ=255)
                                                                                                                                              aaaaaGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAAt                                                                          >  3:50485/1‑80 (MQ=255)
                                                                                                                                                aaaGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGAGATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAACCGGGTGGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGTGGTTGATAGTGAAAACTg   >  2:247168/1‑149 (MQ=255)
                                                                                                                                                 aaGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGt  <  1:342585/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCGCCAGGTGTGAAGCTGGAAAAATTAACCCAACTGATTTCCGCTGCGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGT  >  NC_000913/1806790‑1807083

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