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,7720TC98.2% 194.7 / ‑7.2 57F236F (TTT→TTC) amiBN‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (27/29);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCG  >  NC_000913/4396533‑4397015
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
gcgcgcAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCCGGCGCTATCGGCCCCGGGGGTACGCGGGGGAAAAATATCACCATCGCCATCGCGCGGAAATTGCGTACTTTGCTCAATGGCGATCCGCTGTTCAAAAGCGtttt                                                                                                                                                                                                                                          >  2:335788/1‑251 (MQ=255)
        cccGTTTAAAACGGAAAGTAACCGCACTACGGGTTTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTg                                                                                                                                                                                                                                  <  2:380471/250‑1 (MQ=255)
         cagTTTAAAACGGAAAGTAACCGCCCTACTGGTGTTATCAGCAGGAATACGGGAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTCCTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGa                                                                                                                                                                                                                                 <  2:224694/249‑1 (MQ=255)
           gTTTAAAACGGAAAGTAACCGCACGACGGGTGTGATCAGCAGTAAGACGGGAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATGATCATCGCTATTGATGCCGGACACGGCGGGCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACg                                                                                                                                                                                                                               <  2:230835/251‑1 (MQ=255)
                  aCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGGCCCTGGCGCTATCGGCCCCGGTGGGACGCGGGGGAAAAATGGCACCATCGCCATCGCGCGCAAAATGCGTACTTTGCTCAATTACGCTCCGATGTTCAAAGGCGTTTTAACCCGTg                                                                                                                                                                                                                                  >  2:369696/1‑241 (MQ=255)
                  aCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTg                                                                                                                                                                                                                                  <  1:369696/241‑1 (MQ=255)
                   cGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTAc                                                                                                                                                                                                                       >  1:147619/1‑251 (MQ=255)
                                    aCGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATggg                                                                                                                                                                                                      <  1:265352/251‑1 (MQ=255)
                                    aCGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATggg                                                                                                                                                                                                      <  1:179697/251‑1 (MQ=255)
                                          ggtaacagcatTAATACGGTAACGCGTCCGGCAGCGGCCGGGACGGCTAACCCTGGGGTGAAAATTATCATCGCTATGGATGCCGGACACGGCGGTCAGGACCCGGGCGCTAGCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCAGCGCCATCTCGCGTAAATTGCGTACTTTGCTCAATGTCTAGCCGATGTTCAAAGGCGGTTTAACCCTTGTCGGGTATTACTTTATCTCGGGGATGGTGCGCAg                                                                                                                                                                                                <  2:393592/246‑1 (MQ=255)
                                             aTCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTg                                                                                                                                                                                              <  1:132754/250‑1 (MQ=255)
                                                agcagTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGAtgt                                                                                                                                                                                          >  1:341157/1‑251 (MQ=255)
                                                 gcagTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCCGGACCCTGGCGCTATCGGCCACGGTGGGACGCGGGAGAAAAATGTCACCATCGCCATCGGGCGTAAATTGCGTACTTTGCTCAAATACGATCCCATGTTCAAAGGCGTTTTAACCcggggcggggg                                                                                                                                                                                                                           >  1:173890/1‑216 (MQ=255)
                                                 gcagTAATACGGTAACGCGGCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAAGGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCGCAATGACGATCCGATGTTCAAAGGCGTGTTATCCCGTGACGGGGa                                                                                                                                                                                                                           <  2:173890/217‑1 (MQ=255)
                                                        tACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCg                                                                                                                                                                                                                                              <  1:217753/191‑1 (MQ=255)
                                                        tACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCATAGCGCGTAAATTGCGTACTTTGCTCAATGACGGTCCGATGATTAAAGGCg                                                                                                                                                                                                                                              >  2:217753/1‑191 (MQ=255)
                                                                 gcgTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACg                                                                                                                                                                         <  1:335788/251‑1 (MQ=255)
                                                                  cgTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGGACCTGGCGCTATCGGCCCCGGTGGTACGCGGGGGAAAAATGTCACCCTCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACCATCCGCTGTTCAAAGGCGGTTTAACCCGTGACGGGGGTTACTTTATCTCGGTGATGGGGCGCAGGGGTGTGGGACGGAAGCaaaaagc                                                                                                                                                                        >  2:146689/1‑248 (MQ=255)
                                                                  cgTCAGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGACACGGTGGTACGCGGGAGAAAAATGTCACCATCGCCCTCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGCTGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCAAGTAAGCaaacagc                                                                                                                                                                        >  1:46154/1‑247 (MQ=255)
                                                                      cGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCGCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAAt                                                                                                                                                                    <  1:29237/251‑1 (MQ=255)
                                                                          aGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGACCCGGTGGTACGCGGGCGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATACGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTACGCAAAACGCCAATTTcc                                                                                                                                                                >  1:415105/1‑251 (MQ=255)
                                                                            cgcgcgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGGACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTc                                                                                                                                                              >  1:270180/1‑251 (MQ=255)
                                                                              cgcgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCgt                                                                                                                                                            >  1:114253/1‑251 (MQ=255)
                                                                              cgcgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCAACGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGGTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATTTCGGGGATGGGGGGCAGTGATGGGGCACGTAAGCAAAACGCCAATTTccccgg                                                                                                                                                            >  1:282032/1‑247 (MQ=255)
                                                                                cgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCgtgt                                                                                                                                                          <  1:140234/251‑1 (MQ=255)
                                                                                        cTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGGGAAAAATGTCACCATCGCCATCGCGCGGAAATTGCGTACTTTGCTCAATGACGATCCGGTGTTCAAAGGCGTTTTAAACCGTGTCGGGGATTACTTTTTCTCGGGGTTGGGGCGCAGAGATGTGGCACGTAAGCAAAACGCCAATTTTCTCGTGTCGGTTccc                                                                                                                                                  >  2:336598/1‑249 (MQ=255)
                                                                                           acacTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCTGTGGTACGCGGGAGAAAAATGTCACCAGCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCt                                                                                                                                               <  2:361019/251‑1 (MQ=255)
                                                                                                  ggaaaaaaTTATCACCGCTATTGATGCCTGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCg                                                                                                                                        <  2:84502/247‑1 (MQ=255)
                                                                                                   gATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGGTTTAACCCGTGACGGGGGTTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACACTGATGCCGc                                                                                                                                       >  1:320425/1‑251 (MQ=255)
                                                                                                    aTAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGcc                                                                                                                                                                       <  2:351264/218‑1 (MQ=255)
                                                                                                    aTAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGGGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGcc                                                                                                                                                                       >  1:351264/1‑218 (MQ=255)
                                                                                                                               gACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTcc                                                                                                           <  2:114253/251‑1 (MQ=255)
                                                                                                                                  acGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGGACGCGGGAGAAAAATGTCACCCTAGCCATCCCGCGTAAATTGCATACTTTGCTCCATGACGGTCCGACGTTCAAAGGCGTTTTAACCCGTGACGGGGGTTACTTTATCTCGGTGATGGGGCGCCGTGGTGGGGCACGTAAGCAAAACGCCAATTTCCTCGGGTTCATTCACCCTGGTGCCCCACCGGACCGCAGCGCGACCGGCGCTTCCGtt                                                                                                        >  1:234160/1‑250 (MQ=255)
                                                                                                                                   cggcggTCAGGACCCTGGCGCTATCGGCCCCGTGGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTAt                                                                                                       <  1:344481/251‑1 (MQ=255)
                                                                                                                                          cAGGACCCGGGCGTTATCGGCCCCGGTGGTACGCGGGAGAAAAAGGTCACCATCGCCATCGCGCGTAAATTGCGTACGTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGc                                                                                                 <  2:18672/250‑1 (MQ=255)
                                                                                                                                          cAGGAACCTGGCGCTATCGGCCACGGGGGTACGCGGGAGAAAAATGGCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGACGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCCGTGATGAGGCACGTAAGCAAAACGCCACTTTCCTCGTGTCGATTAACGCTGATGCCGCACCGAACAGCAGCGCGACTGGCGCTTCCggagtggggcc                                                                                                >  2:371120/1‑241 (MQ=255)
                                                                                                                                           aGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAACCGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGctc                                                                                               >  1:368551/1‑251 (MQ=255)
                                                                                                                                            ggACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTAACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGctct                                                                                              <  1:228915/251‑1 (MQ=255)
                                                                                                                                                      gcTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTcg                                                                                    <  1:126608/251‑1 (MQ=255)
                                                                                                                                                            ggCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTACCCGTCGCGCCac                                                                              >  2:250422/1‑250 (MQ=255)
                                                                                                                                                               cccGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAg                                                                           <  1:336598/251‑1 (MQ=255)
                                                                                                                                                                  ggtggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGACGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTg                                                                         >  1:164977/1‑250 (MQ=255)
                                                                                                                                                                          gcgGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGccc                                                                >  2:101407/1‑250 (MQ=255)
                                                                                                                                                                                  aaaaTGTCCCCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctggctg                                                        <  2:341157/251‑1 (MQ=255)
                                                                                                                                                                                  aaaaTGTCACCATCGCCAACGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctggctg                                                        <  1:161500/251‑1 (MQ=255)
                                                                                                                                                                                          aCCATCGCCAGCGCGCGGAAATGTCGTACTTTGCTCATTGACGACCCGATGTTCAAAGGCGTTGTAACCCGTGACGGTGATTACTTTAGCTCGGTGATGGGGCGCAGTGTAGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAgcagca                                                <  2:228997/251‑1 (MQ=255)
                                                                                                                                                                                           ccatcgccatcgCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCGCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCAc                                               <  1:300864/251‑1 (MQ=255)
                                                                                                                                                                                             atcgccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACg                                              >  1:220570/1‑250 (MQ=255)
                                                                                                                                                                                                gccatcgcGGGGAAATTGCGTACTTTGCTAAATGACGATACGCTGTACAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGAGATGGGGCGCTGTGATGTGGCACGGACGAAAAACGCCAATTTACTCGAGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGACGCTTCCGTATGGGTGCTCTATAAACGTCACGCCAACAGGGAAATGGCCAtctggctagggcccaccggaa                                           >  2:64074/1‑236 (MQ=255)
                                                                                                                                                                                                   atcgcGCGTAAATTGCGTACTTTGCTCAATGACGCTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAAc                                        >  1:127786/1‑250 (MQ=255)
                                                                                                                                                                                                      ccgcgTAAATGGCGTACTTTGCTCAATGAAGAGCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGAGTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTc                                    <  1:360260/250‑1 (MQ=255)
                                                                                                                                                                                                       cgcgTAACTTGAGTACTTTGCTCAATGACGCTCCGATGTTCAAAGGCGGTTTAAACCGTGACGGGGCTTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGTTTCACGCTGCGGCCGCCCCGAACCGCAGTGCGACTGGCGCGTCCGTATGGGGTCTCTCTACCCCTCGCGCCCACAGTGAAATGGCCAGCTGTCTGGagacgcgagaggaacagtc                                    >  2:360260/1‑233 (MQ=255)
                                                                                                                                                                                                           tAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCGGGCTGGAGCAGCACGAGAAACAGTCGGAgc                               >  2:349212/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTggggg                  >  2:173303/1‑250 (MQ=255)
                                                                                                                                                                                                                                     aTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGgcg                                                                                                                >  2:274113/1‑144 (MQ=255)
                                                                                                                                                                                                                                     aTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGgcg                                                                                                                <  1:274113/144‑1 (MQ=255)
                                                                                                                                                                                                                                        cGATGTTCAAAGGCGTTTTAACCCGGGACGGGGATTACGGTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCg  <  2:40880/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
GCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCG  >  NC_000913/4396533‑4397015

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