Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,807,102 C→T 33.3% P245S (CCG→TCG)  pfkB → 6‑phosphofructokinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,807,1020CT33.3% 79.1 / 55.7 66P245S (CCG→TCG) pfkB6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base T (9/13);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

ATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACA  >  NC_000913/1806956‑1807245
                                                                                                                                                  |                                                                                                                                               
aTCGCGAACTCACCCGGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCg                                                                                                                                               <  6:176203/149‑1 (MQ=255)
 tCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCCGGTGGTGCCACCACCgg                                                                                                                                              >  5:314290/1‑149 (MQ=255)
     gAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGa                                                                                                                                           >  8:234541/1‑148 (MQ=255)
       aCTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGaa                                                                                                                                          >  3:169118/1‑147 (MQ=255)
           ccccAGCCGGACGATGTCCGCAAAGACGCGAAGGAAATCGTTAATAGCGGCAAGGCCAACCTGGTTGGCGTGTCGCTGGGTCCCCAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCCCCACCGGTGAAACGc                                                                                                                                      >  3:35567/2‑147 (MQ=255)
           aCCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGc                                                                                                                                      >  1:145555/1‑147 (MQ=255)
           aCCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGc                                                                                                                                      >  8:83037/1‑147 (MQ=255)
           aCCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGc                                                                                                                                      >  1:316127/1‑147 (MQ=255)
             ccAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCa                                                                                                                                    >  1:110640/1‑147 (MQ=255)
             ccAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTAAGGTTGTGCCACAACCGGTGAAACGCCcga                                                                                                                                  >  5:313175/1‑149 (MQ=255)
                   ggACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCg                                                                                                                            >  6:221824/1‑149 (MQ=255)
                     aCGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGtt                                                                                                                          <  3:128155/149‑1 (MQ=255)
                      cGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTg                                                                                                                         <  8:179305/149‑1 (MQ=255)
                        aTGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCa                                                                                                                                    <  5:329040/136‑1 (MQ=255)
                        aTGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCa                                                                                                                                    >  6:329040/1‑136 (MQ=255)
                        aTGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGgc                                                                                                                       >  5:173080/1‑149 (MQ=255)
                         tGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAACCGGGGTGTCGTTTCCCTGGGTCCACACGGCGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGgcg                                                                                                                      >  2:75858/1‑149 (MQ=255)
                            ccGCAAAGCCGCGCAGGAAATCGTTCATAGCGGCAAGGCCAAACGGGTTGTCGGTTCCCTGGGTCCACAATGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTTCCACCATCGGTGAAAAGCCAGAGTACCGTTGgcgcgg                                                                                                                   <  7:228516/149‑3 (MQ=255)
                              gCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGt                                                                                                                 >  3:163707/1‑149 (MQ=255)
                              gCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCGGGGTGGTGATAGTGAAAACGGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGt                                                                                                                 >  6:29446/1‑149 (MQ=255)
                              gCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGGCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTAGGCGCTGGt                                                                                                                 >  6:235850/1‑149 (MQ=255)
                                aaaGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGa                                                                                                               <  7:130607/149‑1 (MQ=255)
                                            gAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCa                                                                                                          >  5:44772/1‑142 (MQ=255)
                                            gAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCa                                                                                                          <  6:44772/142‑1 (MQ=255)
                                            gAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCgg                                                                                                   >  6:10404/1‑149 (MQ=255)
                                             aaaTCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGgc                                                                                                  <  2:145555/149‑1 (MQ=255)
                                              aaTCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGgcg                                                                                                 <  4:248445/149‑1 (MQ=255)
                                                  gTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGAt                                                                                             <  7:83037/149‑1 (MQ=255)
                                                  gTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGa                                                                                              <  5:221824/148‑1 (MQ=255)
                                                  gTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGa                                                                                              <  8:144282/148‑1 (MQ=255)
                                                  gTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGGCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGa                                                                                              <  5:296227/148‑1 (MQ=255)
                                                   ttAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGAt                                                                                             >  5:179106/1‑148 (MQ=255)
                                                   ttAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGAt                                                                                             >  1:19869/1‑148 (MQ=255)
                                                     aaTAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTGGGCGCTGGTGACAGCATGGTCGGCGCGATGa                                                                                           <  6:255222/148‑1 (MQ=255)
                                                     aaTAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGa                                                                                           <  6:173080/148‑1 (MQ=255)
                                                      aTAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGaca                                                                                         <  4:163707/149‑1 (MQ=255)
                                                                 gCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACAATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTgg                                                                               >  8:66047/1‑148 (MQ=255)
                                                                  ccAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCa                                                                             <  3:277921/149‑1 (MQ=255)
                                                                          gTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGc                                                                     <  2:113510/149‑1 (MQ=255)
                                                                                 tttCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAACTGCCTCTCtt                                                              >  4:201818/1‑149 (MQ=255)
                                                                                     ccTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGaa                                                           <  4:169118/148‑1 (MQ=255)
                                                                                      cTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGaa                                                           >  7:97421/1‑147 (MQ=255)
                                                                                      cTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGaa                                                           >  7:101996/1‑147 (MQ=255)
                                                                                              acaAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTc                                                 <  2:110640/149‑1 (MQ=255)
                                                                                                           ggTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTa                                      <  5:192731/147‑1 (MQ=255)
                                                                                                            gtgtTGATATTTTAAAGTGTATTGAGGTGGTGGGTCCTTCGGTGTTTTGCGTGAGTAGCGTTGGCGCTGGTGTGAGCATGGTCGGCGCGATGTCACTGATACTGGCAGAAAATGCCTCTCTTGAAGTGATGGTTCGTTTTGGCGTAGCt                                   <  6:246729/149‑1 (MQ=255)
                                                                                                              gtTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGc                                                                     <  6:280325/113‑1 (MQ=255)
                                                                                                               ttGATAGTGAAAACTGTATTCAGGTGGTGCCACAACCGGTGAAAAGCCAGATTACCGTTGTCGCTGGTGACATCATGGTCTGCGCGATGACACTGAAACTGGAAGAAAATGCGTATCTTTAAGAGATGGTTCGATTTGGCGTAGCTGCg                                >  1:57170/1‑149 (MQ=255)
                                                                                                                  aTAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCggg                              <  2:316127/148‑1 (MQ=255)
                                                                                                                  aTAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgg                               <  1:318579/147‑1 (MQ=255)
                                                                                                                  aTAGTGAAAACTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCgg                               >  2:318579/1‑147 (MQ=255)
                                                                                                                   tAGTGAAAACTGTATTCAGGGGGTGCCTCCACCGGTGAAAAGCCAGTGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTGCGTTTTGGCGTAGCTGCGGGGa                            <  6:313175/149‑1 (MQ=255)
                                                                                                                    aGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATCGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAg                           >  2:5826/1‑149 (MQ=255)
                                                                                                                      tGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGa                                                                                           <  7:289344/83‑1 (MQ=255)
                                                                                                                      tGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGa                                                                                           >  8:289344/1‑83 (MQ=255)
                                                                                                                      tGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTg                         >  3:213142/1‑149 (MQ=255)
                                                                                                                        aaaaCTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAg                                                                            >  6:91093/1‑96 (MQ=255)
                                                                                                                        aaaaCTGTATTCAGGTGGTGCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAg                                                                            <  5:91093/96‑1 (MQ=255)
                                                                                                                          atcTGTATTCATGTGGTGCCACCATGGTTGAAAAGCCAGAGTACCGATGGCGCTGTTGACAGCATGGTCGGCGCGATGAGACTGAACCTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTTGCTGCGGGGAGTGCAg                      <  8:303195/146‑1 (MQ=255)
                                                                                                                            cTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCca                   <  5:10404/149‑1 (MQ=255)
                                                                                                                                 ttCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACt               >  1:210688/1‑148 (MQ=255)
                                                                                                                                    aGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGCCAGCATGGGCGGCCCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCaa            <  5:235850/148‑1 (MQ=255)
                                                                                                                                         gtgCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAgg       >  4:98136/1‑148 (MQ=255)
                                                                                                                                         gtgCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAgg       <  7:234541/148‑1 (MQ=255)
                                                                                                                                           gCCACCATCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGaa    >  4:148194/1‑149 (MQ=255)
                                                                                                                                            ccaccaCCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACATCATGGTCTGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCTGGTAGTGCAGCCACACTCAATCAGGGAac   <  1:75858/149‑1 (MQ=255)
                                                                                                                                             caccacCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACCCTCAATCAGGGAaca  >  1:264956/1‑149 (MQ=255)
                                                                                                                                             caccacCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAaca  >  1:260526/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATCGCGAACTCACCCAGCCGGACGATGTCCGCAAAGCCGCGCAGGAAATCGTTAATAGCGGCAAGGCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACA  >  NC_000913/1806956‑1807245

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: