Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,588,915 (G)5→6 coding (1321/3114 nt) acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,588,9101.G96.1% 273.2 / 4.1 78coding (1316/3114 nt)acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (38/36);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑G‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCA  >  NC_000913/2588811‑2589000
                                                                                                    |                                                                                           
tggATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCa                                                                                              <  1:342617/100‑1 (MQ=255)
tggATGACGCCATCGTGGTGGGGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCa                                                                                              >  2:2047568/1‑100 (MQ=255)
 ggATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑g                                                                                            >  2:2635736/1‑100 (MQ=255)
 ggATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑g                                                                                            >  1:1646941/1‑100 (MQ=255)
     gACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑g‑gggg                                                                                       <  1:25760/100‑1 (MQ=255)
      aCGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑g‑ggggg                                                                                      <  1:2249206/100‑2 (MQ=255)
       cGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑g‑gggggc                                                                                     >  1:537361/1‑98 (MQ=255)
       cGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑g‑gggggc                                                                                     >  1:1874568/1‑98 (MQ=255)
         ccATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑g‑gggggcac                                                                                   <  1:2681050/100‑5 (MQ=255)
         ccATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCAATGGGGCAGATCCA‑g‑gggggcac                                                                                   >  1:214849/1‑96 (MQ=255)
          cATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑g‑gggggcact                                                                                  <  2:104810/100‑6 (MQ=255)
          cATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑g‑gggggcact                                                                                  >  2:2500755/1‑95 (MQ=255)
                 gtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCggg                                                                            <  1:676778/100‑1 (MQ=255)
                 gtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCggg                                                                            >  1:1308588/1‑100 (MQ=255)
                  tggtggAAAACGTCGAACGTATTATGAGTGAGGCAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCggg                                                                            >  1:2286117/1‑99 (MQ=255)
                   ggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGAt                                                                          >  2:414940/1‑100 (MQ=255)
                    gtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGAtt                                                                         >  2:2321916/1‑100 (MQ=255)
                     tggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTg                                                                        <  2:1146240/100‑1 (MQ=255)
                         aaaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGAt                                                                    >  1:833835/1‑100 (MQ=255)
                         aaaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGAt                                                                    <  1:2500770/100‑1 (MQ=255)
                           aCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATgg                                                                  <  1:657344/100‑1 (MQ=255)
                            cGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATgg                                                                  >  2:74321/1‑99 (MQ=255)
                            cGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGt                                                                 >  1:2800609/1‑100 (MQ=255)
                            cGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGt                                                                 >  2:952217/1‑100 (MQ=255)
                              tCGACCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTc                                                               >  1:1409874/1‑100 (MQ=255)
                               cGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCt                                                              >  1:217730/1‑100 (MQ=255)
                                gAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCtt                                                             >  2:391935/1‑100 (MQ=255)
                                 aaCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCttt                                                            <  1:19271/100‑1 (MQ=255)
                                  aCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTc                                                           <  2:2746334/100‑1 (MQ=255)
                                   cGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTcg                                                          <  2:1409866/100‑1 (MQ=255)
                                     tattatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTcggc                                                        <  1:2173397/100‑1 (MQ=255)
                                     tattatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTcggc                                                        <  2:179987/100‑1 (MQ=255)
                                      attatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTcggcg                                                       >  1:2987837/1‑100 (MQ=255)
                                      attatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTcggc                                                        <  2:2295599/99‑1 (MQ=255)
                                        tatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGt                                                     <  2:2476212/100‑1 (MQ=255)
                                         atGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTa                                                    >  1:835261/1‑100 (MQ=255)
                                          tgagtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTAt                                                   >  1:668848/1‑100 (MQ=255)
                                           gagtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTAtt                                                  <  1:2614004/100‑1 (MQ=255)
                                           gagtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCAATGGTTCTTTCGGCGGTAt                                                   >  2:1259452/1‑99 (MQ=255)
                                            agtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTAttt                                                 >  2:2565473/1‑100 (MQ=255)
                                             gtgagGCAGGACTCCCTCCTCGCGAGGCCACACGTAAATCGATGGGGCAGAGCCA‑GAGGGGCACTGGTCGGGAGTGCGATGGTTCTTGCGGCGGTATTTg                                                <  2:2286105/100‑1 (MQ=255)
                                              tgagGAAGGACTCACTCCTCGCGAAGCCACACTTAAATCGATGTGGCAGATCCATG‑GGGGCACTGGTCGGGATGGCGATGGTTCTTTCGGCGGTATTTGt                                               >  2:597803/1‑100 (MQ=255)
                                              tgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGt                                               <  1:211907/100‑1 (MQ=255)
                                               gagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTa                                              <  2:251909/100‑1 (MQ=255)
                                               gagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTa                                              <  1:1114091/100‑1 (MQ=255)
                                                 ggaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            >  1:970660/1‑100 (MQ=255)
                                                 ggaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            <  1:2244219/100‑1 (MQ=255)
                                                 ggaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            >  2:2725106/1‑100 (MQ=255)
                                                 ggaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            >  1:1342191/1‑100 (MQ=255)
                                                 ggaaggaCGCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            <  2:2576822/100‑1 (MQ=255)
                                                    agggCTCACTCCTCGGGGAGCCACACGGAAATCGATGGGGCAGGTCCAGG‑GGGGCACTGGTCGGGAATGCGATGGTTCTTTCGGCGGTATTTGTACCaaa                                         <  2:1874557/100‑2 (MQ=255)
                                                    aggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAAt                                         >  1:2490654/1‑100 (MQ=255)
                                                      gaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATgg                                       <  1:1841578/100‑1 (MQ=255)
                                                      gaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATgg                                       <  2:1263842/100‑1 (MQ=255)
                                                        ctcactcCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGcc                                     >  1:1930526/1‑100 (MQ=255)
                                                        ctaactcCTCGGGAAGACACACGTAAATCGATGGGGCAGATCCAGG‑GGGGGACTGTTCGGGATTGCGATGGTTATTTCGGCAGTATTTGTACCAATGGcc                                     <  2:1468193/97‑1 (MQ=255)
                                                              cctcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCttcttc                               <  2:537358/100‑1 (MQ=255)
                                                              cctcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCttcttc                               >  1:2262249/1‑100 (MQ=255)
                                                               ctcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcg                              <  2:200028/100‑1 (MQ=255)
                                                                tcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcgg                             >  1:1548299/1‑100 (MQ=255)
                                                                 cgcgAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggc                            >  1:1969295/1‑100 (MQ=255)
                                                                 cgcgAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggc                            <  2:605253/100‑1 (MQ=255)
                                                                 cgcgAAGCCACACGTAAATCGATGGGGCAGATCCA‑G‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggcg                           >  2:48399/1‑100 (MQ=255)
                                                                  gcgAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggcg                           <  1:1437026/100‑1 (MQ=255)
                                                                   cgAAGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggcgg                          <  1:2913104/100‑1 (MQ=255)
                                                                     aaGCCACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGca                        <  1:414942/100‑1 (MQ=255)
                                                                        ccACACGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGcacca                     <  1:1163197/100‑1 (MQ=255)
                                                                         cacaCGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGcaccac                    >  1:2184490/1‑100 (MQ=255)
                                                                         cacaCGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGACTTCTTCGGCGGcacaa                     >  2:1451914/1‑97 (MQ=255)
                                                                           cacGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCg                  >  1:1504535/1‑100 (MQ=255)
                                                                             cGTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGt                >  2:145879/1‑100 (MQ=255)
                                                                              gTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTg               <  1:2565487/100‑1 (MQ=255)
                                                                              gTAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTg               <  2:888087/100‑1 (MQ=255)
                                                                               tAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGc              >  1:2513393/1‑100 (MQ=255)
                                                                               tAAATCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGc              <  1:675281/100‑1 (MQ=255)
                                                                                 aaTCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCa            >  2:1075075/1‑100 (MQ=255)
                                                                                  aTCGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCAt           <  2:1308581/100‑1 (MQ=255)
                                                                                    cGATGGGGCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCt         >  2:1864903/1‑100 (MQ=255)
                                                                                     gaggtcgCAGAGCCA‑GTGGTGCACTGGTCGTGATTTCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGTGGCACCAGTGGTGAAATCTa        <  2:2789934/94‑1 (MQ=255)
                                                                                        ggggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCg     <  2:214847/100‑1 (MQ=255)
                                                                                        ggggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCg     <  1:1549997/100‑1 (MQ=255)
                                                                                        ggggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCg     >  2:1050593/1‑100 (MQ=255)
                                                                                        ggggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCg     >  1:209657/1‑100 (MQ=255)
                                                                                         nggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGc    >  1:2707222/2‑100 (MQ=255)
                                                                                         gggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGc    >  2:1154297/1‑100 (MQ=255)
                                                                                         gggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGc    >  1:727194/1‑100 (MQ=255)
                                                                                         gggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGc    <  1:391937/100‑1 (MQ=255)
                                                                                          ggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGcc   >  1:2323902/1‑100 (MQ=255)
                                                                                          ggCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGcc   >  1:2323832/1‑100 (MQ=255)
                                                                                           gCAGATCCAGG‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCa  <  2:2947625/100‑1 (MQ=255)
                                                                                                    |                                                                                           
TGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCA‑G‑GGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCA  >  NC_000913/2588811‑2589000

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