Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,336,974 A→T 8.3% S20T (TCG→ACG)  adiC ← arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,336,9740AT 8.3% 438.3 / 25.2 133S20T (TCG→ACG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (63/59);  new base T (7/4);  total (70/63)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-02

Only 100 of 138 total aligned reads displayed.

TACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAA  >  NC_000913/4336705‑4337243
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                             
tACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCa                                                                                                                                                                                                                                                                          >  2:26073/1‑275 (MQ=255)
            caccaTGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgg                                                                                                                                                                                                                                                              >  2:478553/1‑275 (MQ=255)
                         ttGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACACCAGAATATTCCCCGACCCCATCAGGGGTACGGGGATTAAAcccc                                                                                                                                                                                                                                                  >  4:296083/1‑273 (MQ=255)
                                     cagGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                                                                                                                                                                                                                     >  2:70434/1‑275 (MQ=255)
                                                                 gATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGCCATAATTAAACTCCTGc                                                                                                                                                                                                        >  2:289521/1‑275 (MQ=255)
                                                                       cGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTCAGAACCACCAGGACTGGGGGCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCATTGCACCGATCATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTACAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCCACACCCTCAGGTTGACGGGGATTCAGCCCACTTTGTGAGCATCAACCTCCCAAACCATATATAAACACCAGcacagg                                                                                                                                                                                                   >  4:515345/1‑270 (MQ=255)
                                                                       cGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAgg                                                                                                                                                                                                   >  4:37168/1‑275 (MQ=255)
                                                                                 gCCAAAGCAGCGGCGAGCGGAAGCTTAAGAACCACCAGGGTTTGGTTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                                                                                                                                                         <  2:533120/275‑1 (MQ=255)
                                                                                     aagagggggcgggcgAAAGGGCACAAACCACACGGACTTGGAGCGTGAAACGACATTTTGCCGTAAACAATCAAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                                                                                                                                                     <  4:255872/268‑1 (MQ=255)
                                                                                                ggcgTAAGCGTAAAAACCACCCGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTAACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                                                                                                                                                           <  1:93342/274‑1 (MQ=255)
                                                                                                  cgtaagcgtaagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAGCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                                                                                                                                                        <  2:310443/275‑1 (MQ=255)
                                                                                                          taagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                                                                                                                                                                                                                             >  2:313849/1‑214 (MQ=255)
                                                                                                             aaaCCCCCCGGCCTTGGTTCGAGGAACGTAATTTTGGCGGATACAATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAACATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                                                                                                                                              <  4:425450/273‑1 (MQ=255)
                                                                                                                      gggctggggTCGAGGAACAACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGGCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                                                                                                                                                    <  2:278330/269‑1 (MQ=255)
                                                                                                                        aCTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCACCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATCCGCCTGCATTCGATa                                                                                                                                                  >  2:372136/1‑275 (MQ=255)
                                                                                                                              gtggAGGAACGACATTTTGGCTTATACCATCGAAAACCCCAGTGCACCGGTAATCGTCACCAACCAGCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACACCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                                                                                                                                            <  2:515346/272‑1 (MQ=255)
                                                                                                                              gtcggggacctccATTTTGGCTTAAACCCCCGAGAGCCCCAGTGAACCGATACTCGTCAGCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                                                                                                                                            <  4:483483/263‑1 (MQ=255)
                                                                                                                                  aGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCCGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGCCATAATTAAACTCCTACGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat                                                                                                                                        >  1:94854/1‑275 (MQ=255)
                                                                                                                                  aGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCAccctga                                                                                                                                                                                                                                                                       >  8:1035649/1‑144 (MQ=255)
                                                                                                                                    gAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac                                                                                                                                                                                                                                                                              <  7:992259/139‑1 (MQ=255)
                                                                                                                                         aCATTTTGGCGTATACCATCGAAAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAc                                                                                                                                 <  2:397769/275‑1 (MQ=255)
                                                                                                                                                          aTCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACtt                                                                                                                                                                                                                                               >  6:208470/1‑148 (MQ=255)
                                                                                                                                                           tCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                                                                                                                                                                                                                            >  8:182556/1‑134 (MQ=255)
                                                                                                                                                                agCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat                                                                                                          >  2:23983/1‑275 (MQ=255)
                                                                                                                                                                agCCCAAGTGCACCGGTAATCGTCACCAACCATCCATAAAGAGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGtta                                                                                                           <  2:512432/274‑1 (MQ=255)
                                                                                                                                                                 gCCCAAGTGCACCGATAATCGTCACCAACCCGCCAGAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat                                                                                                          <  4:70433/274‑1 (MQ=255)
                                                                                                                                                                                 aaccgtgcccaaccaTCCTAAACTGCCACTCCCGCCAGTAGAGGCCAGGTTTGCATGTACCAGAACAACGCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGATCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAAGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGtt                                                                                         <  4:444624/267‑1 (MQ=255)
                                                                                                                                                                                 aaccGTCCCCACCCCTCCACAAAAGCCAACCCCGCCAGTAGAGGCCCGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGtt                                                                                         <  4:23983/272‑1 (MQ=255)
                                                                                                                                                                                    cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgggg                                                                                                                                                                                                                                                            <  5:377542/109‑1 (MQ=255)
                                                                                                                                                                                    cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgggg                                                                                                                                                                                                                                                            >  6:377542/1‑109 (MQ=255)
                                                                                                                                                                                    cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                                                                                                                                     <  5:410024/148‑1 (MQ=255)
                                                                                                                                                                                     gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCg                                                                                                                                                                                                                                                               >  5:785728/1‑105 (MQ=255)
                                                                                                                                                                                     gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCg                                                                                                                                                                                                                                                               <  6:785728/105‑1 (MQ=255)
                                                                                                                                                                                      tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCg                                                                                                                                                                                                                                                               >  7:485879/1‑104 (MQ=255)
                                                                                                                                                                                      tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgg                                                                                                                                                                                                                                                              <  7:657052/105‑1 (MQ=255)
                                                                                                                                                                                      tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgg                                                                                                                                                                                                                                                              >  8:657052/1‑105 (MQ=255)
                                                                                                                                                                                      tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAg                                                                                                                                                                                                                                                      >  8:800626/1‑113 (MQ=255)
                                                                                                                                                                                      tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGCGATTAAg                                                                                                                                                                                                                                                      <  7:800626/113‑1 (MQ=255)
                                                                                                                                                                                            accaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACg                                                                                                                                                                                                                                                               <  6:51273/98‑1 (MQ=255)
                                                                                                                                                                                                tccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac                                                                                                                                                                                                                                                                            <  7:422495/81‑1 (MQ=255)
                                                                                                                                                                                                tccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                                                                                                                                                                                                                          >  8:1049377/1‑99 (MQ=255)
                                                                                                                                                                                                 ccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAg                                                                          >  2:191698/1‑274 (MQ=255)
                                                                                                                                                                                                    tAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgg                                                                                                                                                                                                                                                              >  6:523073/1‑91 (MQ=38)
                                                                                                                                                                                                    tAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAg                                                                          >  2:343472/1‑271 (MQ=255)
                                                                                                                                                                                                        tGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgg                                                                                                                                                                                                                                                              >  7:503661/1‑87 (MQ=38)
                                                                                                                                                                                                           cAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTg                                                                 >  2:357303/1‑273 (MQ=255)
                                                                                                                                                                                                            aaTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCa                                                                                                                                                                                                                                                                       >  6:418609/1‑74 (MQ=255)
                                                                                                                                                                                                             aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                                                                                                                                                      >  5:247601/1‑74 (MQ=255)
                                                                                                                                                                                                             aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                                                                                                                                                                                                                          <  7:708752/86‑1 (MQ=255)
                                                                                                                                                                                                              tCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                                                                                                                                                     <  1:134737/154‑1 (MQ=255)
                                                                                                                                                                                                               cccGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtg                                                                                                                                                                                                                                             >  8:258144/1‑97 (MQ=255)
                                                                                                                                                                                                               cccGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtg                                                                                                                                                                                                                                             <  7:258144/97‑1 (MQ=255)
                                                                                                                                                                                                                 cGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                                                                                                                                                                                                                            <  5:436256/80‑1 (MQ=255)
                                                                                                                                                                                                                 cGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                                                                                                                                                                                                                            >  6:436256/1‑80 (MQ=255)
                                                                                                                                                                                                                    cAGTAGAGGCCAGTTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCtgtg                                                      <  1:93014/275‑1 (MQ=255)
                                                                                                                                                                                                                    cAGTAGAGGCCAGGGTTGCAGGTAACAGAAAACCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGGGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCtgtg                                                      <  4:355160/275‑1 (MQ=255)
                                                                                                                                                                                                                     aGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                                                                                                                                                                                                                           <  8:587351/93‑1 (MQ=255)
                                                                                                                                                                                                                         gagGCCAGGTTTTCAGGTAACAGAAAAACACCTTACCCCACAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGAc                                                 <  4:26072/275‑1 (MQ=255)
                                                                                                                                                                                                                            gCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                                                                                                                                             >  8:176260/1‑148 (MQ=255)
                                                                                                                                                                                                                              ctggTTTGCAGGTAACAAAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATTAGCATCCGAAGACATAATTGAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                                                                                                                                                           >  8:376775/3‑148 (MQ=255)
                                                                                                                                                                                                                              cAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAAc                                             >  2:457312/1‑274 (MQ=255)
                                                                                                                                                                                                                                 gTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCAGCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACAGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCAACAGTGCGTTCTCTTTGCTGcataacataa                                         >  4:120066/1‑275 (MQ=255)
                                                                                                                                                                                                                                      cAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGATACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGANAATAAATTACGTAATCTCATACGCCTGCATTCGATAGTTTACTTATTATCGACCGGGTTATCGCTGTTGNCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGAAAANANGNTATCTGCGCTGACTAATATAAtaaaa                                    >  3:54572/1‑271 (MQ=255)
                                                                                                                                                                                                                                      cAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATCAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCCACAGTGCGTTCTCTGTTCTGACTAACATACTCCAg                                    >  4:11804/1‑275 (MQ=255)
                                                                                                                                                                                                                                       aGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcat                                                                                                                                                                  <  5:208470/148‑1 (MQ=255)
                                                                                                                                                                                                                                         gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGa                                                                                                                                                                                                <  5:462479/116‑1 (MQ=255)
                                                                                                                                                                                                                                           aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAgcggcg                               <  4:457311/275‑1 (MQ=255)
                                                                                                                                                                                                                                               gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                                                                                                                                                                                                                 <  4:338824/93‑1 (MQ=255)
                                                                                                                                                                                                                                               gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTATGAGCATCAGCATCCGCAGACATAAt                                                                                                                                                                                                                   >  5:135349/1‑91 (MQ=255)
                                                                                                                                                                                                                                                aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGa                                                                                                                                                                                                                            <  6:866447/81‑1 (MQ=255)
                                                                                                                                                                                                                                                aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGa                                                                                                                                                                                                                            >  5:866447/1‑81 (MQ=255)
                                                                                                                                                                                                                                                aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                                                                                                                                     <  7:971679/88‑1 (MQ=255)
                                                                                                                                                                                                                                                aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                                                                                                                                     >  8:971679/1‑88 (MQ=255)
                                                                                                                                                                                                                                                 aaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCGCTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAt                                                                                                                                                        <  7:176260/148‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAg                    <  2:296084/273‑1 (MQ=255)
                                                                                                                                                                                                                                                          gACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                                                                                                                                                >  5:638223/1‑147 (MQ=255)
                                                                                                                                                                                                                                                           aCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                                                                                                                                                      <  7:981404/92‑1 (MQ=255)
                                                                                                                                                                                                                                                           aCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                                                                                                                                                      >  8:981404/1‑92 (MQ=255)
                                                                                                                                                                                                                                                              ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                                                                                                                                                        <  5:531420/103‑1 (MQ=255)
                                                                                                                                                                                                                                                              ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                                                                                                                                                        >  6:531420/1‑103 (MQ=255)
                                                                                                                                                                                                                                                               cATAATATTCCCCGACACCATCAGGGTGACGGGGATTa                                                                                                                                                                                                                                                        <  6:918376/38‑1 (MQ=255)
                                                                                                                                                                                                                                                               cATAATATTCCCCGACACCATCAGGGTGACGGGGATTa                                                                                                                                                                                                                                                        >  5:918376/1‑38 (MQ=255)
                                                                                                                                                                                                                                                               cATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaa                                                                                                                                                                                                                <  7:257897/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                                                                                                                                                                                                                          <  7:134759/67‑1 (MQ=255)
                                                                                                                                                                                                                                                                ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                                                                                                                                                                                                                          >  8:134759/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACAt                                                                                                                                                                                                                      <  7:1055777/71‑1 (MQ=255)
                                                                                                                                                                                                                                                                ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc                                                                                                                                                                                                            >  7:808455/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc                                                                                                                                                                                                            <  8:808455/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                                                                                                                                                       <  5:853461/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                 taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                                                                                                                                                      >  6:827488/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                 taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                                                                                                                                             >  7:886930/1‑111 (MQ=255)
                                                                                                                                                                                                                                                                 taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                                                                                                                                             <  8:886930/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                    taatCCCCGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                                                                                                                                                      >  8:307170/4‑99 (MQ=37)
                                                                                                                                                                                                                                                                       tCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCt                                                                                                                     >  3:16927/1‑161 (MQ=255)
                                                                                                                                                                                                                                                                         cccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                                                                                                                                                          >  1:138384/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                          ccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagc                                                                                                                                                                                                                                  >  7:1019739/1‑49 (MQ=255)
                                                                                                                                                                                                                                                                           cGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATg                                                                                                                 <  1:16927/161‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                                                                                                                                                          <  6:22930/87‑1 (MQ=255)
                                                                                                                                                                                                                                                                            gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATACTGCAGCGGCTTCATACCGTAAAAGTCTTTAAAGACAGaa  >  4:41272/1‑271 (MQ=255)
                                                                                                                                                                                                                                                                             acacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                                                                                                                                                      <  3:138384/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                             
TACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAA  >  NC_000913/4336705‑4337243

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