Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 905,720 C→A 100% L270M (CTG→ATG)  amiD → 1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase, Zn‑dependent; OM lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913905,7200CA98.3% 388.2 / ‑2.6 120L270M (CTG→ATG) amiD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase, Zn‑dependent; OM lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (52/66);  minor base T (2/0);  total (54/66)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 120 total aligned reads displayed.

GCAGCGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCT  >  NC_000913/905623‑905816
                                                                                                 |                                                                                                
gcagcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATg                                                                                                <  2:1407684/100‑1 (MQ=255)
 ctgcggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATgg                                                                                               <  2:2895453/98‑1 (MQ=255)
 cagcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATgg                                                                                               <  1:1773529/100‑1 (MQ=255)
  agcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGTCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATgg                                                                                               <  1:2676062/99‑1 (MQ=255)
  agcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGga                                                                                              <  1:1151922/100‑1 (MQ=255)
  agcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGga                                                                                              >  2:2356488/1‑100 (MQ=255)
   gcggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGgag                                                                                             <  2:265643/100‑1 (MQ=255)
    cggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGgaga                                                                                            <  2:1550898/100‑1 (MQ=255)
     ggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGaa                                                                                           >  1:848316/1‑100 (MQ=255)
     ggcgCGTGATTATGGCATTCCAGATGCATTTACGCCCGACGTTATATAAAGGCGAAGCGGATGCATAAACTCAGGCGATTGGCGAAGCATTGATGGAGaa                                                                                           >  2:523617/1‑100 (MQ=255)
      gcgcGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGaaa                                                                                          >  2:2160221/1‑100 (MQ=255)
      gcgcGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGaaa                                                                                          >  2:1859234/1‑100 (MQ=255)
      gcgcGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGaaa                                                                                          <  2:1202854/100‑1 (MQ=255)
        gcgTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATa                                                                                        <  1:3030679/100‑1 (MQ=255)
          gTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATAc                                                                                       >  2:1395239/1‑99 (MQ=255)
            gATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACggg                                                                                    >  2:1220102/1‑100 (MQ=255)
            gATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACggg                                                                                    >  2:2951417/1‑100 (MQ=255)
              ttATGGCATTCCAGATGTATTTCCGATCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCa                                                                                  <  1:577497/100‑1 (MQ=255)
               tATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAg                                                                                 >  1:48039/1‑100 (MQ=255)
                aTGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAgg                                                                                >  1:2401768/1‑100 (MQ=255)
                 tGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGa                                                                               >  1:144439/1‑100 (MQ=255)
                 tGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGa                                                                               >  1:2470451/1‑100 (MQ=255)
                 tGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGa                                                                               <  1:2003004/100‑1 (MQ=255)
                   gCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGAtt                                                                             <  2:1318788/100‑1 (MQ=255)
                      ttCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgc                                                                          <  2:2378310/100‑1 (MQ=255)
                       tCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgcg                                                                         >  2:1192260/1‑100 (MQ=255)
                        ccAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgcgc                                                                        <  2:2681639/100‑1 (MQ=255)
                        ccAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgcgc                                                                        <  1:437384/100‑1 (MQ=255)
                        ccAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACGCAGGCGATTGCCGAAGCATTAATGGAGAAATAAGGGCAGGATTAgcg                                                                         >  1:1925371/1‑99 (MQ=255)
                           gATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGc                                                                     <  2:1932227/100‑1 (MQ=255)
                            aTGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCa                                                                    <  2:923509/100‑1 (MQ=255)
                                ntttCCGCCCGACGTTATATAACGGCGAAGCGGATGAAGAAACACAGGCGATTGCCGAAGCATTGATGGAGAAAGACGGGCGGGATTAGCGCGGCAGtt                                                                 >  2:728824/2‑99 (MQ=255)
                                 tttCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGtttt                                                               >  1:1415027/1‑100 (MQ=255)
                                 nttCCGCCCGACGTTATATAACGGCGAAGCGGCTGCAGAAACTCAGGCGATTGCCGAAGCATTGAAGGAGAAATACGGGCAGGATTAGCGCGGCAGttt                                                                >  1:1231285/2‑99 (MQ=255)
                                   tCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTcc                                                             >  1:2609071/1‑100 (MQ=255)
                                    ccGCCCGACGTTATATATCGGCGAAGCGGATGCAGAAACTCAGGGGAATGCCGAAGCATTGATGGGGAAAAAAGGGCAGGAGTAGCGCGGCAGTTTTc                                                              <  1:2549250/98‑1 (MQ=255)
                                    ccGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCg                                                            <  1:1348231/100‑1 (MQ=255)
                                    ccGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCg                                                            >  2:2481260/1‑100 (MQ=255)
                                     cGCCCGACGTTATATACCGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGTGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCg                                                            <  2:1639359/99‑1 (MQ=255)
                                         cGACGTTATATAACGTCGAAGCGGAAGCAGAAACTCAGGCGATTGCCTAAGCATTGATGGAGAACTAAGGGCAGGATTAGCGCGGCAGTTTTCCg                                                            >  2:443087/1‑95 (MQ=255)
                                          gACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTcg                                                      >  2:268310/1‑100 (MQ=255)
                                            cGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGATAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTcgcg                                                    <  2:1405067/100‑1 (MQ=255)
                                            cGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTcgcg                                                    <  2:2341794/100‑1 (MQ=255)
                                             gTTATATAACGGCGAAGCGGATGCAGAATCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGGAGTTTTCCGTGGTCGCGt                                                   <  1:2550732/100‑1 (MQ=255)
                                             gTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGt                                                   <  1:2026931/100‑1 (MQ=255)
                                              ttATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTa                                                  >  2:1581008/1‑100 (MQ=255)
                                                atataACGGCGACGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGc                                                >  1:1990365/1‑100 (MQ=255)
                                                atataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGc                                                <  1:2431902/100‑1 (MQ=255)
                                                  ataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCa                                              >  2:1057867/1‑100 (MQ=255)
                                                   taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAg                                             >  2:2308860/1‑100 (MQ=255)
                                                   taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAg                                             <  1:1232275/100‑1 (MQ=255)
                                                   taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAg                                             >  2:2310433/1‑100 (MQ=255)
                                                   taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAg                                             >  1:2191545/1‑100 (MQ=255)
                                                    aaCGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAgg                                            >  1:1687132/1‑100 (MQ=255)
                                                     aCGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAagc                                           >  1:1031447/1‑100 (MQ=255)
                                                       ggCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAggcgg                                         <  1:2707710/100‑1 (MQ=255)
                                                        gCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAggcggc                                        <  2:922941/100‑1 (MQ=255)
                                                         cGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCa                                       >  2:1029659/1‑100 (MQ=255)
                                                         cGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCa                                       >  2:1651020/1‑100 (MQ=255)
                                                          gAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAg                                      <  1:173750/100‑1 (MQ=255)
                                                          gAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAg                                      <  1:173488/100‑1 (MQ=255)
                                                          gAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAg                                      <  2:426911/100‑1 (MQ=255)
                                                               ggATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtc                                 <  2:504320/100‑1 (MQ=255)
                                                                 aTGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtctc                               <  2:1667946/100‑1 (MQ=255)
                                                                 aTGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATCAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtctc                               <  2:1676427/100‑1 (MQ=255)
                                                                 aTGCAAAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACTGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCTCGTAGCCAGGCGGCAGTTTtctc                               <  2:349697/100‑1 (MQ=255)
                                                                  tGCAGAAACTCAGGCGATTGCCGTAGCAGTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtctc                               <  1:2647965/99‑1 (MQ=255)
                                                                  tGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCg                              >  2:2124905/1‑100 (MQ=255)
                                                                  tGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCg                              <  1:2625324/100‑1 (MQ=255)
                                                                    cAGAAACTCCGGCGATTGCCGCAGCATTGATCGCGGAATACTGGCAGGATTAGCGCGCCAGTTTTCCGTCGTCGCGTAGcc                                               <  2:3047692/81‑1 (MQ=255)
                                                                    cAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                            >  2:913588/1‑100 (MQ=255)
                                                                    cAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                            <  1:485117/100‑1 (MQ=255)
                                                                    cAGAAACTCAGGCGATTGACGAAGCATTGATGGAGAAAGACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                            <  2:1420789/100‑1 (MQ=255)
                                                                      gAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATAc                          >  2:2061723/1‑100 (MQ=255)
                                                                      gAAACTCAGGCGATCGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATa                           <  2:684423/99‑1 (MQ=255)
                                                                       aaaCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATAcc                         <  2:3037197/100‑1 (MQ=255)
                                                                       aaaCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATAcc                         <  1:1948792/100‑1 (MQ=255)
                                                                        aaCTTACGCGATTGCCGAAGCATTGATGGGTAAAGACGTGCAGGTTTAGTGCGGCGGTTTTCCGTCGTCGCGTAGCTTGGCGGAAGTTTTCTCGATAcc                         <  2:1214904/99‑1 (MQ=255)
                                                                        aaCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCt                        <  1:2308886/100‑1 (MQ=255)
                                                                        aaCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCt                        <  1:2310459/100‑1 (MQ=255)
                                                                            cAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCAt                    >  2:953726/1‑100 (MQ=255)
                                                                            cAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCAt                    >  2:947884/1‑100 (MQ=255)
                                                                            cAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCAt                    >  2:1435989/1‑100 (MQ=255)
                                                                              ggCGATTGCTGAAGCATTGATTGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATc                   <  1:1412016/99‑1 (MQ=255)
                                                                              ggCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATcc                  >  2:2013830/1‑100 (MQ=255)
                                                                                cGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAg                <  1:1826055/100‑1 (MQ=255)
                                                                                  aTTGCCGAAGCATTGATGGAGAAATACGGTCAGGATTAGCGCGGCAGTTTTGCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAg                <  1:603300/98‑1 (MQ=255)
                                                                                  aTTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAggg              <  1:1478402/100‑1 (MQ=255)
                                                                                     gCCGAAGCATTGTTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtga           >  2:2642741/1‑100 (MQ=255)
                                                                                     gCCGAAGCATTGTTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtg            >  2:614498/1‑99 (MQ=255)
                                                                                     gCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtga           >  2:2811423/1‑100 (MQ=255)
                                                                                     ctcGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtga           <  1:620413/98‑1 (MQ=255)
                                                                                      ccGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgat          >  2:1502389/1‑100 (MQ=255)
                                                                                       cGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatg         <  2:1990789/100‑1 (MQ=255)
                                                                                       cGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgaag         >  2:620409/1‑98 (MQ=255)
                                                                                           gCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGAcc      <  1:1651039/99‑1 (MQ=255)
                                                                                           gCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCg     <  2:1743131/100‑1 (MQ=255)
                                                                                            cATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCgg    >  1:2354797/1‑100 (MQ=255)
                                                                                              ttGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCt  >  2:666899/1‑100 (MQ=255)
                                                                                              ttGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCt  >  1:2142493/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GCAGCGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCT  >  NC_000913/905623‑905816

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