Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 884,410 C→G 5.4% I246M (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,4100CG 5.4% 235.8 / 3.1 74I246M (ATC→ATGmdfAmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base G (2/2);  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 = 9.27e-01

TATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTAT  >  NC_000913/884266‑884559
                                                                                                                                                |                                                                                                                                                     
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:1026419/1‑150 (MQ=255)
                 cACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGagca                                                                                                                                 <  1:1460985/150‑1 (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:308365/150‑1 (MQ=255)
                                            aGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGc                                                                                                      <  1:82871/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                                                                   <  1:585681/150‑1 (MQ=255)
                                                                               gcgctggcgctgGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTg                                                                   <  2:962019/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:176404/1‑150 (MQ=255)
                                                                                            gATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGa                                                      >  1:1526312/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                     >  2:367847/6‑150 (MQ=255)
                                                                                                                             actgctCGCCGGTAATGATGATAATCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGAttatta                     >  1:560954/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)
                                                                                                                                              atcatTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCAACGTACGTTCGCTGATTATTATGGGCGGCTGGCCGAtt    >  2:1161940/1‑150 (MQ=255)
                                                                                                                                               tcatTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTa   >  1:1457684/1‑150 (MQ=255)
                                                                                                                                                catTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTAt  >  2:911214/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTAT  >  NC_000913/884266‑884559

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: