Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 76,727 G→A H191H (CAC→CAT sgrR ← DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091376,7270GA97.9% 133.2 / ‑4.2 48H191H (CAC→CATsgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (17/30);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GAATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACATCGTCCATTTCCAGTTCACGACCATGGTGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGAAATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCGATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGA  >  NC_000913/76477‑76973
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
ggaTCACCGCATACGGTCCGGTGCCGATGGGATAGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGGCCATTTCCAGTTCACGACCATgg                                                                                                                                                                                                                                                        >  1:390768/3‑251 (MQ=255)
          aTACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGa                                                                                                                                                                                                                                              <  1:55675/251‑1 (MQ=255)
            aCGGACCGGTGCCGGTGGTATTGCTGGCAAAGTTACTGAGGGTTTCCCAATCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACTTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACt                                                                                                                                                                                                                                            <  1:22310/251‑1 (MQ=255)
                  cGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGa                                                                                                                                                                                                                                      <  1:15395/251‑1 (MQ=255)
                        cGATGGGATGTCTGTAAAAGTTACTGTGGGTATCCCCTTCGCGGTGCAGGAAACTTGCCAGATCTTGCCCCTGCATTAACTGTAACCAGCGGTCCGGTTGCGTGATATGGATATCCAGCGTCCAGGGCGTCGGCGACACATTTTCAGCAATCTGCGAATAGATCGTCAGCGTATTGATTCGTTTTAAAGAGGCGATCTCGTCGTCCATTTCCAGTTCACGACCATGATGAAAATTGAcctttgtacgcaaa                                                                                                                                                                                                                                <  2:358390/251‑14 (MQ=255)
                          atgggatggCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGCTTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCaaaa                                                                                                                                                                                                                               >  1:3380/1‑250 (MQ=255)
                            gggatggCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGa                                                                                                                                                                                                                            <  1:279018/251‑1 (MQ=255)
                               atggCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCCGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGCTATCCAGCGCCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGAGATCCCGTCGTCCCTTTCCAGTTCACGACCATAATGAAAAAGGAATCCTTGACGCAAAAAGGAAc                                                                                                                                                                                                                         >  2:331804/1‑251 (MQ=255)
                                                tGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGat                                                                                                                                                                                                        >  1:361705/1‑250 (MQ=255)
                                                  aGGGTTTCCCATTCACGCGGCAGGATCATCGCCCGAACTTGCCCCAGCAATAACGGTAACCAGCGTTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGCCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGTACCCGTCGACCGTTTCACGTTCACGACCTTGATTAAAATGGAATACTGGACGCAAACAGCAAAGCCAGTGAAGCggagtaat                                                                                                                                                                                                      >  2:391418/1‑245 (MQ=255)
                                                         cccATTCGCGCGGCAGGATCATCGCCGAAACCTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGGATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGACAGTGAAGCGGAGATATTtgctgc                                                                                                                                                                                               >  1:130720/1‑251 (MQ=255)
                                                          ccATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTATCCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGcc                                                                                                                                                                                              <  1:470950/251‑1 (MQ=255)
                                                              gcgcgcgGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAgt                                                                                                                                                                                           <  1:419916/249‑1 (MQ=255)
                                                                  gcgGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGtggtgg                                                                                                                                                                                      <  1:393026/251‑1 (MQ=255)
                                                                  gcgGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGtggtgg                                                                                                                                                                                      <  2:361705/251‑1 (MQ=255)
                                                                      cAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACCCAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTTGGCg                                                                                                                                                                                   >  2:283805/1‑250 (MQ=255)
                                                                            caccgccgaaccTTGCCCCATCAATAACGGGAACCAGCGGTCCTGTTGCTTGCGATGGATATCCAGCGTCCAGGGCGTCGGCGACTCAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTGTAAAGAGGCGATCACGTCGTCCATTTCCAGAGCTCGACCATGATGAAAATGGTCGTCTGTACGCAAAAAGAAACGCCGGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTg                                                                                                                                                                            <  2:281737/240‑1 (MQ=255)
                                                                                   ggAACTTGCCCCAGCAGTAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACATCGTCCATTTCCAGTTCACGACCATGGTGAAAATGAACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGAAATTTGCTGCCAGTGGTGGGCGAGGTCTGCTTcc                                                                                                                                                                       >  2:183770/1‑249 (MQ=255)
                                                                                     aaCTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTc                                                                                                                                                                                                                                           >  1:289776/1‑179 (MQ=255)
                                                                                     aaCTTGCCCCAGCAAGAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCATCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTc                                                                                                                                                                                                                                           <  2:289776/179‑1 (MQ=255)
                                                                                           ccccAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCAt                                                                                                                                                             <  1:159434/251‑1 (MQ=255)
                                                                                             ccAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCAttt                                                                                                                                                           <  1:294943/251‑1 (MQ=255)
                                                                                               agcaaTAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTc                                                                                                                                                         >  2:178345/1‑251 (MQ=255)
                                                                                                  aatAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTc                                                                                                                                                      <  1:377362/249‑1 (MQ=255)
                                                                                                           aaCCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATgcg                                                                                                                                             <  1:333349/251‑1 (MQ=255)
                                                                                                           aaCCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATgcg                                                                                                                                             >  1:94464/1‑251 (MQ=255)
                                                                                                                             gTGAGATGTATATCCAGCGTCCAGGTCGTCGGCGACACAATGACAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGAt                                                                                                                           <  1:251846/251‑1 (MQ=255)
                                                                                                                                    ggATATCCAGCGGCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCgg                                                                                                                    <  1:346208/251‑1 (MQ=255)
                                                                                                                                                 tccgggCGTCGGCGACACAATGTCCGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCTTGTCGGCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAATAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATGTATGCGCGTTAGCGAACTGATTATTTGCCGGGCTATATGGGtt                                                                                                        <  2:339837/247‑1 (MQ=255)
                                                                                                                                                 ccAGGGCGTCGGCGACACAATGTCAGCAATATGCGCATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCCGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTACTTCCAGTTCCCCATTTTCCTCCTTTATGCGCGTTAGCGACCTGAAGATTTGCCGGGCAATATGGGttt                                                                                                       >  2:361037/1‑251 (MQ=255)
                                                                                                                                                  cagggcgtagtgcacctcaggtaaGCAATATGCGAATAGAGCGGCAACGTATTGATTCTTTTTAAAGAGGCGACCACGTCGTTCATTCCCAGATCACGACCATGATGAAAATTGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTTCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCCTTTTCCTCATTGATGCGCGTTAGCAAACTGAAGATTTGCCGGCCAATATGGGttgc                                                                                                      <  2:409853/228‑3 (MQ=255)
                                                                                                                                                   aggtcGTCGGCGAAATAATGTAAGCAATATGGGAATAGGGCGGCCGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAATAAACGCCAGTGAATCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTGGCCTGGCAATTTGGGTTTCg                                                                                                     <  2:209247/247‑1 (MQ=255)
                                                                                                                                                        gTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACg                                                                                                <  1:285955/251‑1 (MQ=255)
                                                                                                                                                                    aTGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGGTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGGCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGTCGGGCAATAGGGGTTGCGGAACGGCGCAATGCGCt                                                                                    <  1:363977/251‑1 (MQ=255)
                                                                                                                                                                         aGCAATATGCGCATAGAGCGGCAGCGTATTGATTAGTTTTAACGAGGCGATCACGTCGTCCATTTCCAGTTCACGCCCATGATGAAAATGGTCTCCTGGACGCACAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGCTTTGCCGGGCAATATGGGTTTCGGCACGGCGCACTGCGCTGcc                                                                                 >  2:64232/1‑249 (MQ=255)
                                                                                                                                                                          gCAATATGCGCATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAATTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAAGTCGCTGCCAgg                                                                              >  2:440658/1‑251 (MQ=255)
                                                                                                                                                                          gCAATATGCGAATAGAGCGGCAGCTGATTGATGAGTTTTAAAGAGGCGATCCCGTCGGCCAGTGCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGcctgg                                                                              <  2:3380/251‑4 (MQ=255)
                                                                                                                                                                            aaTATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTa                                                                            <  1:112855/251‑1 (MQ=255)
                                                                                                                                                                                    aaTAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGAtt                                                                     >  1:75692/1‑250 (MQ=255)
                                                                                                                                                                                                       gattgggtttAAAGAGGCGTTCACGTCGTCCTATGCCAGTTTACGACCATGATGAAAAGGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGGGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTGTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTaga                                                 <  2:306916/244‑1 (MQ=255)
                                                                                                                                                                                                          tgtttttAAAGAGGCGATCACGTCGTCCTGTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGc                                               <  2:345590/247‑1 (MQ=255)
                                                                                                                                                                                                             gtttAAAGAGGCGATCAAGTCGTCAATTTCAAGTTCACGACCATGAGGAAAATGGACTCCTGGGCGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCACGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTGGCGAACGGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCGAGGTAGCAGCTAACGCAACGGACGATAATAGAGCAcgc                                           >  2:231564/2‑251 (MQ=255)
                                                                                                                                                                                                                 aaaGAGGCGATCACGTCGTCCATTTCCATTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGTCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCa                                        <  1:168711/250‑1 (MQ=255)
                                                                                                                                                                                                                         gatcaagtCGTCCATTTCAAGTTCACGACCATGATGAATATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTACTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTgcc                               >  1:123501/8‑251 (MQ=255)
                                                                                                                                                                                                                                cgctCCATTTCAAGTTCACTACCATGATGAAAATGGACTCTTGGTCGCAAATAGAAACGCCAGTTAAGCGGGGGTATTGGCGTCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGAGTTAGCGAACTGAAGATTTGCCGGGCAATAGGGGGTTCGGATCTGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCGCGCGCAGGATGTGCCGCCCCt                         <  2:455905/247‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTgg                       >  1:323409/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCg                     <  1:178345/250‑1 (MQ=255)
                                                                                                                                                                                                                                               tCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGcc         >  1:447162/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      cATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTTATATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCATTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCAATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATga  <  1:19379/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
GAATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACATCGTCCATTTCCAGTTCACGACCATGGTGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGAAATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCATTTTCCTCATTTATGCGCGTTAGCGAACTGAAGATTTGCCGGGCGATATGGGTTTCGGAACGGCGCAATGCGCTGCCAGGTAGCAGATTACGCAACGGACGATAGTAGAGCACGCGCAGGATGTGCCGCCCCTGGCGGAAGCTGCGGCCCAGATGA  >  NC_000913/76477‑76973

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: