Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 71,085 T→C S233S (AGT→AGC araC → DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,0850TC88.7% 183.2 / 11.3 72S233S (AGT→AGCaraCDNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base C (32/31);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATC  >  NC_000913/70839‑71322
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
ggCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCTTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAgg                                                                                                                                                                                                                                           >  1:563583/1‑251 (MQ=255)
      aaaTCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAg                                                                                                                                                                                                                                      <  1:202939/250‑1 (MQ=255)
           tGCTTGAGCAATTGTTACTACGGCGCATGGCAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCCGTACATCAGCGATCACCTGGCAGACAGCACTTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCCGCAGTTAGGGATTAGCGACTTAATCTGGCCCGCGGACCAACGTATCAGCCAGGCGCAgctgatt                                                                                                                                                                                                                                >  2:626357/1‑248 (MQ=255)
               ttaGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTgg                                                                                                                                                                                                                            >  2:198541/3‑250 (MQ=255)
                taGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGa                                                                                                                                                                                                                            <  2:383279/249‑1 (MQ=255)
                  gCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                                                                                                                                                                                                          >  1:255519/1‑250 (MQ=255)
                     aTTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAccc                                                                                                                                                                                                                      >  1:313733/1‑249 (MQ=255)
                            cTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTCTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCAGCTTTTGAGCACCACCCGGAt                                                                                                                                                                                                               >  1:24723/1‑251 (MQ=255)
                                    caaggAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGGCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGTATTAGCGTCTTAAGCTGGCGCGAGTACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCGCCCGGATGCCTATCg                                                                                                                                                                                                       <  1:88466/248‑1 (MQ=255)
                                            cGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtc                                                                                                                                                                                               <  1:299749/251‑1 (MQ=255)
                                                     aGTCGCTCCATCCACCGAGGGATATGCGGGTACGCCTGGCTTTTCAGTACATCAGCGATTACCTTGCAGTCAGCAATTTTGATATCGCCAGCGTCGCACAGCTTGTTTGCTTGTCGCCGTCGCGTCTGTCTCATCTTTTCCGCCAGCAGTTAGGGATTAGCGGCTTAAGCTGGCGCGAGGACCAACGCTTTAGTCAGGCGAAGCGGCTTTTGAGCACTACCCGGATGCCTATCGCCTCCGTCGGTCGCAAt                                                                                                                                                                                      <  2:202501/251‑1 (MQ=255)
                                                      gTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATg                                                                                                                                                                                     >  1:572708/1‑251 (MQ=255)
                                                       tCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGt                                                                                                                                                                                    >  2:185233/1‑251 (MQ=255)
                                                        cGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCACTGtt                                                                                                                                                                                   >  2:455265/1‑251 (MQ=255)
                                                        cGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGGTCGGGATTAGCGTCTTAAGCTGGCGCGAGGCCCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGCTGCCTATCGCCCCCGTCGGTCGCTATGtt                                                                                                                                                                                   >  2:8974/1‑251 (MQ=255)
                                                              atccaCCGATGGCTAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCGACTGGCAGACAGCAATTTTGATATCGCCCACGTCGCACAGCATGTTTGCTAGTCGCCGTAGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGCTTAGCGTCTTAAGCTGGCGAGCGGACCAACGGAACAGCCAGGCGAAGCTGCTTTTGACCACCACCCGGAAGCCTATCGcaccaggaagtcccaatgtgtatttt                                                                                                                                                                             >  2:465011/1‑226 (MQ=255)
                                                                     gATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                                                                                                                                                                      >  1:633644/1‑251 (MQ=255)
                                                                     gATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                                                                                                                                                                      >  1:158262/1‑251 (MQ=255)
                                                                      gtGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCTCCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCTAGCAGTTAGGGATTTGCGTCTTAAGCTGGCGCGAGGTCCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCGCCCCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCa                                                                                                                                                                     <  1:408862/250‑1 (MQ=255)
                                                                          aTAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAActc                                                                                                                                                                 <  1:545385/251‑1 (MQ=255)
                                                                               ggggTACGCGAGGCTTGTCAGGACATCAGCGATCACCTGGCAGACAGCAATTGTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCATTACCCGGATTCCTACCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCCAttt                                                                                                                                                            <  2:162183/250‑1 (MQ=255)
                                                                                gagTACGCGAGGCTTGTCAGTACCTCAGCGATCACCTGGCAGCCAGCAATTTTGATATCGCCAGCGTCGCACAGCCTGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCACCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGCCCAACGCATTAGTCAGGCGACGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCCCCGTCGGTCGCAATGTTGGTTTTGGCGCTCCACCCTAtttt                                                                                                                                                           >  2:286131/3‑250 (MQ=255)
                                                                                 ggTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtct                                                                                                                                                          <  1:317415/251‑1 (MQ=255)
                                                                                 ggTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtct                                                                                                                                                          >  1:160645/1‑251 (MQ=255)
                                                                                     cgcgATGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTcgcg                                                                                                                                                      <  1:234421/251‑1 (MQ=255)
                                                                                       cgAGGCTTGTCAGTACATCAGCGATCACCTGGCCGACAGCAATTTTGATATCGCCAGAGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGGCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGATGCTTTTGAGCACCACCCGGCAGCCTATCGCCACCGACGGTCGCAATGATGGTTTTGCCGATCACCTCTATTTCTcgcggg                                                                                                                                                    >  2:161095/1‑249 (MQ=255)
                                                                                            cTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGCGCACCACCCGGCTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTTACGATCAACTCCTTTTCTCGCGGGTAttt                                                                                                                                               >  2:112817/1‑251 (MQ=255)
                                                                                            cTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                                                                                                               >  2:137152/1‑251 (MQ=255)
                                                                                               gTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                                                                                                             <  2:572708/250‑1 (MQ=255)
                                                                                                 cAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                                                                                                          <  1:455265/251‑1 (MQ=255)
                                                                                                  aGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                                                                                                          >  1:531450/1‑250 (MQ=255)
                                                                                                          aGCGATCACCTGGCAGACAGCAATTTTGATTTCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTGTGCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                                                                                                 <  2:527735/251‑1 (MQ=255)
                                                                                                          aGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                                                                                                 >  1:223365/1‑251 (MQ=255)
                                                                                                             gATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcggtcg                                                                                                                                                                                          <  1:151741/191‑1 (MQ=255)
                                                                                                                  ccTGGCAGACAGCAATTTTGATAACGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGGCGCGTCTGTCACAACTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCTCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCCCCGGGGCCAGcc                                                                                                                         <  2:423538/251‑1 (MQ=255)
                                                                                                                     ggCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                                                                                                       >  1:4389/1‑250 (MQ=255)
                                                                                                                        aGACATCAATTTTGGTAGCGACAGTGTCGCACAGCATGGTTGCTTGTCGCGGTTGCGTCTGTCACCTCTTGTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCTACGCGTTAGGGAGGCGAAGCTGCTTTTGAGCTGTATCCGGATTCCTATCGCCATCGGCGGTCGCAATGTTGGTTTTGACGATCAACTCTATGTCTCGCGATTATTTAATGAATTCACCGGGGGCAGTCcgagcg                                                                                                                   <  2:356394/251‑1 (MQ=255)
                                                                                                                           cAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                                                                                                                >  1:507404/1‑251 (MQ=255)
                                                                                                                              caagtatgtgatCGCTAGCGTCGCACAGCATGTTTGTTTGACGCCGTCGCGTCTGTCACAGCTTTTCCGACAGCAGTTAGGGAGTAGCGTCTTAAGCTGGTGCTAGGCCCAACGTATCAGCCAGGCGTAGCTGCTTTTGAGCACCACTCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCATCTCTATTTCGCGCGGGTATTTAAAAAATGCACCGGTGCCAGCCCGAGCGAGttcc                                                                                                             <  2:295647/241‑3 (MQ=255)
                                                                                                                              cAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCCCATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttcc                                                                                                             <  2:507404/251‑3 (MQ=255)
                                                                                                                                      atatCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCggg                                                                                                     >  1:123193/1‑250 (MQ=255)
                                                                                                                                        atCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACAATCACCTCTATTTTTCGTGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGAGTGTTCCGGGCCGGTtg                                                                                                   >  2:614003/1‑251 (MQ=255)
                                                                                                                                          cGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                                                                                                 <  1:533456/251‑1 (MQ=255)
                                                                                                                                              aGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtc                                                                                                                                                           <  1:565473/189‑1 (MQ=255)
                                                                                                                                              aGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtc                                                                                                                                                           >  2:565473/1‑189 (MQ=255)
                                                                                                                                               gCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa                                                                                            <  2:531450/251‑1 (MQ=255)
                                                                                                                                               gCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaaga                                                                                             >  1:310093/1‑250 (MQ=255)
                                                                                                                                               gCGTCGCACAGCATGTTTGCTTGTCGCCGGCGCGTCTGTCACATCTTTTCCGCCAGCAGGTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa                                                                                            <  2:255519/251‑1 (MQ=255)
                                                                                                                                                cGTCGCTCAGCATGTTTGCTTTGCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGAGTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa                                                                                           <  2:530398/251‑1 (MQ=255)
                                                                                                                                                     cacaGCATGTTTGCTTGTAACCGTCGCGTCTGTCACATCTTTTCCGCCGGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCCGGCGAAGCTGCTTTTGCGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCCATGTTGGTTTTGAAGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGCATTCCGTGCCGGTTGTGCAGAAAAAgcg                                                                                      >  2:701454/1‑249 (MQ=255)
                                                                                                                                                        aGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAATATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat                                                                                   <  1:198541/251‑1 (MQ=255)
                                                                                                                                                         gCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAACTGCACCGGGGCCAGCCCGAGCGGGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg                                                                                  >  2:498181/1‑251 (MQ=255)
                                                                                                                                                          cATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTGTTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAggga                                                                                     >  2:570853/1‑244 (MQ=255)
                                                                                                                                                          cATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTGTTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa                                                                                     <  1:570853/247‑1 (MQ=255)
                                                                                                                                                              tttGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACAAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGt                                                                              >  2:609627/1‑250 (MQ=255)
                                                                                                                                                                       cGCCGTCTCGGCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAg                                                                    <  2:310093/251‑1 (MQ=255)
                                                                                                                                                                          cgtcgcgtcgGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTg                                                                 <  2:123193/251‑1 (MQ=255)
                                                                                                                                                                             cgcgTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCa                                                              <  1:667811/251‑1 (MQ=255)
                                                                                                                                                                               cgTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCAGTCAAGTTGTCAt                                                             >  1:4773/1‑250 (MQ=255)
                                                                                                                                                                               cgTCTGGCACATCTTTGCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAACACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATa                                                            <  2:223365/251‑1 (MQ=255)
                                                                                                                                                                                       acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTa                                                     >  1:654740/1‑250 (MQ=255)
                                                                                                                                                                                                ccGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAg                                            <  1:614003/250‑1 (MQ=255)
                                                                                                                                                                                                  cccAGCAGTGAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTGTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACa                                         <  2:672794/250‑1 (MQ=255)
                                                                                                                                                                                                     agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAAtt                                      <  1:585920/251‑1 (MQ=255)
                                                                                                                                                                                                       cagTTAGGGATTTGCGTCTTAAGCTGGCGCAAGGTCCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTGTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAAtcga                                    <  2:563583/251‑4 (MQ=255)
                                                                                                                                                                                                       cagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGa                                    <  1:98586/251‑1 (MQ=255)
                                                                                                                                                                                                          ttAGGGATTACCGTATTAAGCTGGCACGAGGACAAACGTATCAGCAAGGCGAAGCTGTTTTTGAGCACCACACGGAGGACAATCGCCACTGTCGGTCGCAATGTTGGTTTTGACGATCAACGCTATTGATCGCGGGGATTAAACAAATGCAACGGGGTAAGACCGGGCGTGATCCGGGATGGATGTGAAGGAACAGTGCATGATGCAGGCGTCAAGTGGTGATAATTGGGAATGAATCAGACAATTggtgg                                 >  2:108024/1‑247 (MQ=255)
                                                                                                                                                                                                            aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTACTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGc                                <  1:175268/250‑1 (MQ=255)
                                                                                                                                                                                                             gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCtt                              <  1:580660/251‑1 (MQ=255)
                                                                                                                                                                                                                    gCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt                                                                                                  <  1:379362/176‑1 (MQ=255)
                                                                                                                                                                                                                    gCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt                                                                                                  >  2:379362/1‑176 (MQ=255)
                                                                                                                                                                                                                                         aCCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGCGTAGCATAGGGTTTGCAGAATc  >  1:530227/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
GGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATC  >  NC_000913/70839‑71322

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