Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,743,611 C→A S52Y (TCT→TAT)  mdtK → multidrug efflux system transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,743,6110CA98.6% 245.4 / NA 69S52Y (TCT→TAT) mdtKmultidrug efflux system transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (29/39);  minor base T (1/0);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCTGGCTGGCAGGTTTTGTTTCCGTT  >  NC_000913/1743464‑1743756
                                                                                                                                                   |                                                                                                                                                 
aGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTatat                                                                                                                                                >  1:1256341/1‑151 (MQ=255)
aGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTatat                                                                                                                                                >  1:692926/1‑151 (MQ=255)
   atatCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTg                                                                                                                                             <  1:873560/151‑1 (MQ=255)
    tatCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGCACTTTTATcctg                                                                                                                                            >  2:984001/1‑148 (MQ=255)
    tatCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCAACGACATGGCGGCGGGCGCTATCGGGACTTATATCTg                                                                                                                                             >  2:1216182/1‑150 (MQ=255)
     atCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTgg                                                                                                                                            <  1:999905/150‑1 (MQ=255)
       cAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCGCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCt                                                                                                                                          <  1:282162/150‑1 (MQ=255)
       cAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCtt                                                                                                                                         <  2:1417764/151‑1 (MQ=255)
       cAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCt                                                                                                                                          <  2:1280780/150‑1 (MQ=255)
           gAAGCGCGTCTGTTATTAGCATTACCAATCCCGGTGATTCTCGCGCAAATCGCCAAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCTGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCgg                                                                                                                                     <  2:1298188/151‑1 (MQ=255)
                gcgTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATc                                                                                                                                <  2:523710/151‑1 (MQ=255)
                     tGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCtt                                                                                                                           >  1:869131/1‑151 (MQ=255)
                       ttattaGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGGCGGTATCGGGACTTATATCTGGCTTCCGGCGATCCTCttt                                                                                                                          >  2:1392563/1‑150 (MQ=255)
                         attagcattagcAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTgg                                                                                                                        <  1:909650/150‑1 (MQ=255)
                           tagcattagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTc                                                                                                                      <  1:1293915/150‑1 (MQ=255)
                            agcattagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCa                                                                                                                     <  1:773482/150‑1 (MQ=255)
                             gcattagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACg                                                                                                                   >  1:1161113/1‑151 (MQ=255)
                               attagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGa                                                                                                                 <  1:1314155/151‑1 (MQ=255)
                               attagcaATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCCATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCCCgg                                                                                                                  >  2:972332/1‑150 (MQ=255)
                                     aaTCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGActgctg                                                                                                           >  1:1413460/1‑151 (MQ=255)
                                         ccGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTgg                                                                                                       >  1:986671/1‑151 (MQ=255)
                                         ccGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTgg                                                                                                       >  2:435751/1‑151 (MQ=255)
                                         ccGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTgg                                                                                                       >  1:1538492/1‑151 (MQ=255)
                                            gTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTAAAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCAt                                                                                                    >  2:891101/1‑151 (MQ=255)
                                                 tctcGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACg                                                                                               >  1:1356610/1‑151 (MQ=255)
                                                    cgcgCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGc                                                                                              <  1:830538/149‑1 (MQ=255)
                                                      cgcAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCGGGt                                                                                          <  2:1532248/151‑1 (MQ=255)
                                                        cAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATATTGCCACCGAAATGGGGGCGGTAGTTATCGGTAGCTATATCTGGCTTTCCGCGGTACTGTTTTGTTACGGGCTATTGCTGGCATTTATGGCGGtt                                                                                         >  2:1110978/1‑150 (MQ=255)
                                                         aaaTCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTa                                                                                        <  1:275826/150‑1 (MQ=255)
                                                         aaaTCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTAt                                                                                       >  1:742277/1‑151 (MQ=255)
                                                           aTCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATcg                                                                                     <  2:63421/151‑1 (MQ=255)
                                                            tCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGGACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATcg                                                                                     >  2:1122763/1‑150 (MQ=255)
                                                             cGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATcgcg                                                                                   <  2:1501769/151‑1 (MQ=255)
                                                              gCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCgcgc                                                                                  <  1:493567/151‑1 (MQ=255)
                                                                 cAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTTGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAAAGCCGGTTATCGCGCAAt                                                                               >  2:1235841/1‑151 (MQ=255)
                                                                  aaaCTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAAtt                                                                              <  1:403316/151‑1 (MQ=255)
                                                                   aaCTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTa                                                                             <  2:517212/151‑1 (MQ=255)
                                                                          aTGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGtt                                                                      <  1:961751/151‑1 (MQ=255)
                                                                          aTGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGtt                                                                      <  1:1550823/151‑1 (MQ=255)
                                                                           tGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTc                                                                     >  1:883891/1‑151 (MQ=255)
                                                                               ttttGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCgg                                                                  <  1:3595/150‑1 (MQ=255)
                                                                                tttGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGt                                                                 <  1:1125388/150‑1 (MQ=255)
                                                                                 ttGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTc                                                                <  1:352580/150‑1 (MQ=255)
                                                                                 ttGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTc                                                                <  1:262020/150‑1 (MQ=255)
                                                                                   gTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGa                                                              <  2:766605/150‑1 (MQ=255)
                                                                                     cGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGt                                                           <  2:918052/151‑1 (MQ=255)
                                                                                     cGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGt                                                           >  1:508750/1‑151 (MQ=255)
                                                                                      gATACAGTGATGGCGGGCGGCTATAGTGCCACCGACATGGGGGCGGTTGCTATCGGGACTTATATCTGGCTTTCGGCGATCCTCTTTGGTCCCGGGCTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGGCGatggg                                                          >  2:92110/1‑147 (MQ=255)
                                                                                       aTACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGa                                                         <  1:174010/151‑1 (MQ=255)
                                                                                                    gggcggCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCAt                                            <  2:309463/151‑1 (MQ=255)
                                                                                                        ggcTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAgg                                        >  2:1198571/1‑151 (MQ=255)
                                                                                                        ggcTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAgg                                        <  1:491721/151‑1 (MQ=255)
                                                                                                         gcTATAGTGCCATCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGATGGCATTAAAGGCGGATATCGCGCAATTAAATGGTTCCGGTCGATGTGAGCGCATTGCGCATc                                           >  2:1148385/1‑147 (MQ=255)
                                                                                                         gcTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAgg                                        >  2:1260890/1‑150 (MQ=255)
                                                                                                         gcTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGt                                       <  2:408037/151‑1 (MQ=255)
                                                                                                            ataGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTg                                      <  1:851755/149‑1 (MQ=255)
                                                                                                               gTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGAc                                  >  2:1009239/1‑150 (MQ=255)
                                                                                                                   cACCGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGAGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAgg                              >  2:927252/1‑150 (MQ=255)
                                                                                                                      cGACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTc                          >  1:502478/1‑151 (MQ=255)
                                                                                                                       gACATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTct                         <  2:1256323/151‑1 (MQ=255)
                                                                                                                        aCATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGGATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTctg                        <  1:502503/151‑1 (MQ=255)
                                                                                                                        aCATGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTctg                        >  2:1484702/1‑151 (MQ=255)
                                                                                                                          aTGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTctggc                      <  2:84660/151‑1 (MQ=255)
                                                                                                                          aTGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTctggc                      <  2:1228124/151‑1 (MQ=255)
                                                                                                                          aTGGCGGCGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTctgg                       >  1:575121/1‑150 (MQ=255)
                                                                                                                            ggcggcGGTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTctggct                     <  2:822595/150‑1 (MQ=255)
                                                                                                                               ggcggTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCtggctggc                  >  1:82297/1‑150 (MQ=255)
                                                                                                                                   gTCGCTATCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCTGGCTGGCAGGtt             <  1:364490/151‑1 (MQ=255)
                                                                                                                                         aTCGGTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCTGGCTGGCAGGTTTTGttt       >  1:1125589/1‑151 (MQ=255)
                                                                                                                                             gTACTTATATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCTGGCTGGCAGGTTTTGTTTCCGt   >  1:885064/1‑151 (MQ=255)
                                                                                                                                               aCTTATATCTGGCTTCCGGCGATCCTCTTTGGTCGCGGACTGCTGCGGGCATTAGCGCCGGTTATCCCGCAATTAAATGGTTCCGGTCGTGGTGGGCGCATTGCGCATCAGGTGCGACATGGTTTCTGGCTGGCAGGTTTTGTTTCCGtt  <  2:1327055/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGTATATCAGTGAAGCGCGTCTGTTATTAGCATTAGCAATCCCGGTGATTCTCGCGCAAATCGCCCAAACTGCGATGGGTTTTGTCGATACCGTGATGGCGGGCGGCTATAGTGCCACCGACATGGCGGCGGTCGCTATCGGTACTTCTATCTGGCTTCCGGCGATCCTCTTTGGTCACGGACTGCTGCTGGCATTAACGCCGGTTATCGCGCAATTAAATGGTTCCGGTCGACGTGAGCGCATTGCGCATCAGGTGCGACAAGGTTTCTGGCTGGCAGGTTTTGTTTCCGTT  >  NC_000913/1743464‑1743756

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