Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,104,887 G→A S152N (AGC→AAC)  mltC → membrane‑bound lytic murein transglycosylase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,104,8870GA98.1% 165.5 / NA 53S152N (AGC→AAC) mltCmembrane‑bound lytic murein transglycosylase C
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (28/24);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATATCAGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGCCGTTCCGA  >  NC_000913/3104643‑3105118
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
cAATATCAGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACtt                                                                                                                                                                                                                                   <  1:177216/251‑1 (MQ=255)
 aaTATCAGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTc                                                                                                                                                                                                                                  <  1:229644/251‑1 (MQ=255)
  atatCAGACCCGCAGACACATCGACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGCACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGCGACGATACGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGCTTCGCTGGGAAGGTCGCGCAAACAACTTcg                                                                                                                                                                                                                                 >  2:404887/1‑251 (MQ=255)
         aCCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGAGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGGCGCGCAAGCACCTTCGCGGAtt                                                                                                                                                                                                                           >  2:311768/1‑250 (MQ=255)
               aGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATctgct                                                                                                                                                                                                                    >  1:433821/1‑251 (MQ=255)
                    caTCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGCTTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGaaa                                                                                                                                                                                                                >  2:142376/1‑250 (MQ=255)
                                           ttACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAAc                                                                                                                                                                                        >  2:27565/1‑251 (MQ=255)
                                              cTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGCGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACggg                                                                                                                                                                                     >  2:242990/1‑251 (MQ=255)
                                                     aaCCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGt                                                                                                                                                                               >  1:230255/1‑250 (MQ=255)
                                                              cGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAACCCGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACAACGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCGGCGTATCATCTAc                                                                                                                                                                      >  2:202203/1‑250 (MQ=255)
                                                               gggACAGAACCTGCCGCGCCTTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACAGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTatca                                                                                                                                                                           >  2:448502/1‑244 (MQ=255)
                                                               gggACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTatca                                                                                                                                                                           <  1:448502/244‑1 (MQ=255)
                                                                    agagccTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGGTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACGTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCa                                                                                                                                                               <  1:183249/247‑1 (MQ=255)
                                                                     gAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCAc                                                                                                                                                              >  1:421167/1‑251 (MQ=255)
                                                                       aCCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCGTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTTCGTATCATCTACAGCGTCACCa                                                                                                                                                            <  1:399386/251‑1 (MQ=255)
                                                                              gcgcATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACAt                                                                                                                                                     <  1:242990/251‑1 (MQ=255)
                                                                                 cATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTatc                                                                                                                                                                            >  1:54981/1‑225 (MQ=255)
                                                                                 cATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTGACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTatc                                                                                                                                                                            <  2:54981/225‑1 (MQ=255)
                                                                                     tGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGcc                                                                                                                                               <  1:73342/250‑1 (MQ=255)
                                                                                           gggCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAaccac                                                                                                                                        <  1:116965/251‑1 (MQ=255)
                                                                                             gCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCt                                                                                                                                      <  1:322026/251‑1 (MQ=255)
                                                                                               aaTTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGATAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCGCCTTg                                                                                                                                    <  2:433821/251‑1 (MQ=255)
                                                                                               aaTTATCAAAACGTTATGGCTGGGTGACGATCCGAGTTCGGTCGATCTCTATTCAGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACAACGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACAGGCTTCGTATCATCTAAAGCGTCACCATTAACATGGTTCCGAACCACCTTg                                                                                                                                    >  2:96708/1‑251 (MQ=255)
                                                                                                 ttATCAAAATGTTATTGAAGGGTGCCGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCACACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACCTGGTGCCGAACCACCTTGAt                                                                                                                                  >  2:85915/1‑251 (MQ=255)
                                                                                                       aaaCGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGt                                                                                                                            >  1:329557/1‑251 (MQ=255)
                                                                                                             tATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGACCCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCCCCTTGATAAACGTGCGca                                                                                                                       <  2:80739/250‑1 (MQ=255)
                                                                                                                     ggTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATAtct                                                                                                              <  1:99993/251‑1 (MQ=255)
                                                                                                                      gTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGTCAGCCGATTCGCGGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGGCGAACCACCTTGATAAACGTGCGCACAAATAtctc                                                                                                             <  2:32543/251‑1 (MQ=255)
                                                                                                                           gATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAACCGTGCGCACAAATCTCTCGGCAt                                                                                                        >  1:84570/1‑251 (MQ=255)
                                                                                                                              ccGAGTTCGGTCGATCGCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGGGCGCACAAATAGCTCTGCATgg                                                                                                      <  2:309250/250‑1 (MQ=255)
                                                                                                                               cGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTc                                                                                                    <  1:27565/251‑1 (MQ=255)
                                                                                                                               cGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTc                                                                                                    >  1:351921/1‑251 (MQ=255)
                                                                                                                                         cGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGt                                                                                          >  1:56472/1‑251 (MQ=255)
                                                                                                                                              tctATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCt                                                                                                                                                                                                        >  1:313832/1‑136 (MQ=255)
                                                                                                                                              tctATGCCGACGTTGATGAGATTACGATTTCGAAAGAACCTTTCCGTGACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCt                                                                                                                                                                                                        <  2:313832/136‑1 (MQ=255)
                                                                                                                                               ctATTCCGACGTTGATGATATTACGATTTCGAAATAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCAc                                                                                       <  2:329557/248‑1 (MQ=255)
                                                                                                                                                tATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGa                                                                                    >  1:54318/1‑250 (MQ=255)
                                                                                                                                                      gACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGACCGCCAGGCGTCACGGAAATATgg                                                                             >  1:424692/1‑251 (MQ=255)
                                                                                                                                                        cGTTGATTATATTACTATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCTCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTTAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCg                                                                           <  1:167513/251‑1 (MQ=255)
                                                                                                                                                             atgatATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGCACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCCGGCGTCACGGAAATATGGCGTtgat                                                                      >  2:107522/1‑251 (MQ=255)
                                                                                                                                                                  atTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCTTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATACCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTc                                                                 <  1:280899/251‑1 (MQ=255)
                                                                                                                                                                   ttACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGCTTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACAGAAATATGGCGTTGATGAGTCg                                                                >  2:132912/1‑251 (MQ=255)
                                                                                                                                                                      cGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTg                                                             >  1:246775/1‑251 (MQ=255)
                                                                                                                                                                          ttCGAAAGATCCTTTCCTGTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGAtt                                                          <  2:19423/250‑1 (MQ=255)
                                                                                                                                                                           tCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAACAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCt                                                        >  2:58805/1‑251 (MQ=255)
                                                                                                                                                                                  aaCCTTTCCTTTACGGTCAGGTGGTGGACAACCCCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAAtt                                                 <  1:58805/251‑1 (MQ=255)
                                                                                                                                                                                        tcctttACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTAGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGACGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGc                                             >  1:421590/1‑249 (MQ=255)
                                                                                                                                                                                            ttACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCGCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGAc                                         <  1:202203/249‑1 (MQ=255)
                                                                                                                                                                                                       tggtggACAACACCGGGCAGCCGATTCGCTGGGCAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCttt                            >  1:158634/1‑251 (MQ=255)
                                                                                                                                                                                                        ggtggACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTa                           <  2:246775/251‑1 (MQ=255)
                                                                                                                                                                                                               aacaCCGGTCAGCCGATTCGCTGGGAAGGTCGCGCACACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACCGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACCAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGACTCTTCCTTTAACCCGt                     >  2:94795/1‑250 (MQ=255)
                                                                                                                                                                                                                    cGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAACAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCgg               >  2:255134/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gCTGGGAAGGTCGCGCCAACAACTTCGCGGATTATCTGCTGACAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAACGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAACTATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTCTGCAGACCGAATCTTCCTTGAACCCGTATGCGGTCAGCCGTTCCGa  >  2:182710/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
CAATATCAGACCCGCAGCCACATCAACTTCGATGACGGTACGATTACTATCGAAACCATCGCCGGGACAGAACCTGCCGCGCATTTGCGCCGGGCAATTATCAAAACGTTATTGATGGGTGACGATCCGAGTTCGGTCGATCTCTATTCCGACGTTGATGATATTACGATTTCGAAAGAACCTTTCCTTTACGGTCAGGTGGTGGACAACACCGGGCAGCCGATTCGCTGGGAAGGTCGCGCAAGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGCAGACCGAATCTTCCTTTAACCCGTATGCGGTCAGCCGTTCCGA  >  NC_000913/3104643‑3105118

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