Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,559 A→G 5.4% intergenic (‑382/+15) asnB ← / ← umpH asparagine synthetase B/UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,5590AG 5.4% 239.2 / 3.6 74intergenic (‑382/+15)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base G (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-02

ATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCG  >  NC_000913/699303‑699826
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                           
aTTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAg                                                                                                                                                                                                                                                           >  2:302110/1‑275 (MQ=255)
                             atatTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTACTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTa                                                                                                                                                                                                                              >  2:88108/1‑275 (MQ=255)
                                                                                   caCGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATACATTGCAGGAGCTGCGTAGGCTTGATAAGCGTAGCGCATCAGGCAGTTTGTCGTTTGTCATCAGAGTCAACCATGTCCGCAGACGTGGTTCCTATTCAGATAACGTTGATTTTCGCGACCGACGGGTAAATCCAGCTTGGGCCGAAAGGCCTACTGTCGCTATCGTCGACTTACGAAacaac                                                                                                                                                                        >  4:460156/2‑273 (MQ=25)
                                                                                    aCGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGNATCAGGCAGTTTGGCGGTTGTCATCAGAGCCAACCCCGTCCGCAGACGTGGTGGCTATCCAAATAACGTNGATTTCCGCGACTGCCGGGTCAATCAAGCCGGGGCGGACACGCAANCNGTCGACAACGCCGACAGACAAAAAAcc                                                                                                                                                                        >  3:47014/1‑274 (MQ=25)
                                                                                         ggCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGATGCGTAGGCCTGAAAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGACCGCAGAAGTGGTTGCTATTCAGATAACGTCAATTTCAGGGACTGACGGGGAAATCCAGATGGGGCGGAAAGGCATACTGTCTAAATCGTGTAGGAACGaaaaccccaaaa                                                                                                                                                                   >  4:88873/1‑265 (MQ=25)
                                                                                                                  ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTg                                                                                                                                                                                                                                                                       <  7:121240/149‑1 (MQ=255)
                                                                                                                     ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGtt                                                                                                                                                                                                                                                                    <  7:384867/149‑1 (MQ=255)
                                                                                                                         tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTa                                                                                                                                                                                                                                                                <  7:185690/149‑1 (MQ=255)
                                                                                                                          gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGACCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGtt                                                                                                                                                                                                                                                                    >  8:6038/1‑144 (MQ=255)
                                                                                                                          gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGACCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGtt                                                                                                                                                                                                                                                                    <  7:6038/144‑1 (MQ=255)
                                                                                                                           ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTAtt                                                                                                                                                                                                                                                              >  6:293294/1‑149 (MQ=255)
                                                                                                                              gATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACCGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGc                                                                                                                              <  2:88874/274‑1 (MQ=255)
                                                                                                                                     ggttaaagCCTGATCCGGCATAAAAAAGCATGAAAATTTAATAAAGTGCAGGGGCTGCGTAGGGCTGATAAGCGTAGTGCATCAGGAAGTTTGGCGTTTGGCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAatat                                                                                                                      <  1:116363/268‑1 (MQ=255)
                                                                                                                                           cGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGAttt                                                                                                                                                                                                                                              <  5:293294/149‑1 (MQ=255)
                                                                                                                                              cTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCa                                                                                                                                                                                                                                            <  7:464234/148‑1 (MQ=255)
                                                                                                                                              cTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCa                                                                                                             >  1:14507/1‑275 (MQ=255)
                                                                                                                                              cTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAAGCAGTTTGACGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCa                                                                                                                                                                                                                                            >  8:464234/1‑148 (MQ=255)
                                                                                                                                               ttATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGc                                                                                                                                                                                                                                          >  5:128632/1‑149 (MQ=255)
                                                                                                                                                     ggCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGa                                                                                                                                                                                                                                    >  8:138046/1‑149 (MQ=255)
                                                                                                                                                            aaaaGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGc                                                                                                                                                                                                                                          >  5:99760/1‑136 (MQ=255)
                                                                                                                                                            aaaaGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGc                                                                                                                                                                                                                                          <  6:99760/136‑1 (MQ=255)
                                                                                                                                                            aaaaGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTaa                                                                                                                                                                                                                             >  5:398151/1‑149 (MQ=255)
                                                                                                                                                              aaGCATGAAAATTCAATATATTGCAGTAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGCCGATTCCAGCGACTGACGGGTAAATCCAGCTGGGGCGGTAAGGCATACTGTCGATATCGTCCAGCGCCGAAACCCCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAGCCCGGCCCCAATATCGGACCCCCGTTTATCGCCTACaa                                                                                              >  4:117886/1‑274 (MQ=255)
                                                                                                                                                               aGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc                                                                                                                                                                                                                                                             >  6:224731/1‑114 (MQ=255)
                                                                                                                                                               aGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc                                                                                                                                                                                                                                                             <  5:224731/114‑1 (MQ=255)
                                                                                                                                                                  aTGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAg                                                                                                                                                                                                                       >  7:52758/1‑149 (MQ=255)
                                                                                                                                                                    gAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCt                                                                                                                                                                                                                     >  7:41681/1‑149 (MQ=255)
                                                                                                                                                                    gAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCt                                                                                                                                                                                                                     >  7:147637/1‑149 (MQ=255)
                                                                                                                                                                     aaaaTTAAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCg                                                                                                                    <  4:128125/245‑1 (MQ=255)
                                                                                                                                                                            aaTATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCTAGCTGGGGCGGa                                                                                                                                                                                                             >  5:376140/1‑149 (MQ=255)
                                                                                                                                                                                tatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACt                                                                                                                                                                                                                                      <  6:387258/120‑1 (MQ=255)
                                                                                                                                                                                tatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACt                                                                                                                                                                                                                                      >  5:387258/1‑120 (MQ=255)
                                                                                                                                                                                 atTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGc                                                                                                                                                                                                        <  6:398151/149‑1 (MQ=255)
                                                                                                                                                                                 atTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGc                                                                                                                                                                                                        >  6:208944/1‑149 (MQ=255)
                                                                                                                                                                                   tGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACggg                                                                                                                                                                                                                                >  3:34369/1‑123 (MQ=255)
                                                                                                                                                                                    gCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATa                                                                                                                                                                                                     <  5:106746/149‑1 (MQ=255)
                                                                                                                                                                                       ggAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTaaa                                                                                                                                                                                                                            <  1:34369/123‑1 (MQ=255)
                                                                                                                                                                                        gAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAg                                                                                                                                                                                                                       >  5:43591/1‑127 (MQ=255)
                                                                                                                                                                                        gAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAg                                                                                                                                                                                                                       <  6:43591/127‑1 (MQ=255)
                                                                                                                                                                                          gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCAt                                                                  >  2:423098/1‑274 (MQ=255)
                                                                                                                                                                                                  ggCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtc                                                                                                                                                                                       <  5:208944/149‑1 (MQ=255)
                                                                                                                                                                                                   gCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtcg                                                                                                                                                                                      >  8:413067/1‑149 (MQ=255)
                                                                                                                                                                                                      tGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGc                                                                                                                                                                                   <  8:52758/149‑1 (MQ=255)
                                                                                                                                                                                                      tGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGc                                                                                                                                                                                   <  8:147637/149‑1 (MQ=255)
                                                                                                                                                                                                        aTAAGCGTAGCGCCTCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGCAACACCACGACGCACCAGAATCGTCCACAGACCTGCATGGCAGTCGGCCAGACTATCGGTACGCAGGTTATCGCCTCCCATCACCGTTCCTTCCGAATGCTCCTGCATTTAGTTTACtgctc                                                    >  4:77278/1‑273 (MQ=255)
                                                                                                                                                                                                             cGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATTCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCAAGAATAACGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGACTGCATTTTGTTTAATGCTGCGCGga                                              >  2:475594/1‑275 (MQ=255)
                                                                                                                                                                                                               tagtgaaTCAGGCAGTTTGGCGTTGGTCAGCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGgat                                             <  2:257316/268‑1 (MQ=255)
                                                                                                                                                                                                               tAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAaca                                                                                                                                                                          <  6:128632/149‑1 (MQ=255)
                                                                                                                                                                                                                 gcgcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACAcc                                                                                                                                                                        >  8:239094/1‑149 (MQ=255)
                                                                                                                                                                                                                  cgcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCACAACTACGGTAACAACGTTAACGCCTACAATCACCGTTTACTCCGACTCCGCATGCATTTTGTTTACTGCTGGGCggttt                                            >  4:388439/1‑269 (MQ=255)
                                                                                                                                                                                                                   gcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACAACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATcc                                        <  2:460161/275‑1 (MQ=255)
                                                                                                                                                                                                                       cAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGc                                                                                                                                  >  4:318718/1‑181 (MQ=255)
                                                                                                                                                                                                                       cAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACgg                                    >  2:318124/1‑275 (MQ=255)
                                                                                                                                                                                                                        agtgagtttggagtttggCATCAGAGCCAACAACGTCCGCAGACGTGGTTGCTATTCAGATAACGTGGATTTCAGCTACTGACGGGTCAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACggg                                   <  1:47014/271‑1 (MQ=255)
                                                                                                                                                                                                                          gCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCAc                                                                                                                                                               >  7:86932/1‑149 (MQ=255)
                                                                                                                                                                                                                           cAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCAcc                                                                                                                                                              <  7:239094/149‑1 (MQ=255)
                                                                                                                                                                                                                           cAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGc                                                                                                                              <  2:318722/181‑1 (MQ=255)
                                                                                                                                                                                                                           cAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTgcgc                                                 <  2:499097/258‑1 (MQ=255)
                                                                                                                                                                                                                                 ggCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATc                                                                                                                                                        <  8:86932/149‑1 (MQ=255)
                                                                                                                                                                                                                                 ggCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATc                                                                                                                                                        <  8:41681/149‑1 (MQ=255)
                                                                                                                                                                                                                                    gTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGCAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAACCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCa                       <  2:77279/275‑1 (MQ=255)
                                                                                                                                                                                                                                           atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTggg                                                                                                                                                                                                                  <  6:314985/81‑1 (MQ=255)
                                                                                                                                                                                                                                           atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTggg                                                                                                                                                                                                                  >  5:314985/1‑81 (MQ=255)
                                                                                                                                                                                                                                           atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAc                                                                                                                                              >  8:155091/1‑149 (MQ=255)
                                                                                                                                                                                                                                            tcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAcc                                                                                                                                             >  7:63263/1‑149 (MQ=255)
                                                                                                                                                                                                                                            tcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAcc                                                                                                                                             >  5:371254/1‑149 (MQ=255)
                                                                                                                                                                                                                                            gcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGAATACAGCGACTGACGGGTACATCCAGCTGTGGCGGAAAGGCATACTGTCCATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAAACGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGaa               <  2:129664/274‑1 (MQ=255)
                                                                                                                                                                                                                                                 gCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGa                                                                                                                                                                                     >  7:317598/1‑104 (MQ=255)
                                                                                                                                                                                                                                                 gCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGa                                                                                                                                                                                     <  8:317598/104‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGAcctgcctg                                                                                                                                       <  7:413067/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaCCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGc                                                                                                                                                                                   <  8:103435/103‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaCCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGc                                                                                                                                                                                   >  7:103435/1‑103 (MQ=255)
                                                                                                                                                                                                                                                          gTCCGCAGACGTGGTTGCTCTTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATACCGGTACGCAGGTTATAGGCGACACTCACCGTTTCCTCCGACGCCCCCTGCATTTTGTTTAATGCTGCCCGGTGGATCCACGGGATCGGCTTACTCACATACAACGTTTTTCCCCCg  >  4:263519/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               cAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGGCTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATACCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTgcgc                                                 <  2:86899/222‑1 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                           
ATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCAGAATATCGGTACGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCG  >  NC_000913/699303‑699826

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