Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,642,135 C→A A425D (GCT→GAT)  dtpB → dipeptide/tripeptide:H(+) symporter DtpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,642,1350CA97.1% 243.3 / 0.1 70A425D (GCT→GAT) dtpBdipeptide/tripeptide:H(+) symporter DtpB
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (36/32);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GTTTATCGTGCTGGTGTACTTATTCCAGAGCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTG  >  NC_000913/3642001‑3642266
                                                                                                                                      |                                                                                                                                   
gTTTATCGTGCTGGTGTACTTATTCCAGAGCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCTGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCcg                                                                                                                                <  1:456757/140‑1 (MQ=255)
          cTGGTGTACTTATTCCAGAGCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTTGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCt                                                                                                                      >  1:463653/1‑140 (MQ=255)
                tACTTATTCCAGAGCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGggcgg                                                                                                                >  1:512122/1‑140 (MQ=255)
                       tCCAGAGCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTAtgtg                                                                                                         >  2:558865/1‑140 (MQ=255)
                        ccAGAGCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTgg                                                                                                        >  1:115481/1‑140 (MQ=255)
                         cAGAGCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCAGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGc                                                                                                       <  2:223075/140‑1 (MQ=255)
                            agCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGc                                                                                                       <  2:862652/137‑1 (MQ=255)
                            agCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGc                                                                                                       >  1:862652/1‑137 (MQ=255)
                            agCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTAACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAAc                                                                                                    >  2:438086/1‑140 (MQ=255)
                             gCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACa                                                                                                   >  1:546082/1‑140 (MQ=255)
                               ttAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACAtt                                                                                                 <  2:302587/140‑1 (MQ=255)
                                 aGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTa                                                                                               >  2:873147/1‑140 (MQ=255)
                                 aGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTa                                                                                               >  2:567748/1‑140 (MQ=255)
                                 aGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTa                                                                                               >  1:302795/1‑140 (MQ=255)
                                    tGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACcg                                                                                            <  2:138122/140‑1 (MQ=255)
                                      aaCTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACcgcg                                                                                          >  2:96952/1‑140 (MQ=255)
                                           tttATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGcc                                                                                     >  2:570599/1‑140 (MQ=255)
                                           tttATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGcc                                                                                     >  1:179694/1‑140 (MQ=255)
                                            ttattaGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCg                                                                                    <  1:558865/140‑1 (MQ=255)
                                             tattaGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCgg                                                                                   >  2:207442/1‑140 (MQ=255)
                                             tattaGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCgg                                                                                   >  2:278619/1‑140 (MQ=255)
                                             tattaGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCgg                                                                                   <  1:567748/140‑1 (MQ=255)
                                                 aGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACcg                                                                                            >  2:174735/1‑127 (MQ=255)
                                                 aGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACcg                                                                                            <  1:174735/127‑1 (MQ=255)
                                                   cgcCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGacaaca                                                                             >  1:372114/1‑140 (MQ=255)
                                                     cccTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGa                                                                                  <  2:127780/133‑1 (MQ=255)
                                                     cccTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGa                                                                                  >  1:127780/1‑133 (MQ=255)
                                                     cccTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACAtt                                                                           >  1:276088/1‑140 (MQ=255)
                                                     cccTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACAtt                                                                           >  2:348865/1‑140 (MQ=255)
                                                     cccTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACAtt                                                                           >  1:343212/1‑140 (MQ=255)
                                                          ggCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGa                                                                      >  2:854256/1‑140 (MQ=255)
                                                          ggCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGa                                                                      <  2:182420/140‑1 (MQ=255)
                                                                gCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCt                                                                >  1:819595/1‑140 (MQ=255)
                                                                     gatTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTgaga                                                           >  1:595040/1‑140 (MQ=255)
                                                                          ctgcCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTg                                                      >  2:317435/1‑140 (MQ=255)
                                                                          ctgcCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTg                                                      <  2:302795/140‑1 (MQ=255)
                                                                           tgcCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGc                                                     <  2:595040/140‑1 (MQ=255)
                                                                           tgcCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGc                                                     <  2:785741/140‑1 (MQ=255)
                                                                           tgcCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGc                                                     <  2:372114/140‑1 (MQ=255)
                                                                               cTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGt                                                 <  2:196562/140‑1 (MQ=255)
                                                                                  gTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCta                                              >  1:662874/1‑140 (MQ=255)
                                                                                  gTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGAAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCta                                              >  1:192408/1‑140 (MQ=255)
                                                                                      cGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATAcc                                          <  1:514660/140‑1 (MQ=255)
                                                                                       gcagcaTTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGt                                                 >  1:312690/1‑132 (MQ=255)
                                                                                       gcagcaTTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGt                                                 <  2:312690/132‑1 (MQ=255)
                                                                                       gcagcaTTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCa                                         >  1:138446/1‑140 (MQ=255)
                                                                                        cagcaTTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCta                                              <  2:377234/134‑1 (MQ=255)
                                                                                        cagcaTTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCta                                              >  1:377234/1‑134 (MQ=255)
                                                                                        cagcaTTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCaa                                        <  1:442940/140‑1 (MQ=255)
                                                                                        cagcaTTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCaa                                        <  2:372948/140‑1 (MQ=255)
                                                                                            aTTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACgtg                                    <  2:115481/140‑1 (MQ=255)
                                                                                             tttGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACgtgt                                   <  2:512122/140‑1 (MQ=255)
                                                                                                     gCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGttt                                 >  2:440536/1‑134 (MQ=255)
                                                                                                     gCTTTATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGttt                                 <  1:440536/134‑1 (MQ=255)
                                                                                                       tttATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGAt                         >  1:358287/1‑140 (MQ=255)
                                                                                                       tttATTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGAt                         <  2:276088/140‑1 (MQ=255)
                                                                                                          aTTCTCGGGATGTGGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATtgg                      <  1:873147/140‑1 (MQ=255)
                                                                                                                     gtgGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGc           >  1:815358/1‑140 (MQ=255)
                                                                                                                      tgGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATtggt                     >  2:366002/1‑129 (MQ=255)
                                                                                                                      tgGTTCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATtggt                     <  1:366002/129‑1 (MQ=255)
                                                                                                                         ttCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTggg       <  2:819595/140‑1 (MQ=255)
                                                                                                                          tCCTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGc      >  2:228013/1‑140 (MQ=255)
                                                                                                                           ccTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCg     <  1:570599/140‑1 (MQ=255)
                                                                                                                            cTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATtggt                     >  1:284429/1‑123 (MQ=255)
                                                                                                                            cTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATtggt                     <  2:284429/123‑1 (MQ=255)
                                                                                                                            cTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGt    <  1:228013/140‑1 (MQ=255)
                                                                                                                            cTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGt    >  1:861054/1‑140 (MQ=255)
                                                                                                                            cTGACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGt    <  1:130187/140‑1 (MQ=255)
                                                                                                                              gACGCAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTg  <  2:358287/140‑1 (MQ=255)
                                                                                                                              gACGAAGGATGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTg  <  2:192408/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GTTTATCGTGCTGGTGTACTTATTCCAGAGCTTAGGTGAACTGTTTATTAGCGCCCTTGGCCTGGCGATGATTGCTGCCCTGGTGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTG  >  NC_000913/3642001‑3642266

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