Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 10,336 A→G A53A (GCT→GCC)  satP ← acetate/succinate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091310,3360AG97.7% 132.2 / NA 44A53A (GCT→GCC) satPacetate/succinate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (25/18);  total (26/18)
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.91e-01

CAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAGCGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCAG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAG  >  NC_000913/10093‑10580
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
cAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑Caaaaa                                                                                                                                                                                                                                               <  1:423885/251‑1 (MQ=255)
   cacagTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAAttt                                                                                                                                                                                                                                            >  1:473971/5‑251 (MQ=255)
    acaggcAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCAGTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTg                                                                                                                                                                                                                                           <  1:260812/246‑1 (MQ=255)
         cAGGCTAAAGAGAAAGAGATGAAGAACGCGTGCGCCTTTCAGCGTGCAGAAGAACATAACCAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGCCCCAGTTTCGGCATCACCAGAATCGCAACCAACGTACGACAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTATTGCCTTTTTTGTATTCCAGCAGCCCGG‑CAAAAATTTgaccgc                                                                                                                                                                                                                                      >  2:63009/1‑246 (MQ=255)
                         aaTTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCTAGTTTCGGCATCAGCAGAATCGCAACAAACGTCAGCCAGAAAGAACCGTAAGAGGCGAATGCTGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAACATTTGCGCGATGCCGCCGTAGTAAAt                                                                                                                                                                                                                      >  2:30414/1‑251 (MQ=255)
                           tagCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACCTAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAAGTGAATGCGGTTAAACAGAAAGTGTTGCCTTTTTTGTACTCCAGCAAACCGG‑CACAAATTTGCGCTGTGCCGCCGTAgcatatgc                                                                                                                                                                                                                    >  2:93151/3‑244 (MQ=255)
                                 aacgcgtgaacatttaaGCGTGCCGAAGAACATAAACAGCGTAAATACGCCACACAGACCCAGGTAGCCACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGAATCAGCAGAATCGCAAACAACGTCAGCAAGAAAGAACCGTACGAGGTGAATGCGGTTAAAACGAAAGGATTGCCTTTTTTGTCATCCCGCAGACCGA‑CAACAATTTGCGCGCTGCCGCAGTACAAACTGCCAATgg                                                                                                                                                                                                              >  2:451252/17‑251 (MQ=255)
                                          gcCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGACGCCGTAGAACATGCCCATGGCAAGaata                                                                                                                                                                                                      >  1:24349/1‑250 (MQ=255)
                                                cAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGt                                                                                                                                                                                               <  1:362599/251‑1 (MQ=255)
                                                    gTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTcc                                                                                                                                                                                             >  1:242382/1‑249 (MQ=255)
                                                       ccGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGc                                                                                                                                                                                        <  1:350554/251‑1 (MQ=255)
                                                          aagaaCATAAACAGCGTAAATACGCCCCACAGACACAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGAAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCACGAATAATACCGTCCAGAGCGac                                                                                                                                                                                     >  2:451196/1‑250 (MQ=255)
                                                                        cGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACAGTCCAGAGCGAAATAACCCACGTtgg                                                                                                                                                                       >  2:162487/1‑250 (MQ=255)
                                                                             aTACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGGCCAGAGCGAAATAACCCACGTTGTGCAgg                                                                                                                                                                  <  1:406384/251‑1 (MQ=255)
                                                                                         gACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGACAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAAtg                                                                                                                                                      >  1:51742/1‑251 (MQ=255)
                                                                                                   gACACCAAGGAACGGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGc                                                                                                                                            <  2:28052/251‑1 (MQ=255)
                                                                                                    acacCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGATCCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGcc                                                                                                                                           <  1:433393/251‑1 (MQ=255)
                                                                                                        cAAGGAACTGTGCCTTTGGCGCATCGGACAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGa                                                                                                                                         >  1:396767/1‑249 (MQ=255)
                                                                                                         aaGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGc                                                                                                                                      >  1:426794/1‑251 (MQ=255)
                                                                                                                    cATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAACCCATCAGGccc                                                                                                                           >  2:191251/1‑251 (MQ=255)
                                                                                                                                    cAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGAtt                                                                                                             >  1:193961/1‑249 (MQ=255)
                                                                                                                                       aCCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTACGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGGGCCGGATTAGcc                                                                                                         >  2:481006/1‑250 (MQ=255)
                                                                                                                                           aGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTCAGAGGTGAATGCGGTTAAACCGAAAGTGTTGACTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGCATAATACCGTCAAGAGCGAAATAACCCACGTTGGGCAGGTTAAGCAGAATGTGGGTCATGCCGAAGCCCATCAGGCCCAGCGGGGCCGGATTAGCCAACtt                                                                                                    >  2:361302/1‑251 (MQ=255)
                                                                                                                                            gTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACtt                                                                                                    >  2:277159/1‑250 (MQ=255)
                                                                                                                                                   cATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGAGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGc                                                                                            >  1:317076/1‑251 (MQ=255)
                                                                                                                                                     tCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCg                                                                                                                 <  2:417014/228‑1 (MQ=255)
                                                                                                                                                     tCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCg                                                                                                                 >  1:417014/1‑228 (MQ=255)
                                                                                                                                                                gcaacaaacgccAGCCAGAAAGAACCTTAAGAGGTGAATGCGGTTAAACCAAATGGGTTGCCTTTTTTGTCCTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCa                                                                               <  2:266730/240‑1 (MQ=255)
                                                                                                                                                                    ccaacGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCAAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCaaaa                                                                            >  1:136866/5‑250 (MQ=255)
                                                                                                                                                                     caacGTCCGCCACAAAGAACCGGAAGAGGTGAAAGCGGATAATCCGAAAGTGTTGCCTTTTTTGTACCCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATTGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTTGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAat                                                                          <  2:39659/248‑1 (MQ=255)
                                                                                                                                                                      aacGTCAGCCAGAACGAACCGTAAGAGGTGAATGCGGTTCAACAGAAAGTGTTACCTTTTTTGTACTCCAGCAGAC‑AGTCAAAACTTTGCGCGATGCCGCCGTAGAAAATTCCCATGGCAAGAATAATACCGGCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGAGGGCATGACGAAACCCATCAGGCCCAGCGGTGCCGGAATTGCCAACTTAGTGTTGCCCATAATTCCTCAAACatc                                                                         >  2:433007/3‑251 (MQ=255)
                                                                                                                                                                         gccagacAGAAATAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCATTTTTTGTACTCCAGCAGACCGG‑CAAAAATGTGCGCGATGCCGCCGTAGAAAATGCCCATTGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGGTGCCCATAATTCCTCAAAAatcatc                                                                      <  2:66216/245‑1 (MQ=255)
                                                                                                                                                                            aGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAatcatcatc                                                                   >  1:313558/1‑251 (MQ=255)
                                                                                                                                                                               cAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCgaa                                                                <  1:150834/251‑1 (MQ=255)
                                                                                                                                                                                aaaaagaaCAGAAAGAGGTGAAGGCGGTTAAACCGAAAGTGTTGCCTTTTTGGTACTCAACCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCgaat                                                               <  2:23020/249‑1 (MQ=255)
                                                                                                                                                                                       aCCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGt                                                        <  1:214286/251‑1 (MQ=255)
                                                                                                                                                                                                ggTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAAtt                                               >  1:474362/1‑251 (MQ=255)
                                                                                                                                                                                                 gTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAGTGAATGGTGAAATAAttt                                              <  1:277159/251‑1 (MQ=255)
                                                                                                                                                                                                        cGGTTAAACCGAAAGTGGTGCATTTTTTGTCCTCCAGCAGACCGG‑CAAAGATTTGCGCGATGCCGCCGGAGAAAATGCCCAGGACAAGAATAATACCGTCCAGTACGAAATAACCCACGTTGTGCAGGTTAAACTGAATGGGGGACATACCGAAGCCCATCAGTCCCAGCGGGTCCGGATTGGCCTACTGAATGTTGCCCATAATTCCTCAAAAATCATCATCGCATGTATTGTGAAGTAATGTCCCTGa                                        >  1:100298/1‑250 (MQ=255)
                                                                                                                                                                                                             aaaCCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAAc                                   >  1:354937/1‑250 (MQ=255)
                                                                                                                                                                                                                     aGTGTTGCCTTTTTTGTACTCCAGCAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGtt                          >  1:388527/1‑251 (MQ=255)
                                                                                                                                                                                                                               tttttGTACCCCAGCAGACCGG‑CAAAAATTTGCGGGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGgc                  <  2:165705/249‑1 (MQ=255)
                                                                                                                                                                                                                                  tttaaCTCCAGAAGACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGcgcgc              <  2:388527/246‑2 (MQ=255)
                                                                                                                                                                                                                                             cagACCGG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAAGCATCATCGAAGGAATGGTGAAAGAATTTCCCTGAATAACAGTATTGTTTTCAGGGCGCGGCATAATAATCAg  <  1:170770/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
CAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAGCGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCAG‑CAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAG  >  NC_000913/10093‑10580

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