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

TGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCA  >  NC_000913/2588811‑2589000
                                                                                                    |                                                                                           
tggATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCa                                                                                              <  1:342617/100‑1 (MQ=255)
tggATGACGCCATCGTGGTGGGGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCa                                                                                              >  2:2047568/1‑100 (MQ=255)
 ggATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAg                                                                                            >  2:2635736/1‑100 (MQ=255)
 ggATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAg                                                                                            >  1:1646941/1‑100 (MQ=255)
     gACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAggggg                                                                                       <  1:25760/100‑1 (MQ=255)
      aCGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAgggggg                                                                                      <  1:2249206/100‑2 (MQ=255)
       cGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAggggggc                                                                                     >  1:537361/1‑98 (MQ=255)
       cGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAggggggc                                                                                     >  1:1874568/1‑98 (MQ=255)
         ccATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAggggggcac                                                                                   <  1:2681050/100‑5 (MQ=255)
         ccATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCAATGGGGCAGATCCAggggggcac                                                                                   >  1:214849/1‑96 (MQ=255)
          cATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAggggggcact                                                                                  <  2:104810/100‑6 (MQ=255)
          cATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAggggggcact                                                                                  >  2:2500755/1‑95 (MQ=255)
                 gtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCggg                                                                            <  1:676778/100‑1 (MQ=255)
                 gtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCggg                                                                            >  1:1308588/1‑100 (MQ=255)
                  tggtggAAAACGTCGAACGTATTATGAGTGAGGCAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCggg                                                                            >  1:2286117/1‑99 (MQ=255)
                   ggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGAt                                                                          >  2:414940/1‑100 (MQ=255)
                    gtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGAtt                                                                         >  2:2321916/1‑100 (MQ=255)
                     tggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTg                                                                        <  2:1146240/100‑1 (MQ=255)
                         aaaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGAt                                                                    >  1:833835/1‑100 (MQ=255)
                         aaaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGAt                                                                    <  1:2500770/100‑1 (MQ=255)
                           aCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATgg                                                                  <  1:657344/100‑1 (MQ=255)
                            cGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATgg                                                                  >  2:74321/1‑99 (MQ=255)
                            cGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGt                                                                 >  1:2800609/1‑100 (MQ=255)
                            cGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGt                                                                 >  2:952217/1‑100 (MQ=255)
                              tCGACCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTc                                                               >  1:1409874/1‑100 (MQ=255)
                               cGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCt                                                              >  1:217730/1‑100 (MQ=255)
                                gAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCtt                                                             >  2:391935/1‑100 (MQ=255)
                                 aaCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCttt                                                            <  1:19271/100‑1 (MQ=255)
                                  aCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTc                                                           <  2:2746334/100‑1 (MQ=255)
                                   cGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTcg                                                          <  2:1409866/100‑1 (MQ=255)
                                     tattatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTcggc                                                        <  1:2173397/100‑1 (MQ=255)
                                     tattatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTcggc                                                        <  2:179987/100‑1 (MQ=255)
                                      attatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTcggcg                                                       >  1:2987837/1‑100 (MQ=255)
                                      attatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTcggc                                                        <  2:2295599/99‑1 (MQ=255)
                                        tatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGt                                                     <  2:2476212/100‑1 (MQ=255)
                                         atGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTa                                                    >  1:835261/1‑100 (MQ=255)
                                          tgagtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTAt                                                   >  1:668848/1‑100 (MQ=255)
                                           gagtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTAtt                                                  <  1:2614004/100‑1 (MQ=255)
                                           gagtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCAATGGTTCTTTCGGCGGTAt                                                   >  2:1259452/1‑99 (MQ=255)
                                            agtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTAttt                                                 >  2:2565473/1‑100 (MQ=255)
                                             gtgagGCAGGACTCCCTCCTCGCGAGGCCACACGTAAATCGATGGGGCAGAGCCAGAGGGGCACTGGTCGGGAGTGCGATGGTTCTTGCGGCGGTATTTg                                                <  2:2286105/100‑1 (MQ=255)
                                              tgagGAAGGACTCACTCCTCGCGAAGCCACACTTAAATCGATGTGGCAGATCCATGGGGGCACTGGTCGGGATGGCGATGGTTCTTTCGGCGGTATTTGt                                               >  2:597803/1‑100 (MQ=255)
                                              tgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGt                                               <  1:211907/100‑1 (MQ=255)
                                               gagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTa                                              <  2:251909/100‑1 (MQ=255)
                                               gagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTa                                              <  1:1114091/100‑1 (MQ=255)
                                                 ggaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            >  1:970660/1‑100 (MQ=255)
                                                 ggaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            <  1:2244219/100‑1 (MQ=255)
                                                 ggaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            >  2:2725106/1‑100 (MQ=255)
                                                 ggaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            >  1:1342191/1‑100 (MQ=255)
                                                 ggaaggaCGCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                            <  2:2576822/100‑1 (MQ=255)
                                                    agggCTCACTCCTCGGGGAGCCACACGGAAATCGATGGGGCAGGTCCAGGGGGGCACTGGTCGGGAATGCGATGGTTCTTTCGGCGGTATTTGTACCaaa                                         <  2:1874557/100‑2 (MQ=255)
                                                    aggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAAt                                         >  1:2490654/1‑100 (MQ=255)
                                                      gaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATgg                                       <  1:1841578/100‑1 (MQ=255)
                                                      gaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATgg                                       <  2:1263842/100‑1 (MQ=255)
                                                        ctcactcCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGcc                                     >  1:1930526/1‑100 (MQ=255)
                                                        ctaactcCTCGGGAAGACACACGTAAATCGATGGGGCAGATCCAGGGGGGGACTGTTCGGGATTGCGATGGTTATTTCGGCAGTATTTGTACCAATGGcc                                     <  2:1468193/97‑1 (MQ=255)
                                                              cctcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCttcttc                               <  2:537358/100‑1 (MQ=255)
                                                              cctcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCttcttc                               >  1:2262249/1‑100 (MQ=255)
                                                               ctcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcg                              <  2:200028/100‑1 (MQ=255)
                                                                tcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcgg                             >  1:1548299/1‑100 (MQ=255)
                                                                 cgcgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggc                            >  1:1969295/1‑100 (MQ=255)
                                                                 cgcgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggc                            <  2:605253/100‑1 (MQ=255)
                                                                 cgcgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggcg                           >  2:48399/1‑100 (MQ=255)
                                                                  gcgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggcg                           <  1:1437026/100‑1 (MQ=255)
                                                                   cgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggcgg                          <  1:2913104/100‑1 (MQ=255)
                                                                     aaGCCACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGca                        <  1:414942/100‑1 (MQ=255)
                                                                        ccACACGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGcacca                     <  1:1163197/100‑1 (MQ=255)
                                                                         cacaCGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGcaccac                    >  1:2184490/1‑100 (MQ=255)
                                                                         cacaCGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGACTTCTTCGGCGGcacaa                     >  2:1451914/1‑97 (MQ=255)
                                                                           cacGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCg                  >  1:1504535/1‑100 (MQ=255)
                                                                             cGTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGt                >  2:145879/1‑100 (MQ=255)
                                                                              gTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTg               <  1:2565487/100‑1 (MQ=255)
                                                                              gTAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTg               <  2:888087/100‑1 (MQ=255)
                                                                               tAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGc              >  1:2513393/1‑100 (MQ=255)
                                                                               tAAATCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGc              <  1:675281/100‑1 (MQ=255)
                                                                                 aaTCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCa            >  2:1075075/1‑100 (MQ=255)
                                                                                  aTCGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCAt           <  2:1308581/100‑1 (MQ=255)
                                                                                    cGATGGGGCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCt         >  2:1864903/1‑100 (MQ=255)
                                                                                     gaggtcgCAGAGCCAGTGGTGCACTGGTCGTGATTTCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGTGGCACCAGTGGTGAAATCTa        <  2:2789934/94‑1 (MQ=255)
                                                                                        ggggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCg     <  2:214847/100‑1 (MQ=255)
                                                                                        ggggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCg     <  1:1549997/100‑1 (MQ=255)
                                                                                        ggggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCg     >  2:1050593/1‑100 (MQ=255)
                                                                                        ggggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCg     >  1:209657/1‑100 (MQ=255)
                                                                                         nggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGc    >  1:2707222/2‑100 (MQ=255)
                                                                                         gggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGc    >  2:1154297/1‑100 (MQ=255)
                                                                                         gggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGc    >  1:727194/1‑100 (MQ=255)
                                                                                         gggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGc    <  1:391937/100‑1 (MQ=255)
                                                                                          ggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGcc   >  1:2323902/1‑100 (MQ=255)
                                                                                          ggCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGcc   >  1:2323832/1‑100 (MQ=255)
                                                                                           gCAGATCCAGGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCa  <  2:2947625/100‑1 (MQ=255)
                                                                                                    |                                                                                           
TGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCA  >  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: