Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,807,169 C→T 5.7% A267V (GCA→GTA)  pfkB → 6‑phosphofructokinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,807,1690CT 5.7% 228.2 / 2.7 70A267V (GCA→GTA) pfkB6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base C (36/30);  new base T (2/2);  total (38/32)
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.19e-01

GCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGC  >  NC_000913/1807021‑1807317
                                                                                                                                                    |                                                                                                                                                    
gCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGc                                                                                                                                                      <  8:277245/149‑1 (MQ=255)
 ccAAACGGGTTGTCGTTTCCCTTGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGGGCCACCACCGGTGAAAAGCCAGAGTACCGTGGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCa                                                                                                                                                     <  6:178668/149‑1 (MQ=255)
         gTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGc                                                                                                                                             >  3:281930/1‑149 (MQ=255)
         gTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGc                                                                                                                                             >  4:72576/1‑149 (MQ=255)
         gTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGc                                                                                                                                             <  5:18653/149‑1 (MQ=255)
            gTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCct                                                                                                                                           <  6:398859/148‑1 (MQ=255)
             tCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCctct                                                                                                                                         >  1:370883/1‑149 (MQ=255)
                tttCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCCCCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCtt                                                                                                                                      >  1:326835/1‑149 (MQ=255)
                  tCCCTGGGGCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGACACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg                                                                                                                                     <  7:93415/148‑1 (MQ=255)
                    ccTGTGTCCAGAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGTTGGTGGCACCTCCTGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGaa                                                                                                                                   <  6:236465/148‑1 (MQ=255)
                    ccTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTg                                                                                                                                     >  6:293600/1‑146 (MQ=255)
                     cTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAaga                                                                                                                                 <  2:152013/149‑1 (MQ=255)
                            cacaAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtt                                                                                                                          <  5:158078/149‑1 (MQ=255)
                                aGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGtt                                                                                                                      >  2:297156/1‑149 (MQ=255)
                                    gcgcTGTGTGTTGATAGTGAAAACTGTATTCAGGTGTTGCCACCTCCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATTGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATTGTTCGTTTTgg                                                                                                                  <  3:208881/149‑1 (MQ=255)
                                             gtTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGc                                                                                                         <  8:419584/149‑1 (MQ=255)
                                             gtTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGc                                                                                                         >  3:240682/1‑149 (MQ=255)
                                                gATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCggg                                                                                                      >  2:128255/1‑149 (MQ=255)
                                                 aTAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCggg                                                                                                      <  5:293600/148‑1 (MQ=255)
                                                  tAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGa                                                                                                    <  3:351862/149‑1 (MQ=255)
                                                    gTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGt                                                                                                  <  4:240682/149‑1 (MQ=255)
                                                     tGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTg                                                                                                 >  7:27890/1‑149 (MQ=255)
                                                     tGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTg                                                                                                 >  3:181650/1‑149 (MQ=255)
                                                        aaaCTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAg                                                                                              <  1:207362/149‑1 (MQ=255)
                                                        aaaCTGTATTCAGGTGGTGCCACCACAGGTGAAAAGCCAGAGTACCGTTGGCGCAGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCGTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGGGGGGAGTGCAg                                                                                              >  6:134646/1‑149 (MQ=255)
                                                              tATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCcaca                                                                                         >  4:17914/1‑148 (MQ=255)
                                                                  cAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCa                                                                                     <  2:370883/148‑1 (MQ=255)
                                                                   aGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGTAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGa                                                                                                    <  4:15352/132‑1 (MQ=255)
                                                                   aGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGTAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGa                                                                                                    <  4:20930/132‑1 (MQ=255)
                                                                   aGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGTAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGa                                                                                                    >  3:20930/1‑132 (MQ=255)
                                                                   aGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGTAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGa                                                                                                    >  3:15352/1‑132 (MQ=255)
                                                                        gtgCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCca                                                                                           <  8:141971/136‑1 (MQ=255)
                                                                        gtgCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCca                                                                                           >  7:141971/1‑136 (MQ=255)
                                                                        gtgCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAgg                                                                               >  7:23233/1‑148 (MQ=255)
                                                                        gtgCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAgg                                                                               <  1:128255/148‑1 (MQ=255)
                                                                        gtgCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGAGAGTGCAGCCACACTCAATCAgg                                                                               <  3:72576/148‑1 (MQ=255)
                                                                            caccacCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAaca                                                                          <  3:206531/149‑1 (MQ=255)
                                                                              ccaccGGTGAAAAGCCAGAGTACCGTTGGCGCTCGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACg                                                                        >  1:223661/1‑149 (MQ=255)
                                                                               caccGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGt                                                                       >  4:403647/1‑149 (MQ=255)
                                                                                 ccGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAacacagac                                                                     <  4:281930/149‑5 (MQ=255)
                                                                                      gAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGGCACACTCAATCAGGGAACACGTCTGTGCt                                                                >  3:174437/1‑149 (MQ=255)
                                                                                            ccAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGa                                                                                                    <  4:36942/107‑1 (MQ=255)
                                                                                            ccAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGa                                                                                                    >  3:36942/1‑107 (MQ=255)
                                                                                            ccAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                                           >  8:361633/1‑148 (MQ=255)
                                                                                            ccAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                                           <  7:361633/148‑1 (MQ=255)
                                                                                             cAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAaca                                                                          <  6:29137/132‑1 (MQ=255)
                                                                                             cAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAaca                                                                          >  5:29137/1‑132 (MQ=255)
                                                                                              agagTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGAc                                                        >  7:128642/1‑149 (MQ=255)
                                                                                              agagTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGAc                                                        <  3:17914/149‑1 (MQ=255)
                                                                                                agTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGa                                                      >  1:155261/1‑149 (MQ=255)
                                                                                                agTACCGTTGGCGCTGGTGACAGCATGGGCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGGCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGa                                                      <  2:72940/149‑1 (MQ=255)
                                                                                                           cgcTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCaaaaaa                                           <  8:27890/149‑1 (MQ=255)
                                                                                                            gcTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTTGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAt                                          >  2:125867/1‑149 (MQ=255)
                                                                                                                gTGACAGCATGGGCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAttt                                        <  1:297156/147‑1 (MQ=255)
                                                                                                                  gACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGc                                    >  7:336191/1‑149 (MQ=255)
                                                                                                                     aGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTa                                 >  6:402706/1‑149 (MQ=255)
                                                                                                                     aGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTa                                 >  8:136910/1‑149 (MQ=255)
                                                                                                                         tGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCtt                             >  1:134522/1‑149 (MQ=255)
                                                                                                                         tGGTCGGCGCAATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTGCGTTTTGGCGTAGCTGCGGGGAGTGCAGGCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCtt                             >  1:46689/1‑149 (MQ=255)
                                                                                                                          ggTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt                            >  7:138508/1‑149 (MQ=255)
                                                                                                                            tCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACGTTTcc                          >  2:202351/1‑149 (MQ=255)
                                                                                                                               gcgcGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGGCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGc                       <  6:386452/149‑1 (MQ=255)
                                                                                                                                  cgATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTaa                    <  2:134522/149‑1 (MQ=255)
                                                                                                                                      gACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGGCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACaa                 >  4:73641/1‑148 (MQ=255)
                                                                                                                                         acTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTATCaaaaa              >  3:282269/1‑148 (MQ=255)
                                                                                                                                          cTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACa            >  4:43000/1‑149 (MQ=255)
                                                                                                                                           tGAAACTGGCAGAAAATGCATCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTACCCGCTAACAAAAACAt           >  6:406428/1‑149 (MQ=255)
                                                                                                                                            gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                                           <  2:132401/100‑1 (MQ=255)
                                                                                                                                            gAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                                           >  1:132401/1‑100 (MQ=255)
                                                                                                                                             aaaCTGGCAGAAAATGCCTCTCTTGAAGAGCTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATAGGCAAAAAATTTACGCTTACATTTCCCGCTAACAAAAACATTc         >  7:386898/1‑149 (MQ=255)
                                                                                                                                             aaaCTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGGCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACAtt          <  2:46689/148‑1 (MQ=255)
                                                                                                                                                    cAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGc  <  3:43000/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCAAACGGGTTGTCGTTTCCCTGGGTCCACAAGGAGCGCTGGGTGTTGATAGTGAAAACTGTATTCAGGTGGTGCCACCACCGGTGAAAAGCCAGAGTACCGTTGGCGCTGGTGACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGC  >  NC_000913/1807021‑1807317

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