Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,389,396 T→C L156P (CTT→CCT)  ycjG → L‑Ala‑D/L‑Glu epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,389,3960TC99.2% 505.5 / ‑4.5 122L156P (CTT→CCT) ycjGL‑Ala‑D/L‑Glu epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (69/52);  total (69/53)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 122 total aligned reads displayed.

TGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCTTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGAACAACCGCT  >  NC_000913/1389249‑1389543
                                                                                                                                                   |                                                                                                                                                   
tggctgATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTAt                                                                                                                                                  >  2:1911865/1‑151 (MQ=255)
 ggctgATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTa                                                                                                                                                   <  2:1696897/149‑1 (MQ=255)
 ggctgATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTAt                                                                                                                                                  <  1:1745229/150‑1 (MQ=255)
 ggctgATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATc                                                                                                                                                 >  1:906929/1‑151 (MQ=255)
  gctgATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCa                                                                                                                                                <  1:656461/151‑1 (MQ=255)
  gctgATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCa                                                                                                                                                <  1:759576/151‑1 (MQ=255)
    tgATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGt                                                                                                                                              >  1:1736053/1‑151 (MQ=255)
     gATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAg                                                                                                                                               <  1:541479/149‑1 (MQ=255)
        ttGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAg                                                                                                                                           <  1:8771/150‑1 (MQ=255)
          gATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCgg                                                                                                                                        <  2:803813/151‑1 (MQ=255)
          gATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCgg                                                                                                                                        >  2:1622656/1‑151 (MQ=255)
          gATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCgg                                                                                                                                        >  2:632903/1‑151 (MQ=255)
                cATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGAtggtg                                                                                                                                  >  2:1383248/1‑151 (MQ=255)
                 aTAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGAtggtgg                                                                                                                                 >  1:1079968/1‑151 (MQ=255)
                   aaCGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGc                                                                                                                                >  1:1489412/1‑150 (MQ=255)
                    aCGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCa                                                                                                                               >  2:768116/1‑150 (MQ=255)
                      gCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAAtt                                                                                                                            >  2:284122/1‑151 (MQ=255)
                       cTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTc                                                                                                                           >  1:630002/1‑151 (MQ=255)
                         tCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGc                                                                                                                         >  2:1732330/1‑151 (MQ=255)
                           ccGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGcac                                                                                                                       >  1:289493/1‑151 (MQ=255)
                            cGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGcaca                                                                                                                      <  2:906737/151‑1 (MQ=255)
                              agaCAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGc                                                                                                                    <  2:23360/151‑1 (MQ=255)
                                  aGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCtgtg                                                                                                                >  2:420804/1‑151 (MQ=255)
                                  aGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCtgtg                                                                                                                >  1:748411/1‑151 (MQ=255)
                                    ttATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGcc                                                                                                              >  1:632885/1‑151 (MQ=255)
                                     tATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGccc                                                                                                             >  2:598203/1‑151 (MQ=255)
                                     tATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGccc                                                                                                             >  1:682298/1‑151 (MQ=255)
                                     tATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGccc                                                                                                             >  1:354225/1‑151 (MQ=255)
                                      aTCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGccc                                                                                                             <  2:394748/150‑1 (MQ=255)
                                      aTCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCg                                                                                                            >  1:1033046/1‑151 (MQ=255)
                                        cACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGa                                                                                                           <  2:740698/150‑1 (MQ=255)
                                        cACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGAt                                                                                                          <  1:1162052/151‑1 (MQ=255)
                                          cttcacaGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGc                                                                                                        >  2:756797/4‑151 (MQ=255)
                                          cTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGc                                                                                                        >  1:1480665/1‑151 (MQ=255)
                                            gCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGc                                                                                                        <  2:729868/149‑1 (MQ=255)
                                            gCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCg                                                                                                       >  2:132670/1‑150 (MQ=255)
                                            gCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGa                                                                                                      <  2:184386/151‑1 (MQ=255)
                                             cacaGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGAc                                                                                                     <  2:78783/151‑1 (MQ=255)
                                             cacaGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGAc                                                                                                     >  2:205077/1‑151 (MQ=255)
                                                aGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCt                                                                                                  <  2:1059563/151‑1 (MQ=255)
                                                  aCGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGa                                                                                                <  2:494713/151‑1 (MQ=255)
                                                  aCGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGa                                                                                                <  2:623915/151‑1 (MQ=255)
                                                      ttGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGt                                                                                            >  1:1960968/1‑151 (MQ=255)
                                                       tGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGtt                                                                                           >  1:1488228/1‑151 (MQ=255)
                                                         tCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTg                                                                                          >  1:952590/1‑150 (MQ=255)
                                                             cGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGc                                                                                      >  2:62335/1‑150 (MQ=255)
                                                                        aTCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTg                                                                          >  2:1933390/1‑151 (MQ=255)
                                                                          cAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGc                                                                        >  2:141014/1‑151 (MQ=255)
                                                                               ggCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCa                                                                   <  1:1958879/151‑1 (MQ=255)
                                                                                 ccAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGc                                                                    >  2:509182/1‑148 (MQ=255)
                                                                                    aTAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAgg                                                              <  2:150347/151‑1 (MQ=255)
                                                                                    aTAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAgg                                                              <  2:1340499/151‑1 (MQ=255)
                                                                                      aGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAggg                                                             <  1:179818/150‑1 (MQ=255)
                                                                                      aGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGt                                                            >  1:1780391/1‑151 (MQ=255)
                                                                                      aGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGt                                                            >  1:111179/1‑151 (MQ=255)
                                                                                       gTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGtt                                                           <  1:543182/151‑1 (MQ=255)
                                                                                          cATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTggc                                                        >  1:28477/1‑151 (MQ=255)
                                                                                             cAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTggcggc                                                     <  1:1714330/151‑1 (MQ=255)
                                                                                             cAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTggcggc                                                     >  1:305422/1‑151 (MQ=255)
                                                                                              aaCACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTggcggc                                                     >  2:1511787/1‑150 (MQ=255)
                                                                                              aaCACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGgcggcg                                                    <  1:540973/151‑1 (MQ=255)
                                                                                               acacTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGgcgc                                                   >  2:817090/1‑151 (MQ=255)
                                                                                                 acTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGt                                                 >  2:1724819/1‑151 (MQ=255)
                                                                                                 acTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGt                                                 >  1:1616457/1‑151 (MQ=255)
                                                                                                    ctGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGtt                                                <  2:1661165/149‑1 (MQ=255)
                                                                                                     tGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGcc                                             >  2:1585277/1‑151 (MQ=255)
                                                                                                      ggcaggcaGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGcc                                             >  1:725591/1‑150 (MQ=255)
                                                                                                          ggcaggcGCGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTa                                        <  1:450405/151‑1 (MQ=255)
                                                                                                               gcgcGAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGc                                   >  1:467379/1‑151 (MQ=255)
                                                                                                                cgcgAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCg                                  >  2:1064817/1‑151 (MQ=255)
                                                                                                                cgcgAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCg                                  <  1:1938191/151‑1 (MQ=255)
                                                                                                                 gcgAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCgg                                 >  1:14268/1‑151 (MQ=255)
                                                                                                                 gcgAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCgg                                 >  1:1911048/1‑151 (MQ=255)
                                                                                                                 gcgAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCgg                                 >  2:461515/1‑151 (MQ=255)
                                                                                                                  cgAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGa                                <  2:1639848/151‑1 (MQ=255)
                                                                                                                   gAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGAt                               >  2:1661724/1‑151 (MQ=255)
                                                                                                                   gAAATTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGAt                               >  1:835064/1‑151 (MQ=255)
                                                                                                                     aaTTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGAttt                             <  2:1465096/151‑1 (MQ=255)
                                                                                                                      aTTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGAttt                             >  2:157439/1‑150 (MQ=255)
                                                                                                                      aTTACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTa                            <  1:649060/151‑1 (MQ=255)
                                                                                                                        tACTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAgg                          >  1:440802/1‑151 (MQ=255)
                                                                                                                          cTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGc                         <  2:1000434/150‑1 (MQ=255)
                                                                                                                          cTGAAAGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCg                        >  1:1026609/1‑151 (MQ=255)
                                                                                                                              aaGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTg                     >  2:1016192/1‑150 (MQ=255)
                                                                                                                               aGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCg                   <  2:293345/151‑1 (MQ=255)
                                                                                                                               aGTGAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCg                   <  1:1734142/151‑1 (MQ=255)
                                                                                                                                  gAAGCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATg                <  1:1423573/151‑1 (MQ=255)
                                                                                                                                     gCTGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCtt             >  1:336221/1‑151 (MQ=255)
                                                                                                                                       tGGATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTg            >  2:215259/1‑150 (MQ=255)
                                                                                                                                        ggATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGa           <  1:811460/150‑1 (MQ=255)
                                                                                                                                         gATAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGaac         >  2:1370113/1‑151 (MQ=255)
                                                                                                                                          aTAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGaaca        <  1:436699/151‑1 (MQ=255)
                                                                                                                                           tAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGaaca        >  1:874752/1‑150 (MQ=255)
                                                                                                                                           tAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGaaca        >  1:895785/1‑150 (MQ=255)
                                                                                                                                           tAACCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGaaca        >  2:85594/1‑150 (MQ=255)
                                                                                                                                             aCCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGAACAAcc     >  1:652651/1‑151 (MQ=255)
                                                                                                                                             aCCATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGAACAAcc     <  1:1586949/151‑1 (MQ=255)
                                                                                                                                              ccATCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGAACAACCg    <  2:491712/151‑1 (MQ=255)
                                                                                                                                                aTCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGAACAACCGCt  <  1:1544897/151‑1 (MQ=255)
                                                                                                                                                aTCCTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGAACAACCGCt  >  2:1708151/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAATTACTGAAAGTGAAGCTGGATAACCATCTTATCAGTGAGCGGATGGTGGCAATTCGCACAGCTGTGCCCGATGCGACGCTGATCGTTGATGCAAATGAATCCTGGCGTGCAGAAGGGTTGGCGGCGCGTTGCCAGCTATTGGCGGATTTAGGCGTTGCGATGCTTGAACAACCGCT  >  NC_000913/1389249‑1389543

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