Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 70,371 T→C 100% intergenic (‑323/‑16) araB ← / → araC ribulokinase/DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091370,3710TC98.3% 29144.0 / 234.0 8171intergenic (‑323/‑16)araB/araCribulokinase/DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base T (70/76);  new base C (3856/4168);  total (3927/4244)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 8324 total aligned reads displayed.

AGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGG  >  NC_000913/70273‑70468
                                                                                                  |                                                                                                  
  cGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAAtggcgg                                                                                                  <  2:2653726/99‑1 (MQ=255)
    gggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGa                                                                                                >  1:119274/1‑99 (MQ=255)
    gggTTACAGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGTATGGCGGGa                                                                                                <  2:13210/99‑1 (MQ=255)
     ggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAg                                                                                               <  1:3021125/99‑1 (MQ=255)
     ggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAg                                                                                               >  1:813195/1‑99 (MQ=255)
      gTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGt                                                                                              >  2:1682155/1‑99 (MQ=255)
       ttACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTa                                                                                             >  2:856861/1‑99 (MQ=255)
       ttACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTa                                                                                             <  1:1486264/99‑1 (MQ=255)
       ttACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTa                                                                                             >  2:2515530/1‑99 (MQ=255)
       ttACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTa                                                                                             >  2:1824814/1‑99 (MQ=255)
              tttggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAg                                                                                     <  2:2439538/96‑1 (MQ=255)
              tttggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAg                                                                                     <  1:2591425/96‑1 (MQ=255)
                tggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTa                                                                                   <  2:403998/98‑1 (MQ=255)
                tggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTa                                                                                   <  1:1623089/98‑1 (MQ=255)
                tggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTa                                                                                   <  2:2305566/98‑1 (MQ=255)
                tggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTa                                                                                   <  2:3678237/98‑1 (MQ=255)
                 ggTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTAt                                                                                  <  2:1963091/99‑1 (MQ=255)
                   ttAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATgg                                                                                >  1:2942712/1‑99 (MQ=255)
                   ttAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATgg                                                                                >  2:1248706/1‑99 (MQ=255)
                   ttAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATgg                                                                                >  2:3020578/1‑99 (MQ=255)
                         agaagaGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAg                                                                          >  2:482779/1‑99 (MQ=255)
                          gaagaGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgc                                                                         <  1:2065257/99‑1 (MQ=255)
                          gaagaGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgc                                                                         >  1:3098612/1‑99 (MQ=255)
                          gaagaGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgc                                                                         >  2:3149873/1‑99 (MQ=255)
                           aagaGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgcg                                                                        >  1:1779706/1‑99 (MQ=255)
                            agagCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgcgc                                                                       >  1:3227733/1‑99 (MQ=255)
                             gagCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCa                                                                      >  1:1345356/1‑99 (MQ=255)
                              agCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaa                                                                     >  2:925461/1‑99 (MQ=255)
                               gCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaa                                                                    >  2:3285838/1‑99 (MQ=255)
                                ccAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaaa                                                                   <  1:280253/99‑1 (MQ=255)
                                ccAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaaa                                                                   >  1:431779/1‑99 (MQ=255)
                                ccAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaaa                                                                   <  2:1479415/99‑1 (MQ=255)
                                ccAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaaa                                                                   <  1:1553883/99‑1 (MQ=255)
                                ccAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaaa                                                                   <  1:2734441/99‑1 (MQ=255)
                                 cAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAAt                                                                  >  1:2668789/1‑99 (MQ=255)
                                  aGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATg                                                                 <  2:3351816/99‑1 (MQ=255)
                                    tAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGAAAAGTATGGCTGAAGCGCAAAATGAt                                                               >  1:499334/1‑99 (MQ=255)
                                    tAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGAt                                                               >  1:3623560/1‑99 (MQ=255)
                                    tAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGAt                                                               >  1:2410191/1‑99 (MQ=255)
                                       aaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATccc                                                            <  1:3415846/99‑1 (MQ=255)
                                       aaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATccc                                                            <  1:2874387/99‑1 (MQ=255)
                                       aaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATccc                                                            <  2:1167331/99‑1 (MQ=255)
                                       aaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATccc                                                            <  1:2337543/99‑1 (MQ=255)
                                          aCGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCctg                                                         <  2:2372007/99‑1 (MQ=255)
                                           cGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCctgc                                                        <  2:2910739/99‑1 (MQ=255)
                                                tcaCGGCAATGTCTGATGCCATATGGACAATTGGTTTCTTCCCTGATTGGGGGGAGTACGAACAGGATGGCTGAAGTGCATAATGATCCCCTGCTGCCg                                                   >  2:740187/3‑99 (MQ=255)
                                                tGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCg                                                   >  1:676450/1‑99 (MQ=255)
                                                tGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCg                                                   >  1:556456/1‑99 (MQ=255)
                                                tGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCg                                                   >  1:354987/1‑99 (MQ=255)
                                                     gCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATa                                              <  1:377760/99‑1 (MQ=255)
                                                     gCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATa                                              >  2:2721188/1‑99 (MQ=255)
                                                      cAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATAc                                             <  1:1984847/99‑1 (MQ=255)
                                                        aTGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTc                                           >  1:1001523/1‑99 (MQ=255)
                                                        aTGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTc                                           >  2:617580/1‑99 (MQ=255)
                                                         tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCg                                          >  1:3162019/1‑99 (MQ=255)
                                                         tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCg                                          >  1:348344/1‑99 (MQ=255)
                                                         tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCg                                          <  2:3409014/99‑1 (MQ=255)
                                                         tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCg                                          >  2:548182/1‑99 (MQ=255)
                                                           tCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGtt                                        <  2:800723/99‑1 (MQ=255)
                                                           tCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGtt                                        <  1:1120566/99‑1 (MQ=255)
                                                           tCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGtt                                        <  2:3473178/99‑1 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                     <  2:2795370/99‑1 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                     <  1:170414/99‑1 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                     <  2:296980/99‑1 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                     <  2:3217275/99‑1 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                     <  2:1893088/99‑1 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                     <  1:213243/99‑1 (MQ=255)
                                                              gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaa                                     <  2:1610640/99‑1 (MQ=255)
                                                                tGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaat                                    <  2:2098827/99‑2 (MQ=255)
                                                                tGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaatg                                   >  1:1916942/1‑97 (MQ=255)
                                                                 gCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTaatgc                                  >  1:3292096/1‑96 (MQ=255)
                                                                   aaTATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGccc                                >  2:2246211/1‑99 (MQ=255)
                                                                   aaTATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGccc                                >  1:1267371/1‑99 (MQ=255)
                                                                   aaTATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGccc                                >  1:750386/1‑99 (MQ=255)
                                                                    atatGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCa                               >  2:674376/1‑99 (MQ=255)
                                                                     tatGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCAt                              <  2:2026092/99‑1 (MQ=255)
                                                                      atGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATc                             <  2:3753110/99‑1 (MQ=255)
                                                                         gACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtgg                          >  2:2163290/1‑99 (MQ=255)
                                                                          aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggt                         >  1:219873/1‑99 (MQ=255)
                                                                          aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggt                         >  1:959350/1‑99 (MQ=255)
                                                                            aaTTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtgg                       <  2:3611941/99‑1 (MQ=255)
                                                                                ggTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                   <  2:1025778/99‑1 (MQ=255)
                                                                                  tttCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGtt                 <  1:1840356/99‑1 (MQ=255)
                                                                                  tttCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGtt                 >  1:88593/1‑99 (MQ=255)
                                                                                    tcttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTa               >  1:10645/1‑99 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa              <  1:618919/99‑1 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa              <  1:3351622/99‑1 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa              <  2:1402133/99‑1 (MQ=255)
                                                                                     cttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa              <  2:1552207/99‑1 (MQ=255)
                                                                                      ttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAc             <  1:2524750/99‑1 (MQ=255)
                                                                                      ttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAc             >  2:2587625/1‑99 (MQ=255)
                                                                                       tctcTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACg            <  1:2276413/99‑1 (MQ=255)
                                                                                        ctctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGc           <  1:2469746/99‑1 (MQ=255)
                                                                                         tctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc          >  2:3545102/1‑99 (MQ=255)
                                                                                         tctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc          <  2:2850272/99‑1 (MQ=255)
                                                                                           tGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGa        <  2:1102440/99‑1 (MQ=255)
                                                                                           tGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGa        <  2:3078429/99‑1 (MQ=255)
                                                                                           tGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGa        <  2:1745368/99‑1 (MQ=255)
                                                                                             aaTGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt      <  1:3694361/99‑1 (MQ=255)
                                                                                              aTGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTg     >  1:1424853/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
AGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGG  >  NC_000913/70273‑70468

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