Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,806,390 G→A 5.2% L7L (TTG→TTA)  pfkB → 6‑phosphofructokinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,806,3900GA 5.2% 252.6 / 2.3 76L7L (TTG→TTA) pfkB6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base G (51/21);  new base A (2/2);  total (53/23)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TTTTTATTTATATATATTTATCTGCAAAATTTTAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCAC  >  NC_000913/1806250‑1806536
                                                                                                                                            |                                                                                                                                                  
tttttATTTATATATATTTACCTGCAAAATTTTAAATAAAGCTCCAATAAGTCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATTAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTgc                                                                                                                                            >  4:981110/1‑149 (MQ=255)
   ttatttatATATATTTATCTGCAAAATTTTAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGcc                                                                                                                                         <  2:1126646/149‑1 (MQ=255)
      tttatATATATTTATCTGCAAAATTTTAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTAACACTTGCGCCctc                                                                                                                                      >  3:723780/1‑149 (MQ=255)
               atTTATCTGCAAAATTTTAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATGAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAg                                                                                                                             <  4:1111336/149‑1 (MQ=255)
                      tGCAAAATTTTAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGcaaca                                                                                                                      >  1:1467442/1‑149 (MQ=255)
                      tGCAAAATTTTAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGcaaca                                                                                                                      >  4:154751/1‑149 (MQ=255)
                       gCAAAAGTTTAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATGCGGTGCCAGACTGAAATCAGCCTATCGGAGGAAAAGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGcaacaa                                                                                                                     <  1:1151911/149‑1 (MQ=255)
                        cAAAATTTTAACTAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAAt                                                                                                                    >  1:768374/1‑149 (MQ=255)
                         aaaaTTTTAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAAtt                                                                                                                   >  3:67695/1‑149 (MQ=255)
                             ttttAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTAccc                                                                                                               >  3:1073765/1‑149 (MQ=255)
                              tttAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTAcccc                                                                                                              <  4:1050674/149‑1 (MQ=255)
                               ttAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTAACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCg                                                                                                             >  1:664265/1‑149 (MQ=255)
                                taaataaaGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTTAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGc                                                                                                            >  1:277811/1‑149 (MQ=255)
                                  aataaaGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCaa                                                                                                          >  3:358699/1‑149 (MQ=255)
                                      aaGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCtctctc                                                                                                                                  <  4:633627/121‑1 (MQ=255)
                                      aaGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCtctctc                                                                                                                                  >  3:633627/1‑121 (MQ=255)
                                      aaGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAAttt                                                                                                      >  1:706708/1‑149 (MQ=255)
                                       aGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTa                                                                                                     >  2:1126570/1‑149 (MQ=255)
                                            cAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCg                                                                                                >  1:736996/1‑149 (MQ=255)
                                             aaTAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCTCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGa                                                                                               >  4:960160/1‑149 (MQ=255)
                                              aTAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGaa                                                                                              >  2:426021/1‑149 (MQ=255)
                                                 aaTCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATc                                                                                                   <  1:1153292/141‑1 (MQ=255)
                                                 aaTCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATc                                                                                                   >  2:1153292/1‑141 (MQ=255)
                                                    cATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGaaaa                                                                                        <  2:1401776/149‑1 (MQ=255)
                                                     atatTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGAt                                                                                                                               <  2:485686/109‑1 (MQ=255)
                                                     atatTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGAt                                                                                                                               >  1:485686/1‑109 (MQ=255)
                                                     atatTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAAc                                                                                       <  3:923301/149‑1 (MQ=255)
                                                      tatTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACt                                                                                      >  1:779726/1‑149 (MQ=255)
                                                        ttGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTgc                                                                                    >  2:726503/1‑149 (MQ=255)
                                                                   ttcACTTTCCTCTGATTCGGTGCCAGACTTAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCa                                                                         <  1:1251402/149‑1 (MQ=255)
                                                                   ttcACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCa                                                                         >  3:500052/1‑149 (MQ=255)
                                                                     cACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCAcc                                                                       >  3:61973/1‑149 (MQ=255)
                                                                       cTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCgg                                                                     >  3:327843/1‑149 (MQ=255)
                                                                        tttCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGgt                                                                    >  1:334996/1‑149 (MQ=255)
                                                                            cGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTAACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTc                                                                <  2:664265/149‑1 (MQ=255)
                                                                             gCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCg                                                               >  3:252101/1‑149 (MQ=255)
                                                                                    cGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTAACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCgg                                                        <  4:723780/149‑1 (MQ=255)
                                                                                         ccAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCCAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGggcg                                                   >  4:336266/1‑149 (MQ=255)
                                                                                          cAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGggcgg                                                  >  3:1230999/1‑149 (MQ=255)
                                                                                             aCTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGggcggcgg                                               >  2:144786/1‑149 (MQ=255)
                                                                                                gAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCAt                                            >  2:1309644/1‑149 (MQ=255)
                                                                                                  aaTCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCa                                          >  3:1192454/1‑149 (MQ=255)
                                                                                                  aaTCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCa                                          >  1:676960/1‑149 (MQ=255)
                                                                                                   aTCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCaa                                         >  3:1113968/1‑149 (MQ=255)
                                                                                                   aTCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCaa                                         <  1:539250/149‑1 (MQ=255)
                                                                                                             aggaggAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCcgc                               >  1:91427/1‑149 (MQ=255)
                                                                                                               gaggaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCTGAAGGAAAACTgcgc                                                                                  <  3:261126/96‑1 (MQ=255)
                                                                                                               gaggaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCTGAAGGAAAACTgcgc                                                                                  >  4:261126/1‑96 (MQ=255)
                                                                                                                aggaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTTGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGcc                            >  4:870506/1‑149 (MQ=255)
                                                                                                                 ggaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCgggg                                                      >  2:470161/1‑122 (MQ=255)
                                                                                                                 ggaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCgggg                                                      <  1:470161/122‑1 (MQ=255)
                                                                                                                 ggaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGggcggc                                                 <  4:830018/127‑1 (MQ=255)
                                                                                                                 ggaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGggcggc                                                 >  3:830018/1‑127 (MQ=255)
                                                                                                                 ggaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCa                           >  4:272891/1‑149 (MQ=255)
                                                                                                                  gaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCAt                          >  3:803705/1‑149 (MQ=255)
                                                                                                                  gaAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCAt                          <  3:981110/149‑1 (MQ=255)
                                                                                                                    aaTGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTTGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTg                        >  4:151063/1‑149 (MQ=255)
                                                                                                                    aaTGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTg                        >  4:1447991/1‑149 (MQ=255)
                                                                                                                    aaTGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTg                        <  2:1467442/149‑1 (MQ=255)
                                                                                                                    aaTGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTg                        >  4:1191091/1‑149 (MQ=255)
                                                                                                                     atgatgGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGc                       <  2:334996/149‑1 (MQ=255)
                                                                                                                            tACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCtt                >  2:629073/1‑149 (MQ=255)
                                                                                                                            tACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCtt                <  3:154751/149‑1 (MQ=255)
                                                                                                                             aCGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTg               >  2:1108223/1‑149 (MQ=255)
                                                                                                                             aCGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTg               >  2:411702/1‑149 (MQ=255)
                                                                                                                               gTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTTGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGa             <  3:151063/149‑1 (MQ=255)
                                                                                                                               gTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGa             <  4:1192454/149‑1 (MQ=255)
                                                                                                                               gTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGa             >  1:125996/1‑149 (MQ=255)
                                                                                                                                tATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAg            >  2:123699/1‑149 (MQ=255)
                                                                                                                                 aTCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAgg           >  2:507739/1‑149 (MQ=255)
                                                                                                                                  tCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGc          >  1:98112/1‑149 (MQ=255)
                                                                                                                                    tataCGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAg        >  3:591080/1‑149 (MQ=255)
                                                                                                                                        cGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGcc    <  4:1073765/149‑1 (MQ=255)
                                                                                                                                          ttGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCcac  >  4:633001/1‑149 (MQ=255)
                                                                                                                                          ttGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCcac  >  4:829433/1‑149 (MQ=255)
                                                                                                                                          ttGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCcac  >  3:973632/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
TTTTTATTTATATATATTTATCTGCAAAATTTTAAATAAAGCTCCAATAAATCATATTGTTAATTTCTTCACTTTCCGCTGATTCGGTGCCAGACTGAAATCAGCCTATAGGAGGAAATGATGGTACGTATCTATACGTTGACACTTGCGCCCTCTCTCGATAGCGCAACAATTACCCCGCAAATTTATCCCGAAGGAAAACTGCGCTGTACCGCACCGGTGTTCGAACCCGGGGGCGGCGGCATCAACGTCGCCCGCGCCATTGCCCATCTTGGAGGCAGTGCCAC  >  NC_000913/1806250‑1806536

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