Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,131 G→A 28.2% A337T (GCC→ACC)  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,1310GA28.2% 103.5 / 49.1 71A337T (GCC→ACC) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base A (7/12);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

ATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGA  >  NC_000913/1195894‑1196363
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        
aTCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGAATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCGGTTTGAAg                                                                                                                                                                                                                                          >  1:186938/1‑238 (MQ=255)
  gcgcGGGGGCCCGTGGCTGAAAGTTAAAACCCCGAACACGGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAAcc                                                                                                                                                                                                                                        <  2:558982/235‑1 (MQ=255)
   gACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGccac                                                                                                                                                                                                                                      <  1:916125/239‑1 (MQ=255)
     cGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGGGCAAACATCGGTGACG‑AAGGCGCCCTGTTTGAAAccaccc                                                                                                                                                                                                                                    >  1:967880/1‑239 (MQ=255)
      ggTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAAccaccca                                                                                                                                                                                                                                   <  1:394085/239‑1 (MQ=255)
      cgTGGCCCGTGGCTGAAAGTTAAAAACCCGAACTCTGGCAAAGAGATCGTCATTTAAGACGTGATTGCTGATGCCTTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGccaccca                                                                                                                                                                                                                                   <  1:59850/238‑1 (MQ=255)
          gCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAGGCACCCCACgga                                                                                                                                                                                                                               >  1:370263/1‑238 (MQ=255)
             cGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTAc                                                                                                                                                                                                                             <  1:957600/238‑1 (MQ=255)
               tGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGccaccca                                                                                                                                                                                                                                   <  2:199670/230‑1 (MQ=255)
               tGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGccaccca                                                                                                                                                                                                                                   >  1:199670/1‑230 (MQ=255)
               tGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTg                                                                                                                                                                                                                           >  1:708097/1‑238 (MQ=255)
                 gCTGAAATTTTAAAACCCGAACCCTGGGAAAGAGTTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTgcg                                                                                                                                                                                                                         <  2:967880/238‑1 (MQ=255)
                               aCCCCAACACTGGCAAAGGGATTGTAAATAAAGGCGTGATTGCTGATGCACTCCTGCAAATGATCCTGGTGCGTCCGGCGGAATATGAAGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGt                                                                                                                                                                                                          <  2:370263/239‑1 (MQ=255)
                                           gCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGCCTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGAATCGGTATCGCCCCTGGGGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAACCCACCCACTGTACTGCGCCGACATAATCCGGTCAGGGCAAAATa                                                                                                                                                                                              >  2:709401/1‑239 (MQ=255)
                                                  gaTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATcc                                                                                                                                                                                         <  1:135325/237‑1 (MQ=255)
                                                    tCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  2:708097/238‑1 (MQ=255)
                                                    tCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  1:644523/238‑1 (MQ=255)
                                                    tCGTCATTAAACGCGTTATTTTTGGGGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  2:529509/238‑1 (MQ=255)
                                                    tCGGCATTAAAGGCGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  1:549394/238‑1 (MQ=255)
                                                     cGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTGTCGCCTGTATGAACCTGAACGGTGACTACATTTCTGAAGCCCTGGCAGCGCAGGCTGGCGGTATCGGTATCGCCCGTGGGGC‑AACAACCGGGACGAAATCGGCCCTGCTTGAAGCCACCCCCGGTCCTGCCCCCCAAAAATCGGGGCGGGACGAAGTAAAACCCGGct                                                                                                                                                                                    >  2:417934/1‑239 (MQ=255)
                                                     cGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGGGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAAAATGCCGGTCAGGACAAAGAAAAACCCGGc                                                                                                                                                                                     >  2:135325/1‑238 (MQ=255)
                                                     cGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGGGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGAACTGCGCCGAAAAATGCCGGCCAGGACAAAGAAAAACCTgg                                                                                                                                                                                      >  2:71267/1‑237 (MQ=255)
                                                      gTCATTAAAGGCGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATTAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGct                                                                                                                                                                                    <  1:990410/238‑1 (MQ=255)
                                                                  gTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑CATGCGCCCTGTTTGAACCCACCCACGGTACTGCCCCGAAATATGCCGGTCAGGCCAAAGTAAATCCCTGctccttttcacaat                                                                                                                                                                        >  1:15166/1‑227 (MQ=255)
                                                                   tGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATtctctc                                                                                                                                                                       >  1:409262/1‑238 (MQ=255)
                                                                          tGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTgaga                                                                                                                                                               <  1:114620/239‑1 (MQ=255)
                                                                              gCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGGACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTTTCTCCCCTGAgatgat                                                                                                                                                           >  1:113218/1‑239 (MQ=255)
                                                                                  tccttgaaCAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGGACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGAtgctg                                                                                                                                                       <  1:325130/233‑1 (MQ=255)
                                                                                   ccTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCCAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCAACCACGGTACTGCGCCGAAATATGCCGGTCAGGACCAAGTAAATCCTGGCTCTCTTATTCTCTCCGCTGCGatgatgatg                                                                                                                                                       >  1:953330/1‑238 (MQ=255)
                                                                                       cAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  1:18928/203‑1 (MQ=255)
                                                                                        taCAGATCCTGCCGCGTCCGGCTGAATAAGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTAGCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                                                                                                 <  1:179667/238‑1 (MQ=255)
                                                                                         aCAGATCCTCCTGCGTCCGGCTGAATTGGATGTTAGCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcaca                                                                                                                                                <  1:417934/239‑1 (MQ=255)
                                                                                          cAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTTTCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGc                                                                                                                                                                    <  1:273467/218‑1 (MQ=255)
                                                                                              tCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGt                                                                                                                                           >  1:839888/1‑239 (MQ=255)
                                                                                                 tgctgcGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGGACTGCGCCGAAATATGCCGGTCCGGGCAAAGTAAATCCTGGCCCTattatt                                                                                                                                                                            >  2:18928/1‑203 (MQ=255)
                                                                                                    tgcGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGAAGCTGCGCCACATGGGTTGGAcc                                                                                                                                     >  2:177074/1‑239 (MQ=255)
                                                                                                        tCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                                                                                                  >  1:985930/1‑238 (MQ=255)
                                                                                                        tCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa                                                                                                                                  >  1:318699/1‑238 (MQ=255)
                                                                                                                          ttccGCCTTTATGAACCTGAACGGTGCCTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                                                               <  2:839888/236‑1 (MQ=255)
                                                                                                                          taacGCCTGGATGAACCTGAACGGTGACTACATTTCTGACCGCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                                                               <  1:71267/236‑1 (MQ=255)
                                                                                                                          tATCGCCTGTTTGAACCTGAACGGTGACTTCATTTCTGGCGCCCTGGGAGGGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAAGCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                                                               <  2:1418/239‑1 (MQ=255)
                                                                                                                          tATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      >  1:645285/1‑168 (MQ=255)
                                                                                                                          tATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  2:645285/168‑1 (MQ=255)
                                                                                                                          tATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGttc                                                                                                               >  1:625537/1‑238 (MQ=255)
                                                                                                                          tATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                                                               >  1:555929/1‑239 (MQ=255)
                                                                                                                          tATCGCCTGTATGAACCTGAACGGTGACTACATTGCCGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGACCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                                                               <  1:526960/239‑1 (MQ=255)
                                                                                                                                  gTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCACACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGGACTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTTTTTCCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGCCCGAAGCGGCGGACTTAATTTTTAAAgtaatg                                                                                                       >  2:645686/1‑234 (MQ=255)
                                                                                                                                   tATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGCAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGatgatg                                                                                                                                                          >  2:739940/1‑187 (MQ=255)
                                                                                                                                   tATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGatgatg                                                                                                                                                          <  1:739940/187‑1 (MQ=255)
                                                                                                                                   tATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGCTGATGCTGGGCCACATGGGTTGGACGGAAGCGGCGTACTTAATTGTTAAAGGGATgg                                                                                                      >  1:985429/1‑239 (MQ=255)
                                                                                                                                          cTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGCCG‑AATGCGCACTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGCAAATCCAGGCTCCATTATTACTTCCCCTGAGTGAAAGCTTCGTCCAATGGGTGGGACGGACAGGGCTGACTTTATTTTTAAAGATATGGAAAgccg                                                                                               >  2:144495/1‑237 (MQ=255)
                                                                                                                                             aaCGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCGGCGCCACATGGGTTGGACCGAGACGGCTGAGTTAATTGTTAACGGTATGGAAAGGGCaaa                                                                                            >  2:697592/1‑238 (MQ=255)
                                                                                                                                               cGGGGACTACATTTCTGACGCCCTGGCAACGCAGGTTTGCGGTATCGGTATCGGCCCTGGTGCAAACATCGGGGACA‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTTACTTAATTGTTAAAGGTATGGAAGGCGCAATc                                                                                           <  1:645686/238‑1 (MQ=255)
                                                                                                                                                cttgccTACTTTTTTGCCGCCCTGGCAGCGCAGGTTGGGGGTTTCGGTATCGGCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATc                                                                                           <  2:186938/232‑1 (MQ=255)
                                                                                                                                                  tGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACCGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTTTTAAAGGTATGGAAGGCGCAATAAAcc                                                                                       >  1:354011/1‑238 (MQ=255)
                                                                                                                                                     cTACATTTCGGACGCCCTGGCCGCGGAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCGGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGa                                                                                    <  2:722063/239‑1 (MQ=255)
                                                                                                                                                       aCATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCGCAATCAACGCGaaa                                                                                  >  1:51184/1‑239 (MQ=255)
                                                                                                                                                           ttCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGTTAACGCCCCTGGTGCAAACGTCGGGTACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCg                                                                              <  1:709401/239‑1 (MQ=255)
                                                                                                                                                                       ggCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCGGAGCAGATCCTGCGCCAAATGGGTTGGACCGAAGCGGCTGCATTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGAc                                                                   >  1:1055179/1‑238 (MQ=39)
                                                                                                                                                                            cgcAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGATGCGGCTGACTTAATTGTTAAAGGTAAGGAAGGCGCAATCAACGGGAAAACCGAAACCttggctttcgag                                                             >  1:504874/1‑228 (MQ=255)
                                                                                                                                                                             gcAGGTTGGCGGTTTCGGTATCGCCCCTGGTGCAAACATCGGGGACC‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg                                                             <  2:504874/238‑1 (MQ=255)
                                                                                                                                                                                        gtatcggtatcgCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                  <  1:177074/238‑1 (MQ=255)
                                                                                                                                                                                           tcggtatcgCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTa                                                                                                                      <  2:747851/167‑1 (MQ=255)
                                                                                                                                                                                           tcggtatcgCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTa                                                                                                                      >  1:747851/1‑167 (MQ=255)
                                                                                                                                                                                           tcggtatcgCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGACAGCGCAAGCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGAAGGAGGgc                                              >  2:603886/1‑239 (MQ=255)
                                                                                                                                                                                               tctcgCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGAAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGCATGGAAGGCGCAATCAACGCGAAAACCGTACCCTATGACTTCGAGCGTCTGATGAATGGCGCTa                                          >  1:549464/3‑239 (MQ=255)
                                                                                                                                                                                                agcgCCCCTGGTGCAAACCTTGGTGACG‑AATGGGCCCTGTTTGAAGCCCCCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTa                                          <  2:623414/236‑1 (MQ=255)
                                                                                                                                                                                                 tcgCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTaaa                                        <  2:625537/239‑1 (MQ=255)
                                                                                                                                                                                                         ggTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGCCAAAGTAAATCCTgg                                                                                                                                                                                      >  2:313854/1‑89 (MQ=255)
                                                                                                                                                                                                         ggTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  1:313854/89‑1 (MQ=255)
                                                                                                                                                                                                                                  gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGa         >  1:342374/1‑238 (MQ=255)
                                                                                                                                                                                                                                        tttGAAGCCCCCCCCGGTTCTTGGCCGAAATAGGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGa  <  2:985930/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        
ATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACG‑AATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGA  >  NC_000913/1195894‑1196363

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