Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,396,373 C→T T103T (ACC→ACT)  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,3730CT96.2% 167.5 / ‑3.7 55T103T (ACC→ACT) amiBN‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (35/16);  total (36/19)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

CGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGAC  >  NC_000913/4396135‑4396618
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
cGCTCTCTGATATTCAGGTTTCTAATGGTAATCACCAGGCGCGGATAACGTTGAGTTTTCTTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTaa                                                                                                                                                                                                                                            >  2:274815/1‑250 (MQ=255)
 gCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTaaaa                                                                                                                                                                                                                                          <  1:399588/251‑1 (MQ=255)
     tctGATATTCAGGTTTCTACCGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGCTTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACCGCCGTTGTTGTTAAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGAGGTCGATTTTACTGAAAACGGTAAAACCg                                                                                                                                                                                                                                       >  2:272326/1‑250 (MQ=255)
                              aTCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCGCCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGgcagc                                                                                                                                                                                                             <  1:303941/251‑1 (MQ=255)
                                cAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGgcagca                                                                                                                                                                                                            <  2:280015/250‑1 (MQ=255)
                                  aCAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTGCAGCGGCAATAATCTGGTGAAGGCGGTTCGCTGTGGAACGGCTAAAGATGCACAAACGGTGCGGCTGGTGGTCGATCTTACTGAAAGCGGTAATGCCGAAGCGGTGAAGCGGCAGAATGGCAGCAAtt                                                                                                                                                                                                         <  1:180156/251‑1 (MQ=255)
                                        gcgGATACCGTTGAGTTGTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTg                                                                                                                                                                                                   <  2:355378/251‑1 (MQ=255)
                                          ggATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTgta                                                                                                                                                                                                 <  2:92797/251‑2 (MQ=255)
                                          ggATAACGGTGAGTTTTATTGGCGATCCTGATTATGCGTTTCGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGAGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGGCGATCTTACTGCAAACGGTAAAACCGAAGCGGTGAAGCGGCAGACTGGCCGCAATTAAACTgtc                                                                                                                                                                                                 >  2:283797/1‑251 (MQ=255)
                                            aTAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTgtcgt                                                                                                                                                                                               <  1:281720/251‑1 (MQ=255)
                                            aTAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTAAACTgtcgt                                                                                                                                                                                               >  1:392417/1‑251 (MQ=255)
                                                cGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGCCTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTACAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGTCttt                                                                                                                                                                                           >  2:189004/1‑251 (MQ=255)
                                                cGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCttt                                                                                                                                                                                           >  1:407934/1‑251 (MQ=255)
                                                 gTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAACCCGAAGCGGTGACGCGGCAGAATGGCAGCAATTACACTGTCGACtttc                                                                                                                                                                                          >  2:431442/1‑250 (MQ=255)
                                                 gTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGTCAGCAATTACACTGTCGTCttt                                                                                                                                                                                           >  2:104671/1‑250 (MQ=255)
                                                       ttttATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGAtt                                                                                                                                                                                     >  1:180571/1‑250 (MQ=255)
                                                        tttATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAATCCGAAGCGGTGAAGCGGCAGAAGGGCAGCAATTACACTGTCGTCTTTTCGATTa                                                                                                                                                                                    <  2:380480/250‑1 (MQ=255)
                                                          tATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAACGGTGAAGCGGCAGAATGGCAGCACTTACACTGTCGCCTTTACGATTAACg                                                                                                                                                                                 >  2:3456/1‑251 (MQ=255)
                                                           aTTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCCGAATGGCAGCAATTACACTGTCGGCTTTACGATTAACGc                                                                                                                                                                                >  2:167666/1‑251 (MQ=255)
                                                                 gATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGAtgt                                                                                                                                                                          >  1:215111/1‑251 (MQ=255)
                                                                           aTGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTgccgcccccc                                                                                                                                                                >  2:287988/1‑247 (MQ=255)
                                                                           aTGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCg                                                                                                                                                                >  1:199297/1‑251 (MQ=255)
                                                                              cGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGcct                                                                                                                                                             >  2:130727/1‑251 (MQ=255)
                                                                               gTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGCTATCAAACAAACGGGCGTGATTCAGGGCCTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGAGTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGGAAAACCGCAGCGGCGCAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGACGCCACCgacgc                                                                                                                                                            >  2:101495/1‑247 (MQ=255)
                                                                                  tAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAACGGGGAAGCGGCAGAATGGCAGAAATTACACTGTCGTCTTTACGATTACCGCCGATGTGCCGCCACCGCctcctc                                                                                                                                                         >  2:96913/1‑251 (MQ=255)
                                                                                     ccATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTccgc                                                                                                                                                      >  1:163130/1‑251 (MQ=255)
                                                                                      cATCAAAGCAAACGCACCGTGGCGCGCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTccgcc                                                                                                                                                     <  2:135011/251‑1 (MQ=255)
                                                                                         cAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCgccgcc                                                                                                                                                  >  1:359480/1‑251 (MQ=255)
                                                                                                    cACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCCGCGGCAATAATCTGGTGAAGGCGATTCGATCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCg                                                                                                                                        >  1:179567/1‑250 (MQ=255)
                                                                                                        gTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAAcg                                                                                                                                   >  1:50765/1‑251 (MQ=255)
                                                                                                                        aacaaacGGGCGTGATTCAGGGACTGCCGTTGTGGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAGCCGAAGCGGTGAAGCGGCAGAAGGGCAGCAATTACACTGGCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCGCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCg                                                                                                                   <  2:215111/251‑1 (MQ=255)
                                                                                                                            aacGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGGGGTCGATCTTACTGAAAACGGAAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGACGTCTTTACGATTAACGCCGATGTGCCGCCGCCGCCTCATCCGCCGCCCGAGGTTGCGAAACGCGTTGAAACGCCTGCggatg                                                                                                               >  2:223522/1‑248 (MQ=255)
                                                                                                                               gggCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCg                                                                                                            <  1:58649/251‑1 (MQ=255)
                                                                                                                                     gATTCAGGGACTGACGTTGTTGTAAAGCGGCAATAATCTGGCGAAGGCGATCAGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGGGGTCGTTCTTACTGAAAACGGTAAAACAAAAGCGGTGAAGCGCCAGAATGGAAGCACTTACACTGTCGACTTTACGATTCACGCCGACGTGCAGCAACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCATGAGGTTGTCGCCCcgc                                                                                                      >  2:26339/1‑251 (MQ=255)
                                                                                                                                           gggACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGATCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACcgcgcggc                                                                                                 >  1:43909/1‑248 (MQ=255)
                                                                                                                                             gacggcggttgttgttCCGCGGCAATAATCGGGTGAAGGCGATTCTCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACTCTGTCGTCTTTACGATTAACGCCGATGTTCCGCCACCGCCTCCTCCGCCGCCCGTGGGGGCGAAACGCGTTGAAACGCCTTCGGTTGTCGCACCGCGCGTCAGc                                                                                              <  2:277791/244‑1 (MQ=255)
                                                                                                                                                      tgtgtttCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGTTTCACAAACGTTACGGCTGGTGGTCGTTCTTACTGAAAACGGGAAAACCGAAGCGGTGAAGGGGCAGAAGGGCAGCAATTACACTGTGGCCTTTACGATTAACGCCGATGTGCCGCCACCGCGTCCTCCGTCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCTCGGCAGCGAACCGgc                                                                                      <  2:22010/245‑1 (MQ=255)
                                                                                                                                                                                 aGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAAc                                                          <  1:189004/251‑1 (MQ=255)
                                                                                                                                                                                   gCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCg                                                        >  1:309779/1‑251 (MQ=255)
                                                                                                                                                                                     gATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGc                                                       >  1:226493/1‑250 (MQ=255)
                                                                                                                                                                                               ggAACGCCTAAAGATGCCCCAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGTCTTTACGATTAACGCCGATGTGCCGCCTCCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGgt                                              <  2:305281/249‑1 (MQ=255)
                                                                                                                                                                                                taaCGCCTAAAGATGCACAAACGCTACGGCTGGTGGGCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTTCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGACCCGCGCGTCAGCGAACCGGAGCGCAATCCGTTTAAAACGGACAGTAACCGCACGACGGgttat                                           >  2:5890/2‑248 (MQ=255)
                                                                                                                                                                                                  aCGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTAt                                         >  1:209371/1‑251 (MQ=255)
                                                                                                                                                                                                   cGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAAGGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAAGCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATc                                        <  2:163130/251‑1 (MQ=255)
                                                                                                                                                                                                          aGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGGTTAAAACGGAAAGTAAACGCACTACGGGGGTTATCAGCAGt                                  >  2:434455/1‑250 (MQ=255)
                                                                                                                                                                                                           gaggCACAAACGACACGGCTGGTGGTCGATCTTACGGAAAACGGTAACACCGAAGCGGTGAAGCGGCAGAATGGAAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTGGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCTCAATACGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTaa                                <  2:115754/248‑1 (MQ=255)
                                                                                                                                                                                                           gATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTa                                 >  2:297496/1‑250 (MQ=255)
                                                                                                                                                                                                                acaAACGCTCCGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGCTTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGCAACGCCTGCGGTTGTCGCCCCGCGCGTCAGCGAACCGCCGCGCAAGCCGCTTAAAACGGAAAGTAACCGCACCACGGGTGTTATCAGCCGTaaaccgg                           >  2:199323/1‑246 (MQ=255)
                                                                                                                                                                                                                        tACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGt                   >  1:100110/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGGAACCGCACTACGGGTGTTATCAGCAGGAATACGGTAACGCGTCCg                <  2:180571/251‑1 (MQ=255)
                                                                                                                                                                                                                            gcgggtggtCGGTATTACAGAAAACGGTAAAACCGAGGCGGTGAAGCGGCTGAATGGCAGCAATTACCCTGTCTTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAAcgcgacctg               <  2:413923/248‑6 (MQ=255)
                                                                                                                                                                                                                              tggtggTCGATCTTCCTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGAGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCCATCCGTTTGAAACGGAAAGTAACCGCACTACGGGGGTTCTCAGCCGTAATCCGGTAACGCGTCCGGc              >  2:341515/1‑250 (MQ=255)
                                                                                                                                                                                                                                    tCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAgcgcgc       >  1:196322/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAgcgcgcgc     >  1:233185/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cTTACTGAAAACGGAAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGACGACTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTCACGCGACCGGCAGCGCGCGCGAc  >  2:42645/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
CGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGAC  >  NC_000913/4396135‑4396618

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