Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,806,982 G→A 82.5% V205I (GTC→ATC)  pfkB → 6‑phosphofructokinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,806,9820GA82.5% 166.6 / 28.2 75V205I (GTC→ATC) pfkB6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base G (3/10);  new base A (32/29);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 7.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTAC  >  NC_000913/1806836‑1807121
                                                                                                                                                  |                                                                                                                                           
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)
                      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                                                                                                     <  8:356556/148‑1 (MQ=255)
                                       gCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGc                                                                                                     >  1:205565/1‑148 (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                                                                           <  1:52510/149‑1 (MQ=255)
                                                                aaTTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTc                                                                           <  4:376838/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)
                                                                                                    aGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTa                                       >  3:342703/1‑149 (MQ=255)
                                                                                                     gAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTa                                       <  8:287411/148‑1 (MQ=255)
                                                                                                                gcTGGTGAATCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGc                           >  6:138258/1‑149 (MQ=255)
                                                                                                                         tCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTc                                                                           >  2:8664/1‑92 (MQ=255)
                                                                                                                         tCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTc                                                                           <  1:8664/92‑1 (MQ=255)
                                                                                                                         tCGCGAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTGTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCgg                  >  8:298393/1‑149 (MQ=255)
                                                                                                                           gcgAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGt                 >  7:39727/1‑148 (MQ=255)
                                                                                                                           gcgAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTg                <  4:342703/149‑1 (MQ=255)
                                                                                                                            cgAACTCACCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGa               <  7:253964/149‑1 (MQ=255)
                                                                                                                              aaCTCCCCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGCGCGCTGGGTGTTGATAGGGAAAACTGAATTCAGGTGGTGCCACCACCGGGGaa              >  6:355648/1‑148 (MQ=255)
                                                                                                                              aaCTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGaaa             <  8:194976/149‑1 (MQ=255)
                                                                                                                               aCTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCAAAAGGAGCGCTGGGTGGTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGaaaa            <  6:193840/149‑1 (MQ=255)
                                                                                                                                   aCCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACTACCGGTGAAAAGc          <  8:320444/147‑1 (MQ=255)
                                                                                                                                   aCCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAACCGGGGGGTCGTTTCCCTGTGTACACACGGAGGGCTGGCGGGAGAGAGTGAAAACGGTATTCAGGTGGTGCCCGCACCGGTGAAAAGc          >  4:30029/1‑147 (MQ=255)
                                                                                                                                   aCCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGc          >  1:260056/1‑147 (MQ=255)
                                                                                                                                   aCCCAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGGTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGcc         <  1:296238/148‑1 (MQ=255)
                                                                                                                                     ccAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCa                                   <  6:297130/120‑1 (MQ=255)
                                                                                                                                     ccAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCa                                   >  5:297130/1‑120 (MQ=255)
                                                                                                                                     ccAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCaga      <  5:138258/149‑1 (MQ=255)
                                                                                                                                     ccAGCCGGACGATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGCCCAAACGGGTTGTCGATTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCCCCGGTGAAAAGCCa        >  7:96112/1‑147 (MQ=255)
                                                                                                                                         acGGACTATATCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAAAGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGTTGTTGATAGTGAAAACTGTATTCATGTGGTGCCACCACCGGTGAAAAGCCAGAGTAc  <  5:210416/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
CGCAAAAACAAGGGATCCGCTGCATCGTCGACAGTTCTGGCGAAGCGTTAAGTGCAGCACTGGCAATTGGTAACATCGAGTTGGTTAAGCCTAACCAAAAAGAACTCAGTGCGCTGGTGAATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTAC  >  NC_000913/1806836‑1807121

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: ‑