Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,396,772 T→C F236F (TTT→TTC)  amiB → N‑acetylmuramoyl‑L‑alanine amidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,396,7720TC94.8% 197.2 / NA 58F236F (TTT→TTC) amiBN‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (20/35);  minor base A (3/0);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 5.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTG  >  NC_000913/4396525‑4397012
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
gCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCTGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCaaa                                                                                                                                                                                                                                               <  2:174509/251‑1 (MQ=255)
   aaCCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGCCGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCACCGGTGGGACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCCACGACGATCCGATGTTCAAAggg                                                                                                                                                                                                                                            >  1:19811/1‑250 (MQ=255)
   aaCCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGc                                                                                                                                                                                                                                            >  1:218910/1‑251 (MQ=255)
                    tttAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACgg                                                                                                                                                                                                                           >  1:320757/1‑251 (MQ=255)
                           cGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGGACGCGGGAGAAAAATGTCCCCATCCCCCTCGCACGTAAATTGCGTACTTTGCTCGATTGCGAGCCTGTGTTCAAAGGCGGTTTACCCCGTGGCggggtttcc                                                                                                                                                                                                                    >  2:61538/1‑246 (MQ=255)
                                    aCCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCg                                                                                                                                                                                                           <  1:132391/251‑1 (MQ=255)
                                     ccGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCgg                                                                                                                                                                                                          >  1:318121/1‑251 (MQ=255)
                                         aCTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGAt                                                                                                                                                                                                      <  1:335059/251‑1 (MQ=255)
                                            aCGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATgg                                                                                                                                                                                                    <  1:419675/250‑1 (MQ=255)
                                               ggTGTTATCAGCAGTAATACGGTAACGCGGCCGGCAGCGCGCGGGACGGCTAACACTGGCGATAAAAGTATCAGCGCTAGTGAGGCCGGACACGGGGGTCAGGACGCTGGCGCTATCGGCCCCGGGGGTACGCGGGAGAAAAAGGTCACCATCGCCAGCGCGCGTAAAGTGCGTACGTTGCGCAATGACGATACGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGgcg                                                                                                                                                                                                <  1:186936/251‑1 (MQ=255)
                                                            ctAATACGGTAACGCCACCGGCCAAGCGCGCGACGAATCACACTGGAGATAAAATTATCATCGCTATTGATGCCGGACCAGGCGGTCAGGACCCTGGCGCTATTGGCACAGGTGGGACGCGGGGGAAAAAAGTCACCATTGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGCAGTACAGAGGTTCAAAGGCGTTTTAACCCGTGACTGGGAATACTTTGTCTCGGTGATGGCGGGCAGTGAtgcggac                                                                                                                                                                                   >  2:46949/2‑246 (MQ=255)
                                                              aaTACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACg                                                                                                                                                                                 <  1:153702/251‑1 (MQ=255)
                                                                   ggTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCGAGGCGCGGTAGGCGCCGGTGGTACGGGGGAGAAAAATGGCAACAGCGGCAACGCGGGTAAATTGCGGACGTTGCTAAATGAAGATACGAAGTTAAAAGGCGGTTTAACCCGGGAAGGGGATTACTTTATCTCGGGGATGGGGGGCAGTGATGGGGCAGGTa                                                                                                                                                                               >  1:161338/1‑248 (MQ=255)
                                                                   ggTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCCTCGCGCGTAAATTGCGGACTTTGCTCAATGGCGATCCGATGTTCAAAGGCGGTTTAACCCGGGACGGGGGTTACTTTATCTCGGTGATGGGGGGCAGTGATGTGGCACGGAAGc                                                                                                                                                                            >  2:115357/1‑251 (MQ=255)
                                                                   ggTAACGCGGCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTAGCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAAGTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGAAAGc                                                                                                                                                                            <  2:361617/251‑1 (MQ=255)
                                                                               ggCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACTCCAAtt                                                                                                                                                                <  1:332077/251‑1 (MQ=255)
                                                                                 cAGCGCGCGCGGCGGCTAACACTGGCGATAAAATTATCATCGCTATTGAGGCCGGACACGGCGGTCAGGGCCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCAGCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGGAAGCAAAACGCCAATTTc                                                                                                                                                              <  2:401089/251‑1 (MQ=255)
                                                                                  aGCGCGCGCGACGTCTACCACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGATCCTGGAGCTATCGGCCCCTGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTATATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCTTTTTAACCCGTGACGGGGATTACTTTATCTCTGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTcc                                                                                                                                                             <  2:351056/251‑1 (MQ=255)
                                                                                                  aaCACTGGCGATAAAATTATCATCGCTATTGCTGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGGACTTTGCTCAATGACGATCCGATGTTCAAAGGAGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTACTCGTGTCGATTCACGc                                                                                                                                             >  1:87543/1‑251 (MQ=255)
                                                                                                   acacTGGCGATAAAATTATCACCGCTATTGATGCCTGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTGCTTTGCGCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATTGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGAGTCACGCt                                                                                                                                            <  2:37084/251‑1 (MQ=255)
                                                                                                     acTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCCGGACCCTGGCGCTATCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTTACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTgg                                                                                                                                          >  2:78134/1‑250 (MQ=255)
                                                                                                        ggCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGAt                                                                                                                                         <  1:322685/249‑1 (MQ=255)
                                                                                                         gCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGcc                                                                                                                                      <  1:162291/251‑1 (MQ=255)
                                                                                                           gATAAAATTATCATCGCTAGTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTCGCGGCCCCGGTGGTGCGCTGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTGTAACCCGTGACGGGGATTACTTTATCTCGGGGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAAGTTCCTCGTGTCGATTCACGCTGATGCCGc                                                                                                                                    <  2:145092/251‑1 (MQ=255)
                                                                                                            aTAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCCAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGcc                                                                                                                                   >  2:213615/1‑250 (MQ=255)
                                                                                                                 aTTATCATCGCCATTAATGACGGACACAGCGGTCAGGACACTCGAGATACCGGCCCCGGTGGTACACGGGGCAAAAATCGCACCATCGCCATCCCGCGGAAATTGCGTACTTTGCTCAATGACGATCCGATGTTACAAGGCGATTTAACCCGGGACGGGGATTACTCTATCTCGGTGGTGGGGCGCAATGATGTGGCAAGTACGCAAAACGCCAATTCACTCGGGTCGATTTACGCTGAAGCCGCACCGa                                                                                                                               >  2:296030/1‑250 (MQ=255)
                                                                                                                      catcGCTATTGATGCCGTACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCGCAATGACGATCCGATGTTCAAAGGCGTGGTAACCCGGGACGGGGATTACTTGATCTCGGTGATGGGGCGCAGTGATGTGTCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCAGCACCTAACCGCa                                                                                                                         <  2:261327/251‑1 (MQ=255)
                                                                                                                         cGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTg                                                                                                                      <  1:109694/251‑1 (MQ=255)
                                                                                                                               tGATGCCGGACCCGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACt                                                                                                                 <  2:11079/250‑1 (MQ=255)
                                                                                                                                       gacgcggcggTCAGGACCCTGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAATTGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACTATCCGATGTTCAAAGGCGTGTTAACCTGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCTACTGGCGCTTcc                                                                                                        <  2:322463/247‑1 (MQ=255)
                                                                                                                                       gacacgggagtcaatgctcTGGCGCTAGCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTAATTTGCTCAAGGACGATCCGATGTTCAAAGGCGTTTTAACCCGGGACGGGGATTACGTTATCTCGGTGATGGGGCGCAGTGATGTGGCAGGTAAGCAAAACGCCAATTTCCTCGTGTCGAGTCACGCTGATGCCGCTCCGAACCGCAATGCGACTGGCGCTTcc                                                                                                        <  2:276974/233‑1 (MQ=255)
                                                                                                                                            ggcggTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATg                                                                                                   <  2:436726/251‑1 (MQ=255)
                                                                                                                                                        ccTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCCCCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTaa                                                                                       <  1:376445/251‑1 (MQ=255)
                                                                                                                                                        ccTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCCCCGAACCGCAGTGCGCCTGGCGCTTCCGAATGGGTGCTCTCTa                                                                                        >  2:916/1‑250 (MQ=255)
                                                                                                                                                              gcTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTcg                                                                                 <  1:77345/251‑1 (MQ=255)
                                                                                                                                                                    ggCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTAACGCTGATGCCGCACCGAACAGCAGTGCGCCTGGCGCTTCCGTCTGGGTGCTCTCTAACCGTCGCGCCa                                                                            >  2:239986/1‑250 (MQ=255)
                                                                                                                                                                         cGGTGGTACGCGGGAGAAAAATGGCCACATCGCAGTCCCGCGTAAATCGCGTACTTTGCCCTCTGACGACCAGCTGTTCAAAGGCGCTTTAACCCGTGCCGGGGATTCCCTTATCTCGGCGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGACAATTTACTCGTGTCGATTGACGATGATGCCGCACCGAACCGCAGGGAGACTGGCGCTTCCGTATCGTGGCTCTTTTCCCGTCGCGCCAATAGTg                                                                      >  2:63490/1‑251 (MQ=255)
                                                                                                                                                                          ggtggtACGCGGGAGAAAAATGTCACCGTCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCGCTAACCGTCGCGCCAACAGTGa                                                                     <  1:303609/251‑1 (MQ=255)
                                                                                                                                                                                   cgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTTGCGCCAACAGTGAAATGGCCAg                                                            <  2:320757/251‑1 (MQ=255)
                                                                                                                                                                                               gTCACCATCGCCATCGCGCGTAAATTGCGTTCTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCATCTGGCTGGagc                                                 <  2:293513/250‑1 (MQ=255)
                                                                                                                                                                                               gTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACCGGGAAATGGCCAGCTGGCTGGagc                                                 >  2:134818/1‑250 (MQ=255)
                                                                                                                                                                                                 cACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGagcag                                               <  1:213615/250‑1 (MQ=255)
                                                                                                                                                                                                  aCCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGGTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAgcagca                                             <  2:318121/251‑1 (MQ=255)
                                                                                                                                                                                                    catcgccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACg                                           >  1:307529/1‑251 (MQ=255)
                                                                                                                                                                                                     atcgccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGAGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGACGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACgg                                          >  1:286389/1‑250 (MQ=255)
                                                                                                                                                                                                         ccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGaac                                      >  1:338209/1‑250 (MQ=255)
                                                                                                                                                                                                           atcgcGCGTAAATTGCGTACTTTGCTCGATGACGAGCCGATGTTCAAAGGCGTGTTAACCCGTGTCGGGGATTACTTTATCTAGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTGCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCTCGTCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACa                                    <  1:194294/251‑1 (MQ=255)
                                                                                                                                                                                                              gcgcgTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCGTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTc                                 <  1:274176/251‑1 (MQ=255)
                                                                                                                                                                                                                     aaTTGCGTGCTTTGCTCAATGACGATCCGACGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAgctg                          <  2:243359/251‑1 (MQ=255)
                                                                                                                                                                                                                        tGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAgctgct                        <  1:74494/250‑1 (MQ=255)
                                                                                                                                                                                                                           ggacgtggCTCGATGACGAGCCGATGGTCAAAGGCTTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCGCGTGACGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGGTGGTTGTAGCTACACGAGAAACGGTCTGAGGGGCtgggt                    <  2:230315/244‑1 (MQ=255)
                                                                                                                                                                                                                           gTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCtgggt                    <  1:78134/251‑1 (MQ=255)
                                                                                                                                                                                                                             aCTTTGCTCAATGGCGATCCGATGGTCAAAGGCGGTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCtgggtgg                  <  2:103287/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGAAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCgggg           >  1:164725/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTg          >  2:43660/1‑251 (MQ=255)
                                                                                                                                                                                                                                       atgacgcgcagATGTTAAAAGGCGTTTTAACCCGTGCCGGGGCTTACTTTATCTCGGTGATAGGCCGCAGTGCTGTGGCACGTAAGCAAAACGCAAATTTCCTCGTGTCGATTCACGCTGATGCAGCACCGAACCGCAGTGCAACTGGCGCTTCAGTATGGGAGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCCGCTGGCTGGAGCAGCACGAGCAACAGTCGGAGCTGCTTGGTGGGGCggggtat        >  2:199685/11‑247 (MQ=255)
                                                                                                                                                                                                                                            gATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCG‑AGCTGGTGGGGGGGGCggggggtttggcg  >  1:291992/1‑241 (MQ=255)
                                                                                                                                                                                                                                             gtCCGATGGTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTGCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTGCGTATGGGTGCTCTCTAACCGTCGCGCCGACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGCAGCTGGTGGGTGGGGCGTGTGATGTGCTg  <  2:66131/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
GCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTG  >  NC_000913/4396525‑4397012

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑