Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,305 C→T 53.3% R395C (CGT→TGT)  icd → isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,3050CT53.3% ‑2.4 / 66.4 61R395C (CGT→TGT) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (17/11);  new base T (15/15);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

TCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGT  >  NC_000913/1196087‑1196530
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
tccccccTGGTGCAAACAGCGGTGAAGAATGCGCCCTGTTTGAAGCCACCCACGGTAATGCCCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgga                                                                                                                                                                                                                         <  2:173239/226‑1 (MQ=255)
 cgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                                                                                                                                                                                         <  1:673415/228‑1 (MQ=255)
  gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTAGAGTGTCTgctgg                                                                                                                                                                                                                          >  1:220528/1‑226 (MQ=255)
          gtcAAACCTCGGTGACGAATGCGCCCTGGTGTAAGCCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTaaa                                                                                                                                                                                                               <  1:338850/227‑1 (MQ=255)
                   cGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCCATCAACGCGAAAACCGTAAACTATGACTTCGAGTGTCGGATGGTGGGCGCTAAAACGCGGaaa                                                                                                                                                                                                      >  2:636194/1‑229 (MQ=255)
                    ggTGACGAATGCGCCCTGTTTGAAACCACCCGCGGTACTGCGCCGAAATATGCCGGTCCGGACAACGTAACACCTGGCTCTATTCTTCTCTCCGCTGAGCTGATGCTGCGCCACATGGGTTGGACCGCACCGGCTGACTTAGTTGTTAAAGTGATGGCAGGGGCAATCAACGCGAAAACCACAACCTATGACTCCGAGCGTGTTAGGGATGGTGCTAAACAGCTGaaaa                                                                                                                                                                                                     >  2:497512/1‑228 (MQ=255)
                           aaTGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGCCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGACATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAAACCCTGAAATGTTCaga                                                                                                                                                                                              >  1:401861/1‑229 (MQ=255)
                                      tttGAAGCCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                                                                                                                                                                                   <  2:446582/229‑1 (MQ=255)
                                                       ttcTGCGACTAAATTTGCCGGTTAGGACAAAGTAAATCCTGGCTGTATTATTGTCTCCTCTGGGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                                                                                                                                  <  2:220528/227‑1 (MQ=21)
                                                       ttcTGCACCGAAAAAGGGCGGTCGGTACAAAGTAACACCTGGCGCTATTATGCCCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                                                                                                                                  <  1:310670/227‑1 (MQ=18)
                                                       tACTGCCGCGAAATATGCCGGTCAGGACAACGGAAATCCTGGCTCTATTATTCGTTCCGCTGAGATGATGCTGAGTCCCATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGTCGCAATCAACGCGTAAACCGTAACCTATGACTTCGAGCGACTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                                                                                                                                  <  1:351727/229‑1 (MQ=255)
                                                        aCTGCGCGTAAATAGGCCTGTCAGGACAAAGTTAATCCTGGCTCTATTATTCTCTCCGGTGAGATGATGCCGCGCCACATGGGTTGGACCGAAGCGGCTGGCTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                                                                                                                                  <  1:648144/228‑1 (MQ=21)
                                                         cTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTa                                                                                                                                                                <  1:250968/229‑1 (MQ=255)
                                                           gcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTa                                                                                                                                                                <  1:415965/227‑1 (MQ=255)
                                                           gcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGAGGGCGCTAAACGGCTGCAAAGGTCAGAGTTTGGTAGCGCGACCATCGAAAACATGaaaa                                                                                                                                                              >  2:415965/1‑228 (MQ=25)
                                                            cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTCATTGTTAAAGGAATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGAAGGGCGCTAAAATGCTGAAATGTTCAACGTTTGGTGACGCGATCATCGCAAACAGGTAATg                                                                                                                                                             >  1:356737/1‑229 (MQ=255)
                                                            cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGGAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTCATg                                                                                                                                                             >  2:244096/1‑229 (MQ=255)
                                                                  aaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTa                                                                                                                                                        <  1:40540/228‑1 (MQ=255)
                                                                   atatGCCGGTCAGGGCAAAGGAAATCGTCGGCATATTATTCTCTCCGCGGAGATGAGGCTGTGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGATAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGt                                                                                                                                                      <  2:386735/229‑1 (MQ=21)
                                                                        ccGGTCAGGACAAAGTAAATCCTGGCTTTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt                                                                                                                                                 <  2:673497/229‑1 (MQ=255)
                                                                                     aGTAAATCCGGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCGGa                                                                                                                                                                                                        <  1:720295/161‑1 (MQ=255)
                                                                                      gTAAATCCTGGCTCTATTTTTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACgg                                                                                                                                   <  1:703013/229‑1 (MQ=255)
                                                                                       tAAATCCGGGCCCGTTTATCCGTTCCGCTGAGATGATGCTGCGCCACATGGGTCGGACCGAAGCGGCTTACTTAATTGTTAAAGTTATGGAAGGCGCAAACATCCCGAAAGCCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGGTTGTTAAATTTATTAACgg                                                                                                                                   <  1:553635/228‑1 (MQ=255)
                                                                                          aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGGGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGACCATCGAAAACATGTAATCCCTTAGTTTGTTAAATTTTTTTACGGGAg                                                                                                                                >  2:670337/1‑228 (MQ=255)
                                                                                          aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGGGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAg                                                                                                                                >  1:58568/1‑228 (MQ=255)
                                                                                               ggCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGGTCTTCGAAAACATGTAATGCCGAGGTTTGTTAAATTTTTTAACGGGAGGGTaaa                                                                                                                          >  2:703012/1‑228 (MQ=255)
                                                                                                           ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAACTGTTCAGAGTTTGGTGACTCGATCATCGAAAACATGTAATGACGGAGTTTGTTAAATTTATTAAAGGGAGCGTAACGGTCCCgttgt                                                                                                               >  1:27995/1‑228 (MQ=255)
                                                                                                              tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGATGTTGGTGACGCGATCATCGAAACAATGTAAGGCCGTAGTTTGTTAAATTTTTTAACGGGGGCATAAACGTCCCGtttttttt                                                                                                           >  2:709748/1‑229 (MQ=21)
                                                                                                              tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGTAGCGTAACGCACCCgttgt                                                                                                               >  1:753058/1‑225 (MQ=255)
                                                                                                              tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGtttt                                                                                                            >  1:594498/1‑228 (MQ=255)
                                                                                                                cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACgg                                                                                                                                   >  1:655870/1‑203 (MQ=255)
                                                                                                                cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTAACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTACATTTATTGACGGGAGCGTAACGCTCCCGTTGTTTTTTg                                                                                                         >  1:745993/1‑229 (MQ=255)
                                                                                                                 gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCTAGCGTCTGATGGATGGCGCTAAACTGATGAAATGTTCAGAGTTTGTTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGt                                                                                                        <  2:695096/229‑1 (MQ=255)
                                                                                                                     gggtaaagctgcGCCACTGGGGTTTGACCGAAGCGGCTGACTTAATTGATAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAGGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAGCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg                                                                                                    <  1:234804/222‑1 (MQ=255)
                                                                                                                     ggaggaggctggGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg                                                                                                    <  2:747387/222‑1 (MQ=255)
                                                                                                                          atgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctgct                                                                                               >  1:660519/1‑229 (MQ=255)
                                                                                                                           tgctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTa                                                                                              <  1:244096/229‑1 (MQ=255)
                                                                                                                                cgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTGACCTATGACTTCGAGCGTCTGATGGATGGCCCTAAACTGCTAAAATGTTCAGAGTTTGGTGACGCGATGATCGTAAACATTCAATGCCGTGGTTTGTTACATTTATTCACGGGGGCGGAACGCTCCCGCGGTTTTTTGTTAGGCTGCTAACGGt                                                                                         >  2:92542/1‑229 (MQ=255)
                                                                                                                                     caTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGATCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTGGTTAAATTTATTAACGGGAGCGTAACGATCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCa                                                                                    >  1:585727/1‑229 (MQ=255)
                                                                                                                                             ggACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACACGAAAACCGTAACCTATGGGTTCGAGTGTCTGATGGATGGCGCTAACCTGCTGCCATGTTCAGAGTTTGGTGAAGCGATCATCGAAACCATGTACTGCGGTAGTTTGTTAAATTTATTAACTGGCGCGTAACGCTACCGTTGTTTTTCTTTTGTCTGCTAACGGTTATCAAAATTTTa                                                                            >  1:558311/1‑229 (MQ=21)
                                                                                                                                                   aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGACAAACGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTGTGGTGACGCGATCATCGCAAACATGTAATTCCGTAGTTTGTTTAATTTATTAACGGGACCGTAAAGCTCCCGTTTTTTTTTTTTAGGTTTCTAACGGTTTACAAAATTTTTTCaaaa                                                                      >  2:404634/1‑229 (MQ=21)
                                                                                                                                                   aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGGGCGTCTCATGGCTGGCGCTAAACTCCTGAAATTTTCCGAGTTTCGTCACTCGATCGTCGCAAACATGTAATCCCGATGTTTGTTACATTTATTAACGGATGAGTAACGCTACCGTGGTTTTTTGTTATGCTGCTAACGTGTTTCAAAAATTTATCtaaa                                                                      >  2:604275/1‑229 (MQ=21)
                                                                                                                                                   aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCCGATGGATGGCGCTACACTGCCGCCATGTTCAGACTTTGTGTACGCCATCATTGCAAACATGTAATGCCGTAGTCTGTTATATTTATTAACAGGAAACGAACAGTCACGTTGTTTTTTGTTAGTCACCTACCGTTTTTCAAAATTTTacacaaa                                                                      >  2:302361/1‑223 (MQ=16)
                                                                                                                                                   aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGGGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGCGTTTGGTTACGCGATCATGGAAAACCTGTAATGCCGTAGTTTGTTCAATTTATTAACGGGAGCGTAACGCTCACGTTGTTTTTTGTTAGGCTTCTAAGGGTTCTCAAAATTTTAACacaa                                                                      >  2:613672/1‑229 (MQ=25)
                                                                                                                                                    aGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTTTTAGCCTGCTAACGGTTTTCAAAATTTTATCaaaaa                                                                     >  1:176085/1‑229 (MQ=255)
                                                                                                                                                        gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGTGTTTGGTGACGCGATCATCGAAAACATGTAATGCAGTAGTTTGTTAAATTTATTAAAGGGGGCGTAACGCTCCCGTGGTTTTTTGTTAGGCTGCTAACGGTAAACAAAATTTTTTCaaaaaatg                                                                  >  2:590208/1‑226 (MQ=255)
                                                                                                                                                        gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACgg                                                                                                                                   >  2:444531/1‑163 (MQ=255)
                                                                                                                                                        gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACgg                                                                                                                                   <  1:444531/163‑1 (MQ=255)
                                                                                                                                                                      tAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTc                                                   >  1:432161/1‑229 (MQ=255)
                                                                                                                                                                                  aGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGATGGTCTGATGGATGGCGATAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg                                       <  1:351459/229‑1 (MQ=255)
                                                                                                                                                                                    gcgcAATCAACTCGAAAACCGTAACCTATGCATTCGAGCGTCTGATGGATGGCGCTAAACGGCTGAAATGTTCAGAATTTGGTGACGCGATCCTCGAACAAATGTAATACCGTCGTTTGTTAAATTTATTAACGGGAGCGTAACGCTTCCGTTTTTTTTTGTTAGGCTGCTAAACGGTATCAAACTTTTTTCAAAAAAAGTTATCAAAAACCCCCGGCAGTTTTGGGGt                                     >  2:691220/1‑229 (MQ=255)
                                                                                                                                                                                     cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTTAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAACTTTAGAAAACAAATTTATCAAAACCCACCGGAAGTTTTGGGGTa                                    >  2:61125/1‑229 (MQ=255)
                                                                                                                                                                                     cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTa                                    >  1:337567/1‑229 (MQ=255)
                                                                                                                                                                                     cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGTTGGATGGCGCTAAACTGCTGAAATGTTCAGAGGTTGGTGACGCGATCATCGAAAACATGTAATGCCGTCGTTTGTTAAAATTATTAACGGGACCGGAACCCTCACGTTGGTTTTTGGTAGGCTGCTAAGCGTTATCAAAATTTTATCAACAAAATTTATCAAAACCCCCCGGTAGTTTTAGGGTa                                    >  2:583201/1‑229 (MQ=255)
                                                                                                                                                                                             aCGGGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGc                             <  1:670338/228‑1 (MQ=255)
                                                                                                                                                                                                  aaacCCGTACCTTTTGACTTCGAGTTTCTGATGGATGGCGCTAAACTGCTGAAATTTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAGCCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCa                       <  2:344937/229‑1 (MQ=255)
                                                                                                                                                                                                  aaaaCCGTAACCTATTGCTTCGAGTGTTTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACAGGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCa                       <  2:599776/229‑1 (MQ=255)
                                                                                                                                                                                                        gTAACCTTTGACTTCGGGGGTCTGATGGATGGCGATAAACTGAGGAAATGTTCAGAGGTTGGGGACGCGATCTTCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtggt                 <  2:491627/229‑1 (MQ=255)
                                                                                                                                                                                                        gTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGTCGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTACCGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtggt                 <  1:497512/229‑1 (MQ=255)
                                                                                                                                                                                                                      cGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGCTCATCGAACACATGTAATGCCGTAGTTTGTTAAATTTATTACCGGGACCGTAACGCTCCCGTGGTTTTGTGTTAGGATGCTAACGGTTATCAAAATTTTACCCAATAAAGTTTTTAAACACCCTCGGAGGTTTTGGGGTAGGTTGGCCCGTACGGTGGTAGTTctccttcta  >  2:745491/1‑228 (MQ=255)
                                                                                                                                                                                                                       gAGTGTCTGATGGATGGCGCTAAACTTCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCTTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGt  <  1:636195/229‑1 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
TCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGT  >  NC_000913/1196087‑1196530

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