Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,619,743 A→C Y57S (TAC→TCC)  marA → multiple antibiotic resistance transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,619,7430AC97.8% 448.6 / ‑0.4 137Y57S (TAC→TCC) marAmultiple antibiotic resistance transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (80/53);  total (83/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 142 total aligned reads displayed.

CGCCACTGTCACTGGAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATACATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGTCGCA  >  NC_000913/1619644‑1619839
                                                                                                   |                                                                                                
cGCCACTGTCACTGGAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGCTGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATa                                                                                                  >  2:864020/1‑100 (MQ=255)
  ccACTGTCACTGGAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCa                                                                                                <  2:1480442/100‑1 (MQ=255)
   cACTGTCACTGGAGAAAGTGTCAGAGCATTCGGGTTACTCCAAATGGCACCTGCAACGGGTGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCAt                                                                                               >  1:3602378/1‑100 (MQ=255)
    aCTGTCACTGGAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATc                                                                                              <  2:2667151/100‑1 (MQ=255)
     cTGTCACTGGAGAAGGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAatccatcc                                                                                             <  1:3313090/100‑1 (MQ=255)
     cTGTCACTGGAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAATCGGTCATTCATTAGGCCAatccatcc                                                                                             <  2:2856699/100‑1 (MQ=255)
       gTCACTGGAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGc                                                                                           >  1:171598/1‑100 (MQ=255)
       gTCACTGGAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGc                                                                                           >  1:2492935/1‑100 (MQ=255)
           cTGGAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGGTTAAAAAAGACACCGGTCATTCATTATGCCAATCCATCCTCacca                                                                                       >  2:3155478/1‑97 (MQ=255)
             ggAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGt                                                                                     >  1:2059189/1‑100 (MQ=255)
              gagaAAGTAACAAAAAGTACGAGTTACTCCAAATGTCACCTGTAAAAAAAAAAAAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTa                                                                                    <  2:2889433/100‑1 (MQ=255)
                 aaaGTGTCAGAGCGTGTGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGCCCAATCCATCCGCAGCCGTAAGa                                                                                 <  2:3352114/100‑1 (MQ=255)
                  aaGCGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGCCATTCATTAGGCCAATCCATCCGCAGCCGTAACAt                                                                                >  2:3671100/1‑100 (MQ=255)
                    gtgtCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGa                                                                              >  1:1668108/1‑100 (MQ=255)
                    gtgtCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGa                                                                              >  1:448552/1‑100 (MQ=255)
                     tgtCAGAGCGTTCCGGTTACTCCAAATGCCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGAc                                                                             >  1:3590648/1‑100 (MQ=255)
                      gtCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACg                                                                            >  1:2531152/1‑100 (MQ=255)
                       tCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACgg                                                                           <  1:645903/100‑1 (MQ=255)
                       tCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACgg                                                                           <  2:2007063/100‑1 (MQ=255)
                           agCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAAt                                                                       <  2:866779/100‑1 (MQ=255)
                           agCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAAt                                                                       <  1:2964432/100‑1 (MQ=255)
                             cGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATcg                                                                     <  2:3883876/100‑1 (MQ=255)
                             cGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATcg                                                                     <  2:3883830/100‑1 (MQ=255)
                             cGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATACATCCGCAGCCGTAAGATGACGGAAATcg                                                                     <  2:1477338/100‑1 (MQ=255)
                              gTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATcgc                                                                    >  1:3420984/1‑100 (MQ=255)
                                tCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCgcgc                                                                  >  2:3221329/1‑100 (MQ=255)
                                tCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCgcgc                                                                  <  2:3482780/100‑1 (MQ=255)
                                 cGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCa                                                                 >  2:2258171/1‑100 (MQ=255)
                                    gTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCaaaa                                                              <  1:3375187/100‑1 (MQ=255)
                                     ttACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAg                                                             <  2:1124466/100‑1 (MQ=255)
                                      tACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGTCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGc                                                            <  1:2004657/100‑1 (MQ=255)
                                      tACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGc                                                            <  1:3197298/100‑1 (MQ=255)
                                       aCTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATTCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCt                                                           >  1:3725548/1‑100 (MQ=255)
                                        cTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTTATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTg                                                          >  2:1329038/1‑100 (MQ=255)
                                         tCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGa                                                         >  2:1951526/1‑100 (MQ=255)
                                          ccAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGaa                                                        >  2:3366342/1‑100 (MQ=255)
                                          ccAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGaa                                                        >  2:1264231/1‑100 (MQ=255)
                                          ccAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGaa                                                        >  2:633786/1‑100 (MQ=255)
                                          ccAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGaa                                                        >  2:120649/1‑100 (MQ=255)
                                          ccAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGaa                                                        >  1:1852353/1‑100 (MQ=255)
                                          ccAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTAATTCATTAGGCTAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGaa                                                        <  2:3973873/100‑1 (MQ=255)
                                           cAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAg                                                       >  1:1592489/1‑100 (MQ=255)
                                           cAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAg                                                       <  1:1006157/100‑1 (MQ=255)
                                           aacaTGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCTTTAGGCCAATCCTTCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAg                                                       <  1:1610324/97‑1 (MQ=255)
                                            aaaTGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAgg                                                      >  1:976461/1‑100 (MQ=255)
                                              aTGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGaa                                                    >  1:1373055/1‑100 (MQ=255)
                                                ggCACCTGCCACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAg                                                  <  2:3457966/100‑1 (MQ=255)
                                                 gCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGt                                                 <  2:245246/100‑1 (MQ=255)
                                                   aCCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTaa                                               <  2:371560/100‑1 (MQ=255)
                                                    ccTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAAc                                              <  2:3491632/100‑1 (MQ=255)
                                                    ccTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGAAGCCGTAAGATAACGGAAATCACGCAAAAGCTGAAGGAAAGTAAc                                              >  2:1419092/1‑100 (MQ=255)
                                                    ccTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCAGCCGCAGCCGTAAGCTGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAAc                                              >  2:3943639/1‑100 (MQ=255)
                                                    ccTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCACTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAAc                                              <  2:3630500/100‑1 (MQ=255)
                                                      tGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGGCGTAAGATGACGGAAATCGCGCAAAACCTGAAGGAAGGTAACGa                                            >  1:3754062/1‑100 (MQ=255)
                                                      tGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGa                                            >  2:2093115/1‑100 (MQ=255)
                                                      tGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGa                                            >  1:294278/1‑100 (MQ=255)
                                                      tGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGa                                            >  1:1537237/1‑100 (MQ=255)
                                                       gCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAg                                           <  2:292405/100‑1 (MQ=255)
                                                         aaCGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGcc                                         <  1:2860254/100‑1 (MQ=255)
                                                         aaCGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGcc                                         <  1:1823298/100‑1 (MQ=255)
                                                         aaCGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGcc                                         <  1:1784055/100‑1 (MQ=255)
                                                           cGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGa                                       >  1:3939786/1‑100 (MQ=255)
                                                            ggATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGAt                                      <  1:1613347/100‑1 (MQ=255)
                                                                gTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGCAAGTAACGAGCCGATActc                                  >  1:392617/1‑100 (MQ=255)
                                                                gTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATActc                                  <  2:1163006/100‑1 (MQ=255)
                                                                gTTTAAAAAAGAAACCGGTCATTCATTAGGACAATCCATCAGAAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGACGGAAAGTAACGAGCCGATActc                                  >  1:520781/1‑100 (MQ=255)
                                                                 tttAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATActct                                 >  1:3165751/1‑100 (MQ=255)
                                                                 tttAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATActct                                 >  2:551398/1‑100 (MQ=255)
                                                                  ttAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTa                                >  1:309833/1‑100 (MQ=255)
                                                                  ttAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTa                                >  1:310034/1‑100 (MQ=255)
                                                                  ttAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTa                                >  2:1023206/1‑100 (MQ=255)
                                                                  ttAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATAATCTa                                >  1:3288786/1‑100 (MQ=255)
                                                                   tAAAAAAGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTAt                               >  1:3564786/1‑100 (MQ=255)
                                                                       aaaGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTgg                           >  2:1833578/1‑100 (MQ=255)
                                                                        aaGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGc                          <  1:887978/100‑1 (MQ=255)
                                                                        aaGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGATGAAGGAAAGTAACGAGCCGATACTCTATCTGGc                          >  2:2344363/1‑100 (MQ=255)
                                                                         aGAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCa                         <  2:2987358/100‑1 (MQ=255)
                                                                          gAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAg                        <  1:1506612/100‑1 (MQ=255)
                                                                          gAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAg                        <  2:2436267/100‑1 (MQ=255)
                                                                          gAAACCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAg                        >  2:508749/1‑100 (MQ=255)
                                                                             aCCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAAc                     >  2:460773/1‑100 (MQ=255)
                                                                             aCCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAAc                     <  1:4026816/100‑1 (MQ=255)
                                                                             aCCGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAAc                     <  1:710663/100‑1 (MQ=255)
                                                                              ccGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACg                    >  1:2425566/1‑100 (MQ=255)
                                                                               cGGTCATTCATTAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGa                   <  1:784113/100‑1 (MQ=255)
                                                                                       cattAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCtt           >  1:32749/1‑100 (MQ=255)
                                                                                        attAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTc          >  1:1735128/1‑100 (MQ=255)
                                                                                        attAGGCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTc          >  1:2803388/1‑100 (MQ=255)
                                                                                           aGGCCAATACATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAg       >  1:1034302/1‑100 (MQ=255)
                                                                                            ggCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGt      >  2:2062150/1‑100 (MQ=255)
                                                                                            ggCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGt      >  2:3698762/1‑100 (MQ=255)
                                                                                            ggCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGt      >  1:2257851/1‑100 (MQ=255)
                                                                                            ggCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGt      <  1:1089717/100‑1 (MQ=255)
                                                                                            ggCCAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCAATCTGGCAGAACGATATGGCTTCGAGt      >  2:1678549/1‑100 (MQ=255)
                                                                                              cgaaTCCATCCGCGGCAGTAAGATGACGGAAATCGAGCAAAAGCTGAAGGAAGGTAACGGGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGTCg    <  2:2302552/98‑1 (MQ=255)
                                                                                              ccAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGTCg    >  2:1783428/1‑100 (MQ=255)
                                                                                              ccAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGTCg    >  1:1206374/1‑100 (MQ=255)
                                                                                              ccAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGTCg    >  1:1669568/1‑100 (MQ=255)
                                                                                               cAATCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCGTCGGGTCGc   >  2:665600/1‑100 (MQ=255)
                                                                                                aaTCCATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTCACGAGCCGATACTGTATCTGGCAGAACGATATGGCTTCGAGTAGca  >  2:3632674/1‑100 (MQ=255)
                                                                                                   |                                                                                                
CGCCACTGTCACTGGAGAAAGTGTCAGAGCGTTCGGGTTACTCCAAATGGCACCTGCAACGGATGTTTAAAAAAGAAACCGGTCATTCATTAGGCCAATACATCCGCAGCCGTAAGATGACGGAAATCGCGCAAAAGCTGAAGGAAAGTAACGAGCCGATACTCTATCTGGCAGAACGATATGGCTTCGAGTCGCA  >  NC_000913/1619644‑1619839

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