Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 884,407 C→G 6.5% I245M (ATC→ATG mdfA → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913884,4070CG 6.5% 236.3 / 4.9 74I245M (ATC→ATGmdfAmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (33/36);  new base G (2/3);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGAT  >  NC_000913/884258‑884556
                                                                                                                                                     |                                                                                                                                                     
gCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTatc                                                                                                                                                       >  2:848729/1‑150 (MQ=255)
 ccACGCGTATTGGCGAGAAACTGTCACTGAAAGAACTCGTTCGAGACTATAAGCTGGTGCCGAAGACCGCCCGCTATGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTatca                                                                                                                                                      <  1:110781/150‑1 (MQ=255)
  cccgcgTATCGGCGAGACACTGTCCCTGACAGAACTCGGTCGGGACTATAAGCTGCTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGGGGATCGCCCAGTCGCCGATTatcat                                                                                                                                                     <  1:1100024/148‑1 (MQ=255)
        tataGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATATCCCAGTCGCCGATTATCAACATTAc                                                                                                                                               >  2:646485/1‑150 (MQ=255)
         ataGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTAcc                                                                                                                                              >  2:993127/1‑150 (MQ=255)
              cGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCgc                                                                                                                                         >  1:1448568/1‑149 (MQ=255)
               gagaAACTGTCACTGAACGAACTCGTGCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGTCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAg                                                                                                                                        <  2:714696/150‑1 (MQ=255)
                agaAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGc                                                                                                                                       <  1:308108/150‑1 (MQ=255)
                   aaCTGTCACTGAAAGAACTCGGGCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGGTAATGATGATAATCGGCGactggg                                                                                                                                    <  2:1155822/150‑6 (MQ=255)
                     cTGTCACTTACAGAACTCGGTCGTGACTATAAGCTGGCGCTGAAGCCCGGCCGCGTTGGGGCGGGGGCGCTGGCGCGGGGATTCGTTAGTCTGCCGTTGCGGGCGTGGAGCGCCCAGTCGCCGGTAATGATGATAATCGGCGactgggcg                                                                                                                                  <  2:250030/150‑8 (MQ=255)
                      tGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGa                                                                                                                                 <  2:940428/150‑1 (MQ=255)
                       gTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGag                                                                                                                                >  1:1593365/1‑150 (MQ=255)
                         cACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGagca                                                                                                                              <  1:1460985/150‑1 (MQ=255)
                         cACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGagca                                                                                                                              >  1:1026419/1‑150 (MQ=255)
                          aCTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGagcag                                                                                                                             >  2:1627284/1‑150 (MQ=255)
                            tGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCt                                                                                                                           >  1:1344814/1‑150 (MQ=255)
                            tGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCAACATTACCGGCGAGCAGTTGAGCAGCt                                                                                                                           >  1:770544/1‑150 (MQ=255)
                               aaGAACTCGGTCGCGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATg                                                                                                                        <  1:993127/150‑1 (MQ=255)
                                aGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGa                                                                                                                       <  2:1454464/150‑1 (MQ=255)
                                 gAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCGGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGaa                                                                                                                      <  1:595306/150‑1 (MQ=255)
                                      cGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATgg                                                                                                                 >  2:1023240/1‑150 (MQ=255)
                                      cGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGGCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATgg                                                                                                                 <  2:169369/150‑1 (MQ=255)
                                           gTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGAAGCTATGAATATGGCTtgc                                                                                                            >  1:981810/1‑150 (MQ=255)
                                           gTGACTATAAGCTGGCGCTGAAGAGCGGCCGCGTGGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTTTGATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTtgc                                                                                                            <  1:76380/150‑1 (MQ=255)
                                             gACTATAAGCTGGTGCTGAAGAACCGCCGCTTTGTGGCGGGGGCGCGGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTtgctg                                                                                                          <  2:740766/150‑1 (MQ=255)
                                                 ataAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAg                                                                                                      <  1:1345367/150‑1 (MQ=255)
                                                    aGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGc                                                                                                   <  1:82871/150‑1 (MQ=255)
                                                    aGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGc                                                                                                   <  1:308365/150‑1 (MQ=255)
                                                    aGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGc                                                                                                   >  1:725769/1‑150 (MQ=255)
                                                        ggTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCCATTATCATCATTACCGCCGCGCAGTTGAGCTGCTATGAAGATGGCTTGCTGCAAGTGCCGAt                                                                                               >  2:281443/1‑150 (MQ=255)
                                                              gaagaaCGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCgg                                                                                         >  2:1507512/1‑150 (MQ=255)
                                                               aagaaCGGCAGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCTTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCggg                                                                                        <  1:272701/150‑1 (MQ=255)
                                                                    cGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGtt                                                                                   <  2:1615135/150‑1 (MQ=255)
                                                                         ggttgggggCGGGGGCGCGGGCGCTGGGCTTCGGTAGTCTGCCGTTGCTGGCGTGTATCGCGCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAATTGCCTATTTTCGGGGCGTTAATTg                                                                              <  1:830196/143‑1 (MQ=255)
                                                                         gCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTg                                                                              >  1:1117663/1‑150 (MQ=255)
                                                                         gCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGAAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTg                                                                              >  2:1364664/1‑150 (MQ=255)
                                                                            ttGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCgg                                                                           <  1:1677607/150‑1 (MQ=255)
                                                                             tgtgGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCggg                                                                          <  1:787518/150‑1 (MQ=255)
                                                                             tgtgGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCggg                                                                          <  2:638257/150‑1 (MQ=255)
                                                                                  cGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACt                                                                     <  1:305633/150‑1 (MQ=255)
                                                                                     gggCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGc                                                                  >  1:531763/1‑150 (MQ=255)
                                                                                     gggCGCTGGCGCTGGGAATCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTATCGGGGCGTTAATTGCGGGTAACTTGc                                                                  <  1:89549/150‑1 (MQ=255)
                                                                                      ggcgctggcgctTGGAGTGGTTAGTCTGCGGTTGCTGGCGTGGATCGCCCAGTCGCAGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTTTGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCt                                                                 <  1:437489/150‑1 (MQ=255)
                                                                                       gcgctggcgctgGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTg                                                                <  2:962019/150‑1 (MQ=255)
                                                                                       gcgctggcgctgGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTg                                                                <  1:585681/150‑1 (MQ=255)
                                                                                       gcgctcgcgctgGGACTCGTGAGTCTGCACGCCCGCGCGTTGATCGCCCAGGCGCCGATCATCATCCTTACCGGCGAGCAGCCGAGCACCTATGAATATGGCTTGCTGCAAGTGCCTCTTTTCCCGGCCTTAAACGCGGGTAACTTGCTg                                                                <  2:80274/150‑1 (MQ=255)
                                                                                           tggcgctggGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAg                                                            >  2:1606690/1‑150 (MQ=255)
                                                                                           tggcgctggGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAg                                                            >  2:1613093/1‑150 (MQ=255)
                                                                                               gctggGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTATTTGCGGGTAACTTGCTGTTATcgcg                                                        >  2:1726862/1‑150 (MQ=255)
                                                                                                 tggGATTCGGTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTc                                                      <  1:848729/150‑1 (MQ=255)
                                                                                                    gATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGa                                                   >  1:1526312/1‑150 (MQ=255)
                                                                                                    gATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGa                                                   >  1:176404/1‑150 (MQ=255)
                                                                                                     agtCGTTAGGCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCTAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGAc                                                  <  2:1116508/148‑1 (MQ=255)
                                                                                                       tCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGGGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCt                                                <  1:581790/150‑1 (MQ=255)
                                                                                                             gTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGcgcc                                          >  1:1647169/1‑150 (MQ=255)
                                                                                                               cTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGcgccgc                                        >  1:1569308/1‑150 (MQ=255)
                                                                                                                 gCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCAc                                      >  2:1357241/1‑150 (MQ=255)
                                                                                                                    gTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGt                                   <  2:978284/150‑1 (MQ=255)
                                                                                                                       gCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGACGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTCGCGCGTCTGACCTCGCGCCGCACCGGACg                                >  1:999741/1‑150 (MQ=255)
                                                                                                                          ggCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTc                             >  2:73250/1‑150 (MQ=255)
                                                                                                                           gCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCg                            <  1:620156/150‑1 (MQ=255)
                                                                                                                           gCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCg                            <  1:441796/150‑1 (MQ=255)
                                                                                                                               ggATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGa                        >  1:1280104/1‑150 (MQ=255)
                                                                                                                               ggATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGGCTGACCTCGCGCCGCACCGTACGTTCGCTGa                        >  1:56609/1‑150 (MQ=255)
                                                                                                                                gATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGat                       >  1:387889/1‑150 (MQ=255)
                                                                                                                                    gCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGattatt                   >  1:402580/1‑150 (MQ=255)
                                                                                                                                     cccAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGAttatta                  >  1:1508352/1‑150 (MQ=255)
                                                                                                                                     actgctCGCCGGTAATGATGATAATCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGAttatta                  >  1:560954/6‑150 (MQ=255)
                                                                                                                                     actgctCGCCGGTAATGATGATAATCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGAttatta                  >  2:367847/6‑150 (MQ=255)
                                                                                                                                           cGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGGTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTGCGCTGATTATTATGggcg            <  1:995330/150‑1 (MQ=255)
                                                                                                                                            gCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGggcgg           >  1:350720/1‑150 (MQ=255)
                                                                                                                                             ccGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGggcggc          <  1:1627284/150‑1 (MQ=255)
                                                                                                                                             ccGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGggcggc          <  1:407097/150‑1 (MQ=255)
                                                                                                                                             ccGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGggcggc          <  2:1448568/150‑1 (MQ=255)
                                                                                                                                             ccGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGGCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGggcggc          <  2:1321907/150‑1 (MQ=255)
                                                                                                                                                aTTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTgg       >  1:927894/1‑150 (MQ=255)
                                                                                                                                                  tATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGcc     <  2:725769/150‑1 (MQ=255)
                                                                                                                                                   atcatcATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTTCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGAAGGCCg    >  1:224405/1‑150 (MQ=255)
                                                                                                                                                   atcatcATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGAGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCg    >  1:935104/1‑150 (MQ=255)
                                                                                                                                                     catcatTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGAt  >  1:1490457/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGAT  >  NC_000913/884258‑884556

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