Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,482,839 C→G 10.0% L167F (TTG→TTC emrK ← multidrug resistance efflux pump membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,482,8390CG10.0% 394.2 / 27.6 129L167F (TTG→TTCemrKmultidrug resistance efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (71/45);  new base G (8/5);  total (79/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

Only 100 of 130 total aligned reads displayed.

ACCGGACTCTTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTA  >  NC_000913/2482696‑2482988
                                                                                                                                               |                                                                                                                                                     
aCCGGACTCTTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGcc                                                                                                                                                 <  1:66482/150‑1 (MQ=255)
   ggACTCTTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTGAATGCCGCt                                                                                                                                              <  2:590411/150‑1 (MQ=255)
   ggACTCTTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCt                                                                                                                                              >  2:582521/1‑150 (MQ=255)
   ggACTCTTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGCTGTGTTCATTACTAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCtt                                                                                                                                             >  1:1131775/1‑150 (MQ=255)
     aCTCTTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTGAATGCCGCtat                                                                                                                                            <  1:800063/150‑3 (MQ=255)
      ctctTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTGAATGCCGCtatc                                                                                                                                           <  1:1521923/150‑4 (MQ=255)
        ctTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTATATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCt                                                                                                                                         >  2:386598/1‑150 (MQ=255)
        ctTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCGTCAACGCCTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTCTTCGCTTTATAAGCCTGGATAGCGGCATTGAATGCCGCtatcca                                                                                                                                         <  1:1095131/150‑6 (MQ=255)
          tAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTAc                                                                                                                                       <  1:1655966/150‑1 (MQ=255)
             tatCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGCCTTGTGGCGGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACtta                                                                                                                                    >  2:68278/1‑150 (MQ=255)
              atCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACttat                                                                                                                                   >  2:1394862/1‑150 (MQ=255)
                 cGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTaa                                                                                                                                >  2:1425642/1‑150 (MQ=255)
                  gTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAAc                                                                                                                               >  1:1405810/1‑150 (MQ=255)
                  gTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAAc                                                                                                                               <  1:1624089/150‑1 (MQ=255)
                   tACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTATTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACg                                                                                                                              <  1:1557680/150‑1 (MQ=255)
                     cGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTa                                                                                                                            >  1:951622/1‑150 (MQ=255)
                     cGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTAt                                                                                                                           <  1:566622/150‑1 (MQ=255)
                          aaGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCttt                                                                                                                       >  2:1560407/1‑150 (MQ=255)
                            gcgcCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACTTATCTTTgg                                                                                                                     >  1:1477442/1‑150 (MQ=255)
                             cgcCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGt                                                                                                                    >  2:496501/1‑150 (MQ=255)
                              gcCAACCAGGTTTCTTTAGTTGCATCCGACTCTTCAACGACTTGTGGCTCAATGTTTAATGGTGTGTTCATTACTAAAGCTTTAATCTCTTTATAAGACTGGATAGCGTCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTa                                                                                                                   <  1:269776/150‑1 (MQ=255)
                               ccAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTAt                                                                                                                  <  2:1158194/150‑1 (MQ=255)
                                 aaCCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGc                                                                                                                >  1:156145/1‑150 (MQ=255)
                                        cttcttTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCg                                                                                                         >  1:762427/1‑150 (MQ=255)
                                        cttcttTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCg                                                                                                         >  1:1478103/1‑150 (MQ=255)
                                        cttcttTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCg                                                                                                         >  1:570608/1‑150 (MQ=255)
                                        cttcttTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATACAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATTCTCAAGCg                                                                                                         >  2:639868/1‑150 (MQ=255)
                                             ttAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGCCGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCt                                                                                                    >  2:1538018/1‑150 (MQ=255)
                                                gTTGCATCCGCCGCTTCAACGTCTTGTGGCTGACGGTTTTATGGGGTGGTCATTACTAAAGCTTTATTCGCGTTCTAAGCCTGGATAGCGGCATTCACTGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCtttt                                                                                                 <  2:1071905/150‑1 (MQ=255)
                                                gTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCtttt                                                                                                 >  2:705755/1‑150 (MQ=255)
                                                gTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCtttt                                                                                                 <  1:1411529/150‑1 (MQ=255)
                                                 ttGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTg                                                                                                <  2:294268/150‑1 (MQ=255)
                                                 ttGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTg                                                                                                <  2:945303/150‑1 (MQ=255)
                                                 ttGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTg                                                                                                <  2:386165/150‑1 (MQ=255)
                                                  tGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTAGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGa                                                                                               <  2:1457572/150‑1 (MQ=255)
                                                   gCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGaa                                                                                              >  2:421103/1‑150 (MQ=255)
                                                     aTCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAAt                                                                                            >  2:1708219/1‑150 (MQ=255)
                                                        cgccgcTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAAc                                                                                         >  1:227465/1‑150 (MQ=255)
                                                         gccgcTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAAcc                                                                                        >  1:401443/1‑150 (MQ=255)
                                                           cgcTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAACGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTAGAAATACcccc                                                                                      >  1:333089/1‑150 (MQ=255)
                                                           cgcTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGATATAAcccc                                                                                      >  2:1394453/1‑150 (MQ=255)
                                                                cAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATATTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCACCTGCt                                                                                 >  2:1321162/1‑150 (MQ=255)
                                                                 aaCGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCtt                                                                                >  1:931867/1‑150 (MQ=255)
                                                                     aCTTGTGGCTGACGGTTTATTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCt                                                                            <  2:1578473/150‑1 (MQ=255)
                                                                      cTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTa                                                                           >  1:1595119/1‑150 (MQ=255)
                                                                      attGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTa                                                                           <  1:1733080/149‑1 (MQ=255)
                                                                         gtgGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACg                                                                        >  1:186974/1‑150 (MQ=255)
                                                                              tGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACt                                                                   >  2:413434/1‑150 (MQ=255)
                                                                              tGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACt                                                                   >  1:807111/1‑150 (MQ=255)
                                                                                       aaTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGttt                                                          >  2:604034/1‑149 (MQ=255)
                                                                                        aTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTtat                                                         >  1:1268956/1‑150 (MQ=255)
                                                                                         tGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTtata                                                        >  2:590672/1‑150 (MQ=255)
                                                                                            tgtgtTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAAtc                                                     <  1:542693/150‑1 (MQ=255)
                                                                                            tgggtTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTAGGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAAtc                                                     <  2:1091557/147‑1 (MQ=255)
                                                                                             gtgtTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAAtct                                                    >  2:1192459/1‑150 (MQ=255)
                                                                                             gtgtTCATTACTAAAGCTTTATTCGCATTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGATACTTATTAACGAATCTTTGGTATGCTCCAGCGTTTCTTTTGAATAAACCCACTGCTACGCTAACGGCTCTCGACGGTTATAAtct                                                    >  1:943495/1‑150 (MQ=255)
                                                                                                ttCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTAGCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAAtcttct                                                 <  1:248581/150‑1 (MQ=255)
                                                                                                    ttACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAg                                             <  1:1269796/150‑1 (MQ=255)
                                                                                                    ttACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATAAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAg                                             <  1:1058221/150‑1 (MQ=255)
                                                                                                     tACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGa                                            >  2:1358295/1‑150 (MQ=255)
                                                                                                     tACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGa                                            <  2:518356/150‑1 (MQ=255)
                                                                                                       cTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGAtt                                          <  1:1644151/150‑1 (MQ=255)
                                                                                                       cTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCAAACAATGCCGCTTTGCTACTTATTAACGTATCTTTTGTATGCTCCAGCGTTTCTTTAGAAAGAACCCCCAGCTACGCTAACGACACTCGACGGTTATAATCTTCTAAAGAtt                                          >  1:292971/1‑150 (MQ=255)
                                                                                                           aGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGAttgttg                                      <  1:429575/150‑1 (MQ=255)
                                                                                                            gCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGa                                     >  2:33492/1‑150 (MQ=255)
                                                                                                            gCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGa                                     >  1:555818/1‑150 (MQ=255)
                                                                                                             cTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGAt                                    <  2:126417/150‑1 (MQ=255)
                                                                                                             cTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGAt                                    >  1:1715003/1‑150 (MQ=255)
                                                                                                              tttATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATa                                   <  2:1312845/150‑1 (MQ=255)
                                                                                                               ttATTCGCTTTATAAGCATGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATaa                                  >  2:1604241/1‑149 (MQ=255)
                                                                                                                tATTCGCTTTATAATCCTGGATAGCGGCATTCAATGCCGCTTTGCTAATTATTAAAGTATCTTGGTTATGCTCCAGCGTTTCTTTTGAAATAACGCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACt                                 <  2:265658/150‑1 (MQ=255)
                                                                                                                  ttCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTAGTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATTACCCCCTGCTTCGCGAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGa                               >  1:294645/1‑150 (MQ=255)
                                                                                                                  ttCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGa                               >  1:1081676/1‑150 (MQ=255)
                                                                                                                    cGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGaaa                             >  1:1508650/1‑149 (MQ=255)
                                                                                                                    cGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAAt                             >  1:846783/1‑150 (MQ=255)
                                                                                                                    cGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAAt                             <  2:167587/150‑1 (MQ=255)
                                                                                                                     gCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATa                            >  2:954073/1‑150 (MQ=255)
                                                                                                                     gCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATa                            <  1:421103/150‑1 (MQ=255)
                                                                                                                        ttATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGt                         >  1:1357040/1‑150 (MQ=255)
                                                                                                                         tataAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGAGTTTCTTTTTAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTg                        >  1:379332/1‑150 (MQ=255)
                                                                                                                          ataAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGc                       >  2:1587701/1‑150 (MQ=255)
                                                                                                                               gtagcaaaGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCATTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAg                  >  1:673226/8‑150 (MQ=255)
                                                                                                                               gtagcaaaGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCATTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAg                  >  1:667172/8‑150 (MQ=255)
                                                                                                                                tagcaaaGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGc                 >  2:1565089/7‑150 (MQ=255)
                                                                                                                                  ggATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTGCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTTACGGCACTCGTCGGTTATAATCTTGTAAAGATTGTTGATACTGAATACGTGCTGAAGCGa               >  1:212183/1‑150 (MQ=255)
                                                                                                                                  ggATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGa               >  1:170518/1‑150 (MQ=255)
                                                                                                                                   gATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGAc              <  1:1560407/150‑1 (MQ=255)
                                                                                                                                   caaaGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGaa              >  1:190344/4‑149 (MQ=255)
                                                                                                                                   caaaGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCATTCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACTTGATGAAGAGAc              >  2:1503027/4‑150 (MQ=255)
                                                                                                                                    aaaGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACt             >  2:1185002/3‑150 (MQ=255)
                                                                                                                                         ggCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGc        >  1:886556/1‑150 (MQ=255)
                                                                                                                                         ggCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGc        <  1:386598/150‑1 (MQ=255)
                                                                                                                                         ggCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGc        >  2:760246/1‑150 (MQ=255)
                                                                                                                                           cATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCAc      >  2:490824/1‑150 (MQ=255)
                                                                                                                                            aTTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACt     >  2:658942/1‑150 (MQ=255)
                                                                                                                                            aTTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACt     >  1:1584278/1‑150 (MQ=255)
                                                                                                                                            aTTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACt     >  2:763985/1‑150 (MQ=255)
                                                                                                                                            aTTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACt     <  2:951622/150‑1 (MQ=255)
                                                                                                                                            aTTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCCGAAGTCGCTTCAGCAcg     <  1:1039884/150‑2 (MQ=255)
                                                                                                                                             ttCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAATGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTg    <  1:564454/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
ACCGGACTCTTAATATCCGTACGTTTAAGCGCCAACCAGGCTTCTTTAGTTGCATCCGCCGCTTCAACGACTTGTGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATACGTGCTGAAGCGACTTCGGCACTGTA  >  NC_000913/2482696‑2482988

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: