Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,090,229 G→A 5.5% G159S (GGT→AGT)  yggI → Zn‑dependent metalloprotease‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,090,2290GA 5.5% 490.7 / 12.3 144G159S (GGT→AGT) yggIZn‑dependent metalloprotease‑related protein
Reads supporting (aligned to +/- strand):  ref base G (65/71);  new base A (4/4);  total (69/75)
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.49e-01

Only 100 of 152 total aligned reads displayed.

CCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGT  >  NC_000913/3090083‑3090375
                                                                                                                                                  |                                                                                                                                                  
ccgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                   <  4:896144/148‑1 (MQ=255)
ccgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                   >  1:853300/1‑148 (MQ=255)
              cATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTg                                                                                                                                     >  1:326443/1‑148 (MQ=255)
                tCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCg                                                                                                                                   >  2:84280/1‑148 (MQ=255)
                 cAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCTCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                   <  1:847645/131‑1 (MQ=255)
                 cAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGa                                                                                                                                  >  1:596357/1‑148 (MQ=255)
                    ttCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAAt                                                                                                                               >  3:904988/1‑148 (MQ=255)
                     tCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGa                                                                                                                                  >  3:373758/1‑144 (MQ=255)
                     tCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGa                                                                                                                                  <  4:373758/144‑1 (MQ=255)
                         aCTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGt                                                                                                                                       <  4:713149/135‑1 (MQ=255)
                         aCTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGt                                                                                                                                       >  3:713149/1‑135 (MQ=255)
                          cTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCAt                                                                                                                         >  3:503989/1‑148 (MQ=255)
                           tGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATc                                                                                                                        <  2:853300/148‑1 (MQ=255)
                           tGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATc                                                                                                                        <  2:979267/148‑1 (MQ=255)
                           tGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATc                                                                                                                        <  2:326443/148‑1 (MQ=255)
                            gCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCt                                                                                                                                          <  1:264306/129‑1 (MQ=255)
                            gCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCt                                                                                                                                          >  2:264306/1‑129 (MQ=255)
                            gCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCt                                                                                                                       >  4:576159/1‑148 (MQ=255)
                              aaTCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGa                                                                                                                     <  3:576159/148‑1 (MQ=255)
                               aTCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATc                                                                                                                        >  2:710501/1‑144 (MQ=255)
                                tCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACa                                                                                                                                             >  3:953616/1‑122 (MQ=255)
                                tCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACa                                                                                                                                             <  4:953616/122‑1 (MQ=255)
                                 ccGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                   <  1:280477/115‑1 (MQ=255)
                                 ccGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                   >  2:280477/1‑115 (MQ=255)
                                  cGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTa                                                                                                                 <  4:503989/148‑1 (MQ=255)
                                         cGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGt                                                                                                                                       <  2:872688/119‑1 (MQ=255)
                                         cGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGt                                                                                                                                       >  1:872688/1‑119 (MQ=255)
                                            aaCACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAAc                                                                                                                                              <  3:58963/109‑1 (MQ=255)
                                            aaCACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAAc                                                                                                                                              >  4:58963/1‑109 (MQ=255)
                                             acacCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAc                                                                                                                            >  1:5501/1‑126 (MQ=255)
                                             acacCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAc                                                                                                                            <  2:5501/126‑1 (MQ=255)
                                             acacCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACttt                                                                                                      <  1:620250/148‑1 (MQ=255)
                                                  ttCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGaaa                                                                                                                                >  2:123155/1‑117 (MQ=255)
                                                  ttCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGaaa                                                                                                                                <  1:123155/117‑1 (MQ=255)
                                                  ttCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAAt                                                                                                                               <  4:187746/118‑1 (MQ=255)
                                                   tCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAAAGTTCCTGAt                                                                                                <  3:308105/148‑1 (MQ=255)
                                                            gCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGAt                                                                                       <  2:596357/148‑1 (MQ=255)
                                                             cTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTAt                                                                                                                <  4:846604/122‑1 (MQ=255)
                                                             cTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTAt                                                                                                                >  3:846604/1‑122 (MQ=255)
                                                               gCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGc                                                                                    <  4:904988/148‑1 (MQ=255)
                                                                    tGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGc                                                                                    >  1:1053657/1‑143 (MQ=255)
                                                                    tGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGc                                                                                    <  2:1053657/143‑1 (MQ=255)
                                                                        aGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATaa                                                                                                                             <  3:970226/98‑1 (MQ=255)
                                                                          gAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGAt                                                                                       <  1:701711/134‑1 (MQ=255)
                                                                              aTCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCa                                                                                                                          >  2:53340/1‑95 (MQ=255)
                                                                              aTCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCa                                                                                                                          <  1:53340/95‑1 (MQ=255)
                                                                              aTCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAAc                                                                                                                   >  1:127503/1‑102 (MQ=255)
                                                                              aTCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAAc                                                                                                                   <  2:127503/102‑1 (MQ=255)
                                                                               tCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGa                                                                                                           >  2:644716/1‑109 (MQ=255)
                                                                                cAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATc                                                                                                                        >  3:836644/1‑95 (MQ=255)
                                                                                cAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATc                                                                                                                        <  4:836644/95‑1 (MQ=255)
                                                                                   cTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATaa                                                                                                                             >  4:457186/1‑87 (MQ=255)
                                                                                   cTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATaa                                                                                                                             <  3:457186/87‑1 (MQ=255)
                                                                                   cTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                 <  2:1019416/115‑1 (MQ=255)
                                                                                   cTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                 >  1:1019416/1‑115 (MQ=255)
                                                                                    ttACCGTGCGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGa                                                                                                           >  3:1051190/1‑104 (MQ=255)
                                                                                    ttACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGa                                                                                                           <  4:1051190/104‑1 (MQ=255)
                                                                                      aCCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACg                                                             >  1:65257/1‑148 (MQ=255)
                                                                                                      aTCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                                                    <  2:872769/93‑1 (MQ=255)
                                                                                                        cgcgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTgg                                                                                            <  2:451145/99‑1 (MQ=255)
                                                                                                        cgcgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTgg                                                                                            >  1:451145/1‑99 (MQ=255)
                                                                                                        cgcgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAac                                           >  1:69768/1‑148 (MQ=255)
                                                                                                          cgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAc                                                                                                                            <  4:151161/65‑1 (MQ=255)
                                                                                                          cgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAc                                                                                                                            >  3:151161/1‑65 (MQ=255)
                                                                                                          cgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTa                                             >  2:168271/1‑144 (MQ=255)
                                                                                                          cgTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTa                                             <  1:168271/144‑1 (MQ=255)
                                                                                                           gTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTg                                                                                                  >  3:370870/1‑90 (MQ=255)
                                                                                                           gTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTg                                                                                                  <  4:370870/90‑1 (MQ=255)
                                                                                                                        aGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGAt                                                                                       <  1:895426/88‑1 (MQ=255)
                                                                                                                            cGTCTATCGCTGTATTCACTGCGGTGAACAGCTGGTTGCGAAATAACCa                                                                                                                          <  3:454916/49‑1 (MQ=255)
                                                                                                                            cGTCTATCGCTGTATTCACTGCGGTGAACAGCTGGTTGCGAAATAACCa                                                                                                                          >  4:454916/1‑49 (MQ=255)
                                                                                                                                 aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGttt                                               <  4:112579/119‑1 (MQ=255)
                                                                                                                                 aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGttt                                               >  3:112579/1‑119 (MQ=255)
                                                                                                                                  tCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaa                                                                            >  1:671722/1‑89 (MQ=255)
                                                                                                                                  tCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaa                                                                            <  2:671722/89‑1 (MQ=255)
                                                                                                                                  tCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTggat                                                    >  2:842741/1‑111 (MQ=255)
                                                                                                                                  tCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTggat                                                    <  1:842741/113‑3 (MQ=255)
                                                                                                                                   cGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg                                                          >  3:532451/1‑106 (MQ=255)
                                                                                                                                   cGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg                                                          <  4:532451/106‑1 (MQ=255)
                                                                                                                                     cTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGc                                                                                           >  2:952136/1‑71 (MQ=255)
                                                                                                                                     cTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGc                                                                                           <  1:952136/71‑1 (MQ=255)
                                                                                                                                     cTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTg                                                                                         <  2:159611/73‑1 (MQ=255)
                                                                                                                                     cTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTg                                                                                         >  1:159611/1‑73 (MQ=255)
                                                                                                                                       gtgtTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt                                                 <  2:262898/111‑1 (MQ=255)
                                                                                                                                       gtgtTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt                                                 >  1:262898/1‑111 (MQ=255)
                                                                                                                                         gtTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGt                                >  4:623384/1‑126 (MQ=255)
                                                                                                                                         gtTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGt                                <  3:623384/126‑1 (MQ=255)
                                                                                                                                         gtTCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt                                                 <  2:463791/109‑1 (MQ=255)
                                                                                                                                         gtTCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt                                                 >  1:463791/1‑109 (MQ=255)
                                                                                                                                         gtTCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                       <  3:1081879/119‑1 (MQ=255)
                                                                                                                                         gtTCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                       >  4:1081879/1‑119 (MQ=255)
                                                                                                                                           tCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATGTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGtt                                                <  1:730826/108‑1 (MQ=255)
                                                                                                                                           tCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGtt                                                >  2:730826/1‑108 (MQ=255)
                                                                                                                                           tCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                       >  4:708819/1‑117 (MQ=255)
                                                                                                                                           tCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACgg                                       <  3:708819/117‑1 (MQ=255)
                                                                                                                                              cTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt                                                 <  4:393444/104‑1 (MQ=255)
                                                                                                                                              cTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt                                                 >  3:393444/1‑104 (MQ=255)
                                                                                                                                                 cGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCggt  <  2:69768/148‑1 (MQ=255)
                                                                                                                                                  ggTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATa                                                                                 <  3:646928/68‑1 (MQ=255)
                                                                                                                                                  ggTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATa                                                                                 >  4:646928/1‑68 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGT  >  NC_000913/3090083‑3090375

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: