Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,305 C→T 40.1% 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,3050CT40.1% 25.8 / 73.7 74R395C (CGT→TGT) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base T (16/13);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-03

GTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGA  >  NC_000913/1196078‑1196542
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
gtatcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                                                                                                                                                                                                      <  1:177074/238‑1 (MQ=255)
   tcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGACAGCGCAAGCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGAAGGAGGgc                                                                                                                                                                                                                                  >  2:603886/1‑239 (MQ=255)
       tctcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGAAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGCATGGAAGGCGCAATCAACGCGAAAACCGTACCCTATGACTTCGAGCGTCTGATGAATGGCGCTa                                                                                                                                                                                                                              >  1:549464/3‑239 (MQ=255)
        agcgCCCCTGGTGCAAACCTTGGTGACGAATGGGCCCTGTTTGAAGCCCCCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTa                                                                                                                                                                                                                              <  2:623414/236‑1 (MQ=255)
         tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTaaa                                                                                                                                                                                                                            <  2:625537/239‑1 (MQ=255)
                                         gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGa                                                                                                                                                                                             >  1:342374/1‑238 (MQ=255)
                                               tttGAAGCCCCCCCCGGTTCTTGGCCGAAATAGGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGa                                                                                                                                                                                      <  2:985930/239‑1 (MQ=255)
                                                         cccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCCAAGTTTGGGGCAGCGATCATCGaaacaagttaa                                                                                                                                                                            >  2:988611/1‑231 (MQ=255)
                                                         cccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTCCTGAATTGTTCTGAGTTTGGGGACGCGACAATCCAAAAAATGaaa                                                                                                                                                                           >  2:308554/1‑239 (MQ=255)
                                                                 aCTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGCAAACATGTAATGCCGTAg                                                                                                                                                                   >  1:944262/1‑239 (MQ=255)
                                                                    gcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGttt                                                                                                                                                                >  1:227095/1‑239 (MQ=255)
                                                                     cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                                                                                                                    <  1:795029/219‑1 (MQ=255)
                                                                     cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                                                                                                                    >  2:795029/1‑219 (MQ=255)
                                                                     cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACGGCTGAAATGTTCAGAGGTTGGTGACGCGATCATAGAAAACATGTAATGCCGTAGTTTg                                                                                                                                                               >  1:65805/1‑239 (MQ=255)
                                                                       ccGAAATTTGGCGGTAAGGAAAAAGTAAATCCTGGTTCTATTTTTCCTTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt                                                                                                                                                             <  2:113218/239‑1 (MQ=21)
                                                                       ccGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt                                                                                                                                                             <  2:409262/239‑1 (MQ=255)
                                                                         aaaaTATGCCGGTTAGGACAAAGTAAATCCTGGCTCTTTGATTCTCTCCGCTGAGATGATGCCGCGCCACATGGGTTGGACCGAAGCCGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTaa                                                                                                                                                           <  2:65805/238‑1 (MQ=255)
                                                                          aaaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGTAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTCCTGAAATGTTCAGAGTTTGGTGACGCGTTCATCGTAAACATGTAATGCCGTAGTTTGTTaca                                                                                                                                                          >  1:247635/1‑237 (MQ=255)
                                                                          aaaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTaaa                                                                                                                                                          >  1:585362/1‑239 (MQ=255)
                                                                                      ctggACAAAGTAAATCCGTGCTCTTTTTTTCCTCCCGCTGAGAGTATGCTGCGCCAAATGGGTTGGACCGAAGGGGCTTACTTAATTGGTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACgg                                                                                                                                               <  2:985429/236‑1 (MQ=255)
                                                                                      cAGGACCAAGTAAATCCTGGCTCTATTATTCTCTCCGCTTGGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTTACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACgg                                                                                                                                               <  2:227095/238‑1 (MQ=255)
                                                                                              aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACGGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGCGCGTaa                                                                                                                                       >  1:972585/1‑238 (MQ=255)
                                                                                                     ccTGGCTTTTTTATTCTCTCCGCTGGGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCg                                                                                                                               <  2:51184/239‑1 (MQ=255)
                                                                                                                 attCTCTCCGCTGGGATGAAGCAGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAATCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGtt                                                                                                                   <  1:308554/239‑1 (MQ=255)
                                                                                                                    ctctcCGCTTGGATGAGGCTGGGGCAACTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg                                                                                                                <  2:318699/239‑1 (MQ=255)
                                                                                                                    ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcat                                                                                                                                                                                         <  2:153233/167‑1 (MQ=255)
                                                                                                                    ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcat                                                                                                                                                                                         >  1:153233/1‑167 (MQ=255)
                                                                                                                    ctcgcCGCTTAGAGTATGCTGGGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg                                                                                                                <  2:354011/239‑1 (MQ=255)
                                                                                                                       tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                                                                                                                    <  1:761550/169‑1 (MQ=255)
                                                                                                                       tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                                                                                                                    >  2:761550/1‑169 (MQ=255)
                                                                                                                       tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAg                                                                                                                 >  2:851648/1‑235 (MQ=255)
                                                                                                                       tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCGGATGGATGGAGCTAAACTGCTGAAATGTTCAGCGGTTGGTGACGCGATCATCGAAAACATGTAAGGCCGTAGTTTGTTGAATTTATTAACGCGAGCGTAACGCCTCCGTTGTTTTTTGTTATgctg                                                                                                             >  1:385505/1‑239 (MQ=21)
                                                                                                                         cggtgcgatgatgatGCGGCCCATGGGTTGGACCGGAGCGGCTGACTTAATTGGTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACGTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATGATCGAAAACATGTAATGCCGTAGTTTGTTAAAGTTATTAACGGGAGCGTAACGCTCTCGTTGTTTTTTGTTAGgctgct                                                                                                           <  1:862052/233‑1 (MQ=255)
                                                                                                                         cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGa                                                                                                                                                                                             <  1:134365/158‑1 (MQ=255)
                                                                                                                         cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGa                                                                                                                                                                                             >  2:134365/1‑158 (MQ=255)
                                                                                                                                atgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTGTTTTTTTTTTAGGCTGCTAAAGGtt                                                                                                    >  2:866427/1‑239 (MQ=255)
                                                                                                                                          gcCACATGGGTTGGACCGAAGCGGCAGACTTAATTGTTAAAGGGATGGACGGCCCAATCACCGCGAAAACCGGAACCTATGACTTCGAGCGTCTGATGGATGGCGCCAAAATGCTGAATTCTTCAGAGTTTGGTGACGCGACCCTCGAAATATGTTAATGCAGTGGTTTGTTAAATTGTTTACAGGTACCCAAAGGCTCCCGTTGTTTTTTGTTAGTCCGCTAAAGGTTATCAAAAttt                                                                                          >  2:275353/1‑239 (MQ=255)
                                                                                                                                           ccACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAtttt                                                                                         >  2:541374/1‑239 (MQ=255)
                                                                                                                                                   gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTTTCaaaaa                                                                                 >  1:653245/1‑239 (MQ=255)
                                                                                                                                                                               tAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                                                                                                                    <  1:75714/113‑1 (MQ=255)
                                                                                                                                                                               tAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                                                                                                                    >  2:75714/1‑113 (MQ=255)
                                                                                                                                                                                taaGGTTTGGAAGGCCGAATCAACGGGAAAGCCTGAACCTATGAGTGCGGGCGTTGGAGGGAGGGCGCTAAACTGCTGAAATGTTCAGATTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCTCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgg                                                    <  2:972585/238‑1 (MQ=255)
                                                                                                                                                                                aaaGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGCTGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCCCCCGTTGTTTTTTGTTAGGCTGGTAACGGTTTCCAAAATTTTTTCAAAAAAAGTTATACAAACCCCTCGGCGGTTTTgg                                                    >  2:84898/1‑239 (MQ=255)
                                                                                                                                                                                        ggAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCGGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg                                               >  2:302335/1‑236 (MQ=255)
                                                                                                                                                                                        ggAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGCCTGCTAACGGTAATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTCGTTTTTGGGTAg                                               >  1:729912/1‑236 (MQ=255)
                                                                                                                                                                                        ggAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGGAATGTTCAGAGTTTGGTGACGCGATCATCGACAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg                                               >  2:429389/1‑236 (MQ=255)
                                                                                                                                                                                        ggAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGGAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTCTTAACGGGCGCGTAAAGCTCCCGTTGTTTTTTGTTCGGCTGCTAACCGTTATCAAAATTTTTTCAAAAAAATTTTTCAAAACCCCTCGGGAGTTTTGGGGTAGGcg                                            >  2:436456/1‑238 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGACAACCGTAACCTGTGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATACCGTAGTTTGTTAAGTTTATTAACGGTAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAACCCTCGGTGGTTTTTGGGTAGGGTGGCCgg                                      >  2:647666/1‑239 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTTTTAGGCTGCTACCGGTTATCAAAATTTTATCAAAAAAAGTTATAAAAACCCCCCGGTAGTTTTTGGGTAGGCTGGCCg                                       >  2:484028/1‑238 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTTTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTTTCAAAACCCCCGGGGGGTTTTGGGGAAgggggggcgg                                      >  2:1075768/1‑231 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCgg                                                             >  2:902154/1‑216 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCgg                                                             <  1:902154/216‑1 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg                                      >  2:905536/1‑239 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg                                      >  2:993017/1‑239 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTACATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTTGGGTAggagggacg                                       >  2:206879/1‑231 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg                                      >  2:816977/1‑239 (MQ=255)
                                                                                                                                                                                              cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg                                      >  2:716099/1‑239 (MQ=255)
                                                                                                                                                                                               gcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTTTCAAAAAAAGTTATCAAAACCCCTCGGAGGTTTTGGGTGAGGCTGGCgggg                                     >  2:468638/1‑238 (MQ=255)
                                                                                                                                                                                                 aaTCAACGCGGAAACCTTAACCTTTGGCTTCGAGCGTCTGATGGATGGCGCTAAACTGCGTAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCa                                   <  2:247635/239‑1 (MQ=255)
                                                                                                                                                                                                       cgggAAAACCGTTACCTTTGACTTCGGGCGTTTGATGGATGGCCATAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtggt                             <  2:342374/236‑1 (MQ=255)
                                                                                                                                                                                                       cccgAAAACCGTAACCTTTGGCTTCGAGCGGCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtggt                             <  2:555929/237‑1 (MQ=255)
                                                                                                                                                                                                               ccGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTctac                     <  1:697592/239‑1 (MQ=255)
                                                                                                                                                                                                                  tAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGTTGGCCGGCCAGGGGGGAGTTctaatac                  >  2:1007427/1‑239 (MQ=255)
                                                                                                                                                                                                                    aCCTATGTCTGCGAGTGTCTGATGGAGTGCGCTAAACTGGTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTctactacta                <  1:647666/239‑1 (MQ=255)
                                                                                                                                                                                                                      cTATTACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGt              <  1:603886/239‑1 (MQ=255)
                                                                                                                                                                                                                      cTATGGCGTCGAGCGGCTGATGGATGGCGCTAAACTGCTGAAAGGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGt              <  1:436456/239‑1 (MQ=255)
                                                                                                                                                                                                                      cTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGt              <  1:988611/239‑1 (MQ=255)
                                                                                                                                                                                                                      cTATGACTTCGAGCGTCTGATTGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGt              <  1:206879/239‑1 (MQ=255)
                                                                                                                                                                                                                         tGACTTCGAGTGTCTGATTGATGGCGCTAAACTTCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGATCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGtctc           <  1:1075768/239‑1 (MQ=255)
                                                                                                                                                                                                                         tGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCCAAAAAAGTGATCAAAACCCCTCGGAAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGGTctactcctaattt            >  1:201244/1‑234 (MQ=255)
                                                                                                                                                                                                                          tgcTTCGGGGGTTTGAGGGTGGGCGATAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTcc          <  2:1055179/237‑1 (MQ=255)
                                                                                                                                                                                                                          gACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTcc          <  1:816977/239‑1 (MQ=255)
                                                                                                                                                                                                                             ttCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCcac       <  1:798461/239‑1 (MQ=255)
                                                                                                                                                                                                                                gAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGAAGTTCTACTACTAGTCTCCCACAtg    >  1:1031274/1‑238 (MQ=255)
                                                                                                                                                                                                                                   cGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGa  <  1:851648/238‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
GTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAGA  >  NC_000913/1196078‑1196542

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