Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,485,819 C→T A45T (GCC→ACC)  pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,485,8190CT99.2% 462.5 / ‑4.9 123A45T (GCC→ACC) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (68/54);  total (68/55)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 124 total aligned reads displayed.

TGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGA  >  NC_000913/4485671‑4485966
                                                                                                                                                    |                                                                                                                                                   
   aTAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTg                                                                                                                                                 >  3:915568/1‑150 (MQ=255)
     aaCCTGCTTGTACTGACGCTCATCCAGCTCGCGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGAt                                                                                                                                               <  2:1207841/150‑1 (MQ=255)
     aaCCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATg                                                                                                                                              <  2:787402/151‑1 (MQ=255)
          gCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTAccc                                                                                                                                         >  2:8980/1‑151 (MQ=255)
           cTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTAccc                                                                                                                                         >  2:933289/1‑150 (MQ=255)
            ttGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAACGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCAt                                                                                                                                       >  4:312882/1‑151 (MQ=255)
             tGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCAt                                                                                                                                       <  2:1060439/150‑1 (MQ=255)
             tGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCAt                                                                                                                                       <  3:250282/150‑1 (MQ=255)
             tGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCAt                                                                                                                                       <  3:1365784/150‑1 (MQ=255)
             tGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATc                                                                                                                                      >  3:687547/1‑151 (MQ=255)
              gTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCg                                                                                                                                     >  4:113153/1‑151 (MQ=255)
               tACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGc                                                                                                                                    >  3:1154996/1‑151 (MQ=255)
                aCTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCt                                                                                                                                   <  2:443559/151‑1 (MQ=255)
                   gACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGAt                                                                                                                                <  1:129560/151‑1 (MQ=255)
                    aCGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGAtt                                                                                                                               <  3:663396/151‑1 (MQ=255)
                     cGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGAttt                                                                                                                              >  2:818999/1‑151 (MQ=255)
                        tCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTAt                                                                                                                           >  3:294660/1‑151 (MQ=255)
                         cATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATc                                                                                                                          >  2:1180605/1‑151 (MQ=255)
                            ccAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGa                                                                                                                        <  3:394171/150‑1 (MQ=255)
                             cAGCTCACGTTCTTTGCGGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAg                                                                                                                       >  3:643240/1‑150 (MQ=255)
                               gCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCt                                                                                                                     >  4:48868/1‑150 (MQ=255)
                               gCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTg                                                                                                                    <  2:1282589/151‑1 (MQ=255)
                                 tCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGtt                                                                                                                  >  1:124313/1‑151 (MQ=255)
                                    cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTg                                                                                                               <  2:890058/151‑1 (MQ=255)
                                    cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTg                                                                                                               <  1:273273/151‑1 (MQ=255)
                                    cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTg                                                                                                               >  2:928370/1‑151 (MQ=255)
                                    cGTTCTCTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGTAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGt                                                                                                                   <  4:419639/147‑1 (MQ=255)
                                         tttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAAt                                                                                                           >  3:840595/1‑150 (MQ=255)
                                          ttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATc                                                                                                          <  2:848172/150‑1 (MQ=255)
                                          ttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATc                                                                                                          <  2:24064/150‑1 (MQ=255)
                                            gCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                        >  2:173124/1‑150 (MQ=255)
                                            gCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                        >  3:663834/1‑150 (MQ=255)
                                             ccGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGgag                                                                                                      <  1:267775/151‑1 (MQ=255)
                                             ccGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGga                                                                                                       <  1:20165/150‑1 (MQ=255)
                                              cGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCGGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                        >  2:589733/1‑148 (MQ=255)
                                              cGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgggg                                                                                                      >  2:270904/1‑148 (MQ=255)
                                              cGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                        >  3:628343/1‑148 (MQ=255)
                                               gCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                                    <  2:259211/151‑1 (MQ=255)
                                                cAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaaa                                                                                                   <  1:137851/151‑1 (MQ=255)
                                                  aCCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAgg                                                                                                 >  2:812234/1‑151 (MQ=255)
                                                    cAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCg                                                                                               >  4:125156/1‑151 (MQ=255)
                                                    cAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCg                                                                                               <  1:147001/151‑1 (MQ=255)
                                                     aaTAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCg                                                                                               >  2:663686/1‑150 (MQ=255)
                                                          ggAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTgg                                                                                         <  2:1164482/151‑1 (MQ=255)
                                                             gaATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTc                                                                                      >  3:821387/1‑151 (MQ=255)
                                                             gaATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTc                                                                                      >  3:1172988/1‑151 (MQ=255)
                                                              aaTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCg                                                                                     >  1:218787/1‑151 (MQ=255)
                                                                 tCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCg                                                                                     >  4:451152/1‑148 (MQ=255)
                                                                 tCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGa                                                                                    >  2:1057963/1‑149 (MQ=255)
                                                                 tCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGa                                                                                    >  3:269013/1‑149 (MQ=255)
                                                                 tCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGa                                                                                    >  3:1102114/1‑149 (MQ=255)
                                                                  cGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAg                                                                                  <  3:865934/150‑1 (MQ=255)
                                                                   gCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGAc                                                                                >  2:945370/1‑151 (MQ=255)
                                                                         aaaGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCac                                                                          >  2:822812/1‑151 (MQ=255)
                                                                           aGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacga                                                                        <  3:1131251/151‑1 (MQ=255)
                                                                              aCATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacgacg                                                                      >  2:428949/1‑150 (MQ=255)
                                                                               cATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGAt                                                                    >  2:865213/1‑151 (MQ=255)
                                                                                   cGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCa                                                                 <  3:1224590/150‑1 (MQ=255)
                                                                                      aaCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                                              >  2:517545/1‑150 (MQ=255)
                                                                                      aaCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                                              >  1:384019/1‑150 (MQ=255)
                                                                                      aaCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                                              >  1:48704/1‑150 (MQ=255)
                                                                                       aCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCAc                                                            >  1:114756/1‑151 (MQ=255)
                                                                                         aTGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTc                                                          >  2:632726/1‑151 (MQ=255)
                                                                                          tGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTcc                                                         >  2:500037/1‑151 (MQ=255)
                                                                                           ggTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCg                                                        <  3:865465/151‑1 (MQ=255)
                                                                                           ggTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCg                                                        <  3:8981/151‑1 (MQ=255)
                                                                                                 gcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCGGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGtt                                                   <  3:589762/150‑1 (MQ=255)
                                                                                                   aacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTTTTGAAATGTTCTGCAGTCCAAGCGAGGCCACATGGTTGGAGGCCGCCCTCTCCCCTGCAGGCCCTCCTCCGttgc                                                 <  3:365928/150‑3 (MQ=255)
                                                                                                     caaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTcc                                              <  2:317194/151‑1 (MQ=255)
                                                                                                       aTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCgg                                            <  2:36137/151‑1 (MQ=255)
                                                                                                         gTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCggg                                           <  1:65638/150‑1 (MQ=255)
                                                                                                         gTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGc                                          <  3:392709/151‑1 (MQ=255)
                                                                                                            tGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTAc                                       >  2:1302151/1‑151 (MQ=255)
                                                                                                             gCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCt                                         <  4:338939/148‑1 (MQ=255)
                                                                                                             gCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTAcc                                      >  2:103740/1‑151 (MQ=255)
                                                                                                              ccccGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCg                                     <  3:819034/151‑1 (MQ=255)
                                                                                                                 cGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCtt                                  <  4:291623/151‑1 (MQ=255)
                                                                                                                 cGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCtt                                  >  2:119326/1‑151 (MQ=255)
                                                                                                                       tccttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac                             >  2:1372605/1‑150 (MQ=255)
                                                                                                                        ccttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac                           <  2:962493/151‑1 (MQ=255)
                                                                                                                         cttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                          >  3:130699/1‑151 (MQ=255)
                                                                                                                          ttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                          >  2:647385/1‑150 (MQ=255)
                                                                                                                          ttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                          >  2:1289613/1‑150 (MQ=255)
                                                                                                                            ccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa                        <  3:173779/150‑1 (MQ=255)
                                                                                                                            ccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa                         <  2:663812/149‑1 (MQ=255)
                                                                                                                              aGTTCGCCCCGACGTAGCAGGGTGCTGATGTGCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt                     <  1:451156/151‑1 (MQ=255)
                                                                                                                              aGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt                     >  2:666300/1‑151 (MQ=255)
                                                                                                                                 tCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc                   <  4:218785/150‑1 (MQ=255)
                                                                                                                                 tCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc                   >  1:60744/1‑150 (MQ=255)
                                                                                                                                  caccccGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt                 >  2:338943/3‑151 (MQ=255)
                                                                                                                                  cGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt                 <  2:847111/151‑1 (MQ=255)
                                                                                                                                  cGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCAt                 >  4:452401/1‑151 (MQ=255)
                                                                                                                                   gCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg                >  2:610657/1‑151 (MQ=255)
                                                                                                                                    ccccGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGc               >  2:568580/1‑151 (MQ=255)
                                                                                                                                     cccGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGc               >  3:708112/1‑150 (MQ=255)
                                                                                                                                     cccGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCa              <  3:1334848/151‑1 (MQ=255)
                                                                                                                                         aCGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTAc          <  2:294633/151‑1 (MQ=255)
                                                                                                                                             aGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCt       <  3:812268/150‑1 (MQ=255)
                                                                                                                                               cAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCt    <  3:928404/151‑1 (MQ=255)
                                                                                                                                               cAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCt    >  2:335335/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGA  >  NC_000913/4485671‑4485966

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