Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,090,128 A→C 15.3% T126P (ACC→CCC)  yggI → Zn‑dependent metalloprotease‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,090,1280AC15.3% 161.0 / 2.3 76T126P (ACC→CCC) yggIZn‑dependent metalloprotease‑related protein
Reads supporting (aligned to +/- strand):  ref base A (37/22);  new base C (3/12);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 3.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAACGGCCTGGCTGGAAAGCTATGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTA  >  NC_000913/3089893‑3090349
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                             
cgGAACGTCCTGGCTGAAAAGTGATGAATGTCGCCTCAATCCCTGTGTGCTGTGGGAAAACAGGGAAGCTTTTATGGAAGAAGTGGTACCGCACGAACTGGCACATTTGGGGGTAGGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTCCGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCt                                                                                                                                                                                                                            <  1:659432/239‑1 (MQ=255)
cgGAACGGCCTGGCTGGAAAGCTATGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTCCCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCt                                                                                                                                                                                                                            <  2:573648/239‑1 (MQ=255)
cgGAACGGCCTGGCTGGAAAGCTATGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTTGAACTGCAATCCGTGCGTCGCAACAcc                                                                                                                                                                                                                             <  2:78809/238‑1 (MQ=255)
cgGAACGGCCTGGCTGGAAAGCTAAGAAATTCGCCTAAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTCCCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCt                                                                                                                                                                                                                            <  2:830334/239‑1 (MQ=255)
   aaCGGCCTGGCTGGAAAGCTATGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCTCACGAACTGGCACATTTGCTGGTATGGAAAAATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGGGGATGATGGAAAACGAGCTGGGGGTTCCGGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTc                                                                                                                                                                                                                          >  2:509697/1‑238 (MQ=255)
   aaCGGCCTGGCTGGAAAGCTATGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGGGGATGATGGGAAACGGGCTGGGGGTTCCCGCCCGTCGTACGCCTCAGTTCGAACTCAAATCCGGGGGCGGCAAAACCTTc                                                                                                                                                                                                                          >  1:420867/1‑238 (MQ=255)
        ccTGGCTGGAAAGCTATGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGTTGGAAACCGTGCAGGGTGTTCCAGCCCGTCGTACCAATAAGTTCGAACGGCAATCCGTGCGTCGCAACCCCCCCCCCTAc                                                                                                                                                                                                                    >  1:870580/1‑239 (MQ=255)
            gctggAAAGCTATGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGGGGAAGGGGATGATGGAAAACGTGCTGGGGGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGCGCCACAACCTTCCCCAActg                                                                                                                                                                                                                  >  1:558493/1‑235 (MQ=255)
              tggAAAACTATTAAAGTCGCCCTAAACCCGTTTTGCTGTTGGAAAACAGGGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTg                                                                                                                                                                                                               <  1:572270/238‑1 (MQ=255)
               ggAAAGGTTAGAAATTCCCCTCAATCCCTTTTTGCTGGTGAAAAACAGTGAAGGTGTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCa                                                                                                                                                                                                             <  1:549246/239‑1 (MQ=255)
               ggAAAGCTTTGAAATTCGCCCTAATCCCGTTTGGGTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGGACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCa                                                                                                                                                                                                             <  1:706796/239‑1 (MQ=255)
               ggAAAGCTTGGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACATGGAAGCTTTTATTGAAGAAGTGGTCCCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACAGGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCa                                                                                                                                                                                                             <  1:723492/239‑1 (MQ=255)
                 aaaGCTATGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCaa                                                                                                                                                                                                            <  2:709874/238‑1 (MQ=255)
                     cTATGAAATTCGCCTCAAACCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGcc                                                                                                                                                                                                       <  2:2904/239‑1 (MQ=255)
                       tttaaaTTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCa                                                                                                                                                                                                             <  1:105586/228‑1 (MQ=255)
                       aTGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCa                                                                                                                                                                                                             >  2:105586/1‑231 (MQ=255)
                        tGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAgg                                                                                                                                                                                                    <  2:324554/239‑1 (MQ=255)
                             ctCGCCTAAATCCCGTTTTGCTGTTGGAAAACAGTTAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCGCATTTGCTGGTAGGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCAt                                                                                                                                                                                               <  2:852774/238‑1 (MQ=255)
                               cGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGACATGCAATCCGTGCGTGGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCaccc                                                                                                                                                                                             >  2:856571/1‑236 (MQ=255)
                                    cAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGTAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGACCGCGAAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGCAAACGTGCTGGGTGTTCCCGCCCGCGGTACGCATGAGTCCGAACTGCAAGCCGGGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCAAGGAGCATCAGCtt                                                                                                                                                                                         >  2:527638/1‑238 (MQ=255)
                                     aaTCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTa                                                                                                                                                                                        >  2:50857/1‑238 (MQ=255)
                                       tCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCg                                                                                                                                                                                     <  1:840877/239‑1 (MQ=255)
                                            atttGCAGTTGGAAGACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTGTGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGAATCAGTTCGAACTGCAATCCGTGCGTCGCAACCCCTTCCCCTACCCCTGCAAGTGcccgcacactcaccttacccccac                                                                                                                                                                                 >  1:325698/2‑216 (MQ=255)
                                               tGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTAcgccg                                                                                                                                                                              >  2:790984/1‑238 (MQ=255)
                                                     tGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGGGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTCCGTCGCAAACCCTTCCCCTACCCCTCCAAGAGCCAGGAGCATcgccctcacgaccgcccccacaaa                                                                                                                                                                       >  1:707469/1‑216 (MQ=255)
                                                      ggAAAACAGTGAAGCTTTTATTGAAGACGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACTTGCTGGGCGTTCCGGCCCGTCGTACGCATCAGTTCGAACTGAAATCCGCGCGTTGCAACACCTTCCCCTACCCCTGCAAGAGCAAGGCGCCACAGCTTACCTAACACCGCACTAATc                                                                                                                                                                      >  1:527813/1‑239 (MQ=255)
                                                                      tttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGcaa                                                                                                                                                      >  2:188235/1‑237 (MQ=255)
                                                                      tttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCAGCGTAGCGCCCCATGGCAAAGAGTGGACGTGGATGATGGAAAACGTGCTGGGCGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCCGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCTCAATCGCATAGTTCGTGGGGa                                                                                                                                                      >  1:540806/1‑239 (MQ=255)
                                                                       ttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACGCCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTAGGCCGCCATAATCGCGTAGTACCTGGCGAg                                                                                                                                                     >  2:544489/1‑239 (MQ=255)
                                                                       ttATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCACCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGa                                                                                                                                                      >  2:669958/1‑238 (MQ=255)
                                                                       ttATGGAAGGAGTTGTCCCGTACAAACGGGCACTGTTGCTGGTATGGAAACTTTGCGGCTTCTTAGAGCCACGTGGCAAAGGTGGGAAGTGGATGAGGGAACACGTGTGGGGTCTTCCCTCCCGTCGTACCCCACAGTTCGAACTGCAACCCGTGCGTCGCAACACCTCCCCCTACCGCTGCAAGTGCAAGGGGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAg                                                                                                                                                     <  1:39989/239‑1 (MQ=255)
                                                                        tattgaaagggggggaCCGCACGAAGTGGCAATTGTGGGGGTAGGGAAACGTTTGGGCCGGGTTGCGCCAATGGGAAAAGGGTGGAAGTGGATGATGGAAACCGTGCGGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAACCCGTGCGCCGCAACACCTTCCCCTCCCGCGGCAAGTGCCGGGAGAATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAgg                                                                                                                                                    <  1:517293/224‑1 (MQ=255)
                                                                         aTTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGGGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCCACCGCTGCAATGGCCGGGGGCATCAGCTTACCGAACCCCGCCATAAACGCGTAGGTCGGGGTGAggg                                                                                                                                                   >  1:409030/1‑238 (MQ=255)
                                                                          tttaagtatggttcccGCCCGAACGGCTCATTTGGCTGTGATGGAAACTTTTCGTCCGTGTCGGCCCGATTGCCAAAGAGTGGAAGTTGATGGTCAAAATCGTGCTGGGTCTTCCCGCCCGGCGTCCGTATCAGTTCGACCTGCAATCCGTGCGTCGCAACTCCTTCCCCTACCGGTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGcc                                                                                                                                                  <  1:206951/225‑1 (MQ=255)
                                                                          ttGAAGAAGTGGTACCGCACGACCTGGCACATTGGCTGGTATGGAACCATTTCGGCCGCGTAGCGCCACATGGCAAAGAGGGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGcc                                                                                                                                                  <  2:420867/239‑1 (MQ=255)
                                                                               gaagTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGTTGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGGGCATAAGCCTACCGTACGCCCCCCTAAACGCGTAGTTCGTGGCAAGGCCGCCtt                                                                                                                                             >  2:490937/1‑238 (MQ=255)
                                                                                  gTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGACCTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGtctctc                                                                                                                                           >  2:346024/1‑238 (MQ=255)
                                                                                    ggTACCGAACGAACGGGCCCATTTGCTGGTATGGAAACATTTCGGCCGCGTGGCGCCAAATGGCAAAGAGTGGAAGGGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGAAATCCGTGCGTCGCAACACCTTCCCCTACCGCGGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGc                                                                                                                                         <  2:255128/238‑1 (MQ=255)
                                                                                         agggacggcTGGCAATTTTGTTGGTATGGAACCTTTTCGGCCGGGTAGCGCCACATGGCAAGGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTCCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACCCCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGtt                                                                                                                                   <  2:412245/231‑1 (MQ=255)
                                                                                           cACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCGGTGTTc                                                                                                                                  >  2:445535/1‑238 (MQ=255)
                                                                                                aCTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGCGCATCAGCTTACCGTACCCCGCCATAATCGCGACGTTCGTGGGGGGGCGGTCCATCGCtggtttccccggc                                                                                                                            >  1:491571/1‑228 (MQ=255)
                                                                                                aCTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTa                                                                                                                                                                 >  2:754336/1‑202 (MQ=255)
                                                                                                aCTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGCGTTCACTGCg                                                                                                                            >  2:693378/1‑239 (MQ=255)
                                                                                                          ttGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGCCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATc                                                                                                                                                                      >  2:596456/1‑187 (MQ=255)
                                                                                                             gCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGGAGTTCGGGGCGAGGCCGTCTATCGCGGGGTTCACTCCggggcatgtgtct                                                                                                                >  1:412245/1‑228 (MQ=255)
                                                                                                                   aTGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTTCGGTGAACAGCTGGTTGCGa                                                                                                          >  2:307342/1‑239 (MQ=255)
                                                                                                                            tttCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCCCGTAGTTCGTGGCGAGGCCGTCt                                                                                                                                              >  2:62749/1‑194 (MQ=255)
                                                                                                                                 tccctggggcccaacggggaaatgagTGGAAGGGGTTGAGGGAAAGCGTGCTGGGGTTCCCCGCCCGTTGCCCGCCTCGGTCCGAACGGCAATCCGTGCGCCCCAACCCCTCCCCGCTCCGCGGCAATGGCCGGGAGCAGCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGaa                                                                                            <  1:583071/216‑1 (MQ=255)
                                                                                                                                     cgTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTCCGCCGCCATAGTCGCGTAGTTCGTGCCAAGGCCGCTTACCCCTGTGTTCCCTCCGGGGAACAGCTGGTTGCGAAAAACCCCTCTGAACTAt                                                                                        >  1:701465/1‑239 (MQ=255)
                                                                                                                                       tAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCCTCTATCGCTGGTTTAACCCCGGGGAACAACGGGTTGCGAAAAAACAACATGAAcacaaa                                                                                      >  1:852332/1‑234 (MQ=255)
                                                                                                                                               caTGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCCCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACAATCTGAACTTTCAGGAACtt                                                                               >  2:770209/1‑238 (MQ=255)
                                                                                                                                                 gggCAAAGGTGGAGAGTGGATGTGGAAAACCGTGCGGGGTTGCCCCCCCCCTCTTCCGCCCAGTTTCGAAGCTCAATCGGTGCGTCCCAACACCTTCCCCTCCCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTCCGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                            <  1:20947/238‑1 (MQ=255)
                                                                                                                                                           gggaagtggtTGGTGGAAAACGGGCTGGGTGTTCCCCCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCGGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGc                                                                   <  2:719878/237‑1 (MQ=255)
                                                                                                                                                             gaagtggaTGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGACCTATCAGGAACTTTCCTGATCTGGCTGa                                                                >  2:67244/1‑239 (MQ=255)
                                                                                                                                                              aaggggaTGATGGAAACCGGGCTGGGTTTCCCCGCCGGTCGTCCGCATCAGTTCGAACGGCAATCCGGGCGTCGCAACCCCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGa                                                                <  2:852332/238‑1 (MQ=255)
                                                                                                                                                                   ggaTGATGGAAAACGGGCTGGGTGTCCCCGCCCGTCGTACGCAGCAGGTCGAACTGCAATCCGTGCGTCGCACCCCCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCa                                                           <  2:409030/239‑1 (MQ=255)
                                                                                                                                                                    gatgatGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTCGTTGCAt                                                          >  2:336550/1‑239 (MQ=255)
                                                                                                                                                                    gatgatGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCGGAACTATCAGGAACTTTCCTGATCTGGCTGATTGct                                                           >  2:756890/1‑237 (MQ=255)
                                                                                                                                                                         tgGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaa                                                     >  2:708065/1‑239 (MQ=255)
                                                                                                                                                                             aaaaccgtgtGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAACCCGTGCGGCGAAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaaaa                                                   <  2:149697/229‑1 (MQ=255)
                                                                                                                                                                              aaCGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAGTAACCATCTGAACTATAAGGAACTTTACTGATCTGGCTGATTGCATACCAAACCAg                                                >  2:315792/1‑239 (MQ=255)
                                                                                                                                                                              aaCGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGAGCTGGCTGGTTCAGTACCAAAACa                                                 >  2:746097/1‑238 (MQ=255)
                                                                                                                                                                                     cTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCTTACCAAAACAGCTTTc                                           >  2:559604/1‑238 (MQ=255)
                                                                                                                                                                                      tGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGACCTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCg                                          >  2:437670/1‑238 (MQ=255)
                                                                                                                                                                                        ggtttcccccccccgcgtCCGCATCGCTTCGAACCGAAATCCGGGCGCCGCAACACCTCCCCCTCCCGCTGCAAGTGCCAGGAGCATCAGCTTCCCGGCCGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAGACAGCTTTCGCTa                                       <  1:621527/224‑1 (MQ=255)
                                                                                                                                                                                                  cccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCACTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGCGGCCGTCGCTCGCTGTGTTCACTGCGGTGACCAGCATGTTGCGAAATAACCATATGAACTCTCAGGAACTTTCCTACTGTGGGTGATTGCATACCAAAACAACTTTCGCAAcagttctgg                              >  1:817906/1‑230 (MQ=255)
                                                                                                                                                                                                   ccttcctACGCATCAGTTCGAACTGGAATCCGGGCGTCGAAACCCCATCCCCTACCGCTGCAAGTGCCAGGAGAATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCCAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt                            <  2:299470/233‑1 (MQ=255)
                                                                                                                                                                                                   ccgTCGTTCGCATAAGTTCGAACTTCAATCCGTGCTTCGCAACCCCTTCCCCTACCGCTGCAAGTTCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt                            <  2:71196/239‑1 (MQ=255)
                                                                                                                                                                                                   ccctcggtccccatcgttccGAACTGCACTCCGTGCGCCGCACCCCCTTCCCCTCCCCCTGCAAGTGCCAGGGGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc                             <  2:870580/220‑1 (MQ=255)
                                                                                                                                                                                                               tCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGTCCGTTTTAACCCg                >  2:267907/1‑239 (MQ=255)
                                                                                                                                                                                                                    ttcaccctcaactCCGTGCGGCGCAACCCCTCCCCCTCCCGCGGCAAGGGCCAGGGCCACCACTTTCGCGTCCGCCCGCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGt            <  1:693009/227‑1 (MQ=255)
                                                                                                                                                                                                                    ttCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTCCGTTGCTGGCCCGTTTTAACACggcga            >  2:440576/1‑236 (MQ=255)
                                                                                                                                                                                                                    ttCAAACGGCATCCCGGGCGCCCCACCCCCTCCCCCCACCCGGGCAAGGGCCAGGAGCATCAGTTTCCCGTCCGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGt            <  1:839603/239‑1 (MQ=255)
                                                                                                                                                                                                                     tcaaacctcatgcagggCGTCGCAACCCCTTCCCCTACCGCTGCAAGTGCAAGGAGCATCAGCTTACTGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGt            <  2:734175/222‑1 (MQ=255)
                                                                                                                                                                                                                     tCGAACTGCAACCCGTGCGTCCAAACCCTTTCCCCTTCCGGTGAAGGTCCCAGGAGCTTCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGt            <  2:491571/238‑1 (MQ=255)
                                                                                                                                                                                                                              aaTCCGTGCGTCGCACCCCCTTCCCCTGCCGAGGCAATGGCCAGGAGCATCAGCTTACCGTACGCCGCCCAAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTa  <  2:896809/239‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                             
CGGAACGGCCTGGCTGGAAAGCTATGAAATTCGCCTCAATCCCGTTTTGCTGTTGGAAAACAGTGAAGCTTTTATTGAAGAAGTGGTACCGCACGAACTGGCACATTTGCTGGTATGGAAACATTTCGGCCGCGTAGCGCCACATGGCAAAGAGTGGAAGTGGATGATGGAAAACGTGCTGGGTGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTA  >  NC_000913/3089893‑3090349

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