Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 905,720 C→A 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,7200CA96.4% 281.3 / 1.0 85L270M (CTG→ATG) amiD1,6‑anhydro‑N‑acetylmuramyl‑L‑alanine amidase, Zn‑dependent; OM lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (43/38);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CAGCGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCT  >  NC_000913/905624‑905816
                                                                                                |                                                                                                
tccgatctcgTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATgg                                                                                               <  2:311243/92‑1 (MQ=255)
cagcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACTTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATgg                                                                                               <  1:1772980/100‑1 (MQ=255)
cagcGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATgg                                                                                               <  1:1244134/100‑1 (MQ=255)
    ggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGaa                                                                                           <  2:1656200/100‑1 (MQ=255)
    ggcgCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGaa                                                                                           >  1:1754900/1‑100 (MQ=255)
       gcgTGATTATGGCATTCCAGATGCATTTCCACCAGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGGTGGAGAAATa                                                                                        >  1:1276314/1‑100 (MQ=255)
        cgTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGgagatcgga                                                                                       >  1:311241/1‑95 (MQ=255)
        cgTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATAc                                                                                       >  2:1471774/1‑100 (MQ=255)
         gTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACg                                                                                      >  2:1470854/1‑100 (MQ=255)
           gATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACggg                                                                                    <  1:1411386/100‑1 (MQ=255)
              tATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAg                                                                                 >  2:1272995/1‑100 (MQ=255)
               aTGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAgg                                                                                >  2:1684416/1‑100 (MQ=255)
                 ggCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGAt                                                                              >  2:805171/1‑100 (MQ=255)
                   cATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTa                                                                            >  2:1611666/1‑100 (MQ=255)
                       ccAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAgcgc                                                                        <  1:62505/100‑1 (MQ=255)
                         aGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGcgcg                                                                       <  1:9144/99‑1 (MQ=255)
                          gATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGc                                                                     <  1:1611656/100‑1 (MQ=255)
                           aTGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCa                                                                    >  2:590558/1‑100 (MQ=255)
                             gCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGt                                                                  >  1:1320949/1‑100 (MQ=255)
                              cATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGCAACTCAGGCGATTCCCGAAGCATTGATGGAGAAATACGGGCCGGATTAGCGCGGCAGtt                                                                 >  1:693485/1‑100 (MQ=255)
                                tttCCGCCCGACGTTATTTAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGttt                                                                <  2:283809/99‑1 (MQ=255)
                                 ttCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTc                                                              <  2:305483/100‑1 (MQ=255)
                                  tCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTcc                                                             >  2:863488/1‑100 (MQ=255)
                                   ccGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCg                                                            <  1:717831/100‑1 (MQ=255)
                                   ccGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCg                                                            <  1:1177311/100‑1 (MQ=255)
                                    cGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGt                                                           >  2:1004024/1‑100 (MQ=255)
                                     gCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTg                                                          >  1:1226814/1‑100 (MQ=255)
                                       ccGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGt                                                        <  2:690056/100‑1 (MQ=255)
                                        cGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTc                                                       <  2:1809431/100‑1 (MQ=255)
                                        cGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTc                                                       <  2:1809730/100‑1 (MQ=255)
                                        cGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTc                                                       >  1:547265/1‑100 (MQ=255)
                                             ttATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTa                                                  <  1:1489328/100‑1 (MQ=255)
                                              tagataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAg                                                 <  1:532584/97‑1 (MQ=255)
                                                tataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGcc                                               >  2:984309/1‑100 (MQ=255)
                                                tataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGcc                                               >  2:431805/1‑100 (MQ=255)
                                                 ataACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGcc                                               >  2:1189036/1‑99 (MQ=255)
                                                  taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAg                                             <  1:321333/100‑1 (MQ=255)
                                                  taACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAg                                             >  1:1733911/1‑100 (MQ=255)
                                                   aaCGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAgg                                            >  1:630701/1‑100 (MQ=255)
                                                   aaCGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAgg                                            >  1:262999/1‑100 (MQ=255)
                                                   aaCGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAgg                                            >  1:1237842/1‑100 (MQ=255)
                                                     cGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAggcg                                          >  2:543424/1‑100 (MQ=255)
                                                     cGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAggcg                                          <  2:560972/100‑1 (MQ=255)
                                                      ggCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAggcgg                                         >  2:622193/1‑100 (MQ=255)
                                                        cGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCa                                       <  2:1237846/100‑1 (MQ=255)
                                                         gAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAg                                      <  1:1559113/100‑1 (MQ=255)
                                                           aGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGtt                                    >  2:562/1‑100 (MQ=255)
                                                           aGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGtt                                    <  2:833517/100‑1 (MQ=255)
                                                           aGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATCAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGtt                                    >  1:1147547/1‑100 (MQ=255)
                                                            gCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGttt                                   <  2:504029/100‑1 (MQ=255)
                                                             cGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGttt                                   <  2:564960/99‑1 (MQ=255)
                                                              ggATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtc                                 <  1:1898288/100‑1 (MQ=255)
                                                              ccgatctgAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtc                                 <  1:1789277/93‑1 (MQ=255)
                                                               gATGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtct                                >  1:1070488/1‑100 (MQ=255)
                                                                 tGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCg                              <  2:263001/100‑1 (MQ=255)
                                                                 tGCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCg                              <  1:1470846/100‑1 (MQ=255)
                                                                  gCAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGa                             >  1:1625357/1‑100 (MQ=255)
                                                                   cAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                            <  2:533635/100‑1 (MQ=255)
                                                                   cAGAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGATTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGAt                            >  1:1094989/1‑100 (MQ=255)
                                                                     gAAACTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTtcagatcgg                          >  2:1789288/1‑93 (MQ=255)
                                                                       aaCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCt                        >  1:1593616/1‑100 (MQ=255)
                                                                        aCTCAGGCGATTGCCGAAGCATTGATGGTGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCtt                       <  1:974876/100‑1 (MQ=255)
                                                                        aCTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCtt                       <  1:1745734/100‑1 (MQ=255)
                                                                         cTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTc                      <  2:1593627/100‑1 (MQ=255)
                                                                         cTCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTc                      >  1:126686/1‑100 (MQ=255)
                                                                          tCAGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCa                     <  1:704352/100‑1 (MQ=255)
                                                                            aGGCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATc                   <  2:1783141/100‑1 (MQ=255)
                                                                             ggCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATcc                  >  2:1863510/1‑100 (MQ=255)
                                                                              gCGATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCa                 <  1:1544691/100‑1 (MQ=255)
                                                                                gATTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAgg               >  2:732320/1‑100 (MQ=255)
                                                                                 aTTGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAggg              <  1:1471766/100‑1 (MQ=255)
                                                                                   tGCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtg            >  2:792551/1‑100 (MQ=255)
                                                                                    gCCGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtga           >  2:950558/1‑100 (MQ=255)
                                                                                     ccGAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgat          >  1:1733516/1‑100 (MQ=255)
                                                                                       gAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatgn        <  1:451600/100‑2 (MQ=255)
                                                                                       gAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatga        >  1:1061963/1‑100 (MQ=255)
                                                                                       gAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatga        >  1:192038/1‑100 (MQ=255)
                                                                                       gAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatga        >  1:84419/1‑100 (MQ=255)
                                                                                       gAAGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGtgatga        >  2:1866061/1‑100 (MQ=255)
                                                                                         aGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGAcc      <  2:1149957/100‑1 (MQ=255)
                                                                                         aGCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGAcc      >  2:1272274/1‑100 (MQ=255)
                                                                                          gCATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCg     <  2:979915/100‑1 (MQ=255)
                                                                                           cATTGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCgg    >  1:501573/1‑100 (MQ=255)
                                                                                             ttGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCt  <  2:556410/100‑1 (MQ=255)
                                                                                             ttGATGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCt  >  2:1522114/1‑100 (MQ=255)
                                                                                                |                                                                                                
CAGCGGCGCGTGATTATGGCATTCCAGATGCATTTCCGCCCGACGTTATATAACGGCGAAGCGGATGCAGAAACTCAGGCGATTGCCGAAGCATTGCTGGAGAAATACGGGCAGGATTAGCGCGGCAGTTTTCCGTGGTCGCGTAGCCAGGCGGCAGTTTTCTCGATACCTTCATCCAGGGTGATGACCGGCT  >  NC_000913/905624‑905816

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