Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,336,983 G→T 7.4% L17M (CTG→ATG) ‡ adiC ← arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,336,9830GT 7.4% 460.7 / 20.2 134L17M (CTG→ATG) ‡adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base G (62/62);  new base T (6/4);  total (68/66)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-03

Only 100 of 139 total aligned reads displayed.

TGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGT  >  NC_000913/4336713‑4337257
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                  
tGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGt                                                                                                                                                                                                                                                                                 >  2:255873/1‑274 (MQ=255)
                 ttGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACACCAGAATATTCCCCGACCCCATCAGGGGTACGGGGATTAAAcccc                                                                                                                                                                                                                                                                >  4:296083/1‑273 (MQ=255)
                             cagGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                                                                                                                                                                                                                                   >  2:70434/1‑275 (MQ=255)
                              agGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAACTGTCAATCCCCCCAGTAAAGGCCAGGTTTGCAGGCAACAGAAAAACCCCCTACCCCATAATGTTCCCCGACACCCTCAGGGTGACCGGAATTAATCCCAtttatgtacactca                                                                                                                                                                                                                                                  >  4:371112/1‑264 (MQ=255)
                                                         gATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGCCATAATTAAACTCCTGc                                                                                                                                                                                                                      >  2:289521/1‑275 (MQ=255)
                                                               cGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAgg                                                                                                                                                                                                                 >  4:37168/1‑275 (MQ=255)
                                                                             aagagggggcgggcgAAAGGGCACAAACCACACGGACTTGGAGCGTGAAACGACATTTTGCCGTAAACAATCAAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                                                                                                                                                                   <  4:255872/268‑1 (MQ=255)
                                                                                        ggcgTAAGCGTAAAAACCACCCGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTAACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                                                                                                                                                                         <  1:93342/274‑1 (MQ=255)
                                                                                          cataagcggaagACCGCCCAGGACTTGGGTCGAGGAACGACATTTTGGCGTAGACCATCGAGAGCGCCAGTGCACCGGTAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAGTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                                                                                                                                                                      <  2:515371/273‑1 (MQ=255)
                                                                                                  taagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                                                                                                                                                                                                                                           >  2:313849/1‑214 (MQ=255)
                                                                                                     aaaCCCCCCGGCCTTGGTTCGAGGAACGTAATTTTGGCGGATACAATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAACATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                                                                                                                                                            <  4:425450/273‑1 (MQ=255)
                                                                                                              gggctggggTCGAGGAACAACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGGCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                                                                                                                                                                  <  2:278330/269‑1 (MQ=255)
                                                                                                               gACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCAGCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                                                                                                                                                                  <  2:371113/274‑1 (MQ=255)
                                                                                                                   tGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATACTATTCCCCGACACAATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCCGAATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTCCGTCATATCATACGCCTTCCTTCGTTAttt                                                                                                                                                             >  4:278329/1‑275 (MQ=255)
                                                                                                                      gtggAGGAACGACATTTTGGCTTATACCATCGAAAACCCCAGTGCACCGGTAATCGTCACCAACCAGCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACACCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                                                                                                                                                          <  2:515346/272‑1 (MQ=255)
                                                                                                                      gtcggggacctccATTTTGGCTTAAACCCCCGAGAGCCCCAGTGAACCGATACTCGTCAGCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                                                                                                                                                          <  4:483483/263‑1 (MQ=255)
                                                                                                                          aGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCCGACCCCATAATATTCCCCGACACCATAAGGGTGACGGGGATTAAGCCCACTTAGCGAGCATCAGCATCCGAAGATATACTTAAACTAGAGAGAAGGCGCGAACCGTGACAATAAATTAAGTCATATCATACGCCTGAAGGAGCTTGTTGAattaa                                                                                                                                                      >  1:308/1‑274 (MQ=255)
                                                                                                                                 aCATTTTGGCGTATACCATCGAAAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAc                                                                                                                                               <  2:397769/275‑1 (MQ=255)
                                                                                                                                 aCATTTTGGCGGAACCACTCGGAGGCCCCAGGGACCCGATAATGGTCACCACCAACCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAAGATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAGTCGATAGTTGAATTATTATCGAc                                                                                                                                               <  4:424735/275‑1 (MQ=255)
                                                                                                                                  cATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCg                                                                                                                                                                                                                                                                             <  5:1073273/148‑1 (MQ=255)
                                                                                                                                                    tCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                                                                                                                                                                                                                                          <  7:182556/134‑1 (MQ=255)
                                                                                                                                                         agCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                                                                                                                                                          >  4:458846/1‑209 (MQ=255)
                                                                                                                                                         agCCCAAGTGCACCGGTAATCGTCACCAACCATCCATAAAGAGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGtta                                                                                                                         <  2:512432/274‑1 (MQ=255)
                                                                                                                                                          gCCCAAGTGCACCGATAATCGTCACCAACCCGCCAGAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat                                                                                                                        <  4:70433/274‑1 (MQ=255)
                                                                                                                                                            cccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                                                                                                                                                                    >  7:947755/1‑116 (MQ=255)
                                                                                                                                                            cccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                                                                                                                                                                    <  8:947755/116‑1 (MQ=255)
                                                                                                                                                             ccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                                                                                                                                                                      <  2:458847/209‑1 (MQ=255)
                                                                                                                                                                          aaccgtgcccaaccaTCCTAAACTGCCACTCCCGCCAGTAGAGGCCAGGTTTGCATGTACCAGAACAACGCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGATCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAAGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGtt                                                                                                       <  4:444624/267‑1 (MQ=255)
                                                                                                                                                                             cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgggg                                                                                                                                                                                                                                                                          >  6:377542/1‑109 (MQ=255)
                                                                                                                                                                             cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                                                                                                                                                   <  5:410024/148‑1 (MQ=255)
                                                                                                                                                                              gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCg                                                                                                                                                                                                                                                                             >  5:785728/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)
                                                                                                                                                                                         tccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGGCGGGGAt                                                                                                                                                                                                                                                                        <  7:1049377/99‑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)
                                                                                                                                                                                                   gCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                                                                                                                                                                       >  3:134737/1‑154 (MQ=255)
                                                                                                                                                                                                      aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                                                                                                                                                                    >  5:247601/1‑74 (MQ=255)
                                                                                                                                                                                                      aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                                                                                                                                                                    <  6:247601/74‑1 (MQ=255)
                                                                                                                                                                                                      aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                                                                                                                                                                                                                                        <  7:708752/86‑1 (MQ=255)
                                                                                                                                                                                                      aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                                                                                                                                                                                                                                        >  8:708752/1‑86 (MQ=255)
                                                                                                                                                                                                       tCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                                                                                                                                                                   <  1:134737/154‑1 (MQ=255)
                                                                                                                                                                                                        cccGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtg                                                                                                                                                                                                                                                           <  7:258144/97‑1 (MQ=255)
                                                                                                                                                                                                        cccGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtg                                                                                                                                                                                                                                                           >  8:258144/1‑97 (MQ=255)
                                                                                                                                                                                                          cGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                                                                                                                                                                                                                                          <  5:436256/80‑1 (MQ=255)
                                                                                                                                                                                                             cAGTAGAGGCCAGTTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCtgtg                                                                    <  1:93014/275‑1 (MQ=255)
                                                                                                                                                                                                             cAGTAGAGGCCAGGGTTGCAGGTAACAGAAAACCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGGGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCtgtg                                                                    <  4:355160/275‑1 (MQ=255)
                                                                                                                                                                                                              aGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                                                                                                                                                                                                                                         >  7:587351/1‑93 (MQ=255)
                                                                                                                                                                                                              aGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                                                                                                                                                                                                                                         <  8:587351/93‑1 (MQ=255)
                                                                                                                                                                                                                  gagGCCAGGTTTTCAGGTAACAGAAAAACACCTTACCCCACAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGAc                                                               <  4:26072/275‑1 (MQ=255)
                                                                                                                                                                                                                   agGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACGCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGAc                                                               <  2:144805/274‑1 (MQ=255)
                                                                                                                                                                                                                     gCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                                                                                                                                                           >  8:176260/1‑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)
                                                                                                                                                                                                                                  gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGa                                                                                                                                                                                                              >  6:462479/1‑116 (MQ=255)
                                                                                                                                                                                                                                  gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGa                                                                                                                                                                                                              <  5:462479/116‑1 (MQ=255)
                                                                                                                                                                                                                                    aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                                                                                                                                                   >  2:338825/1‑93 (MQ=255)
                                                                                                                                                                                                                                        gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATatt                                                                                                                                                                                                                                 <  6:135349/91‑3 (MQ=255)
                                                                                                                                                                                                                                        gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTATGAGCATCAGCATCCGCAGACATAAt                                                                                                                                                                                                                                 >  5:135349/1‑91 (MQ=255)
                                                                                                                                                                                                                                         aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGa                                                                                                                                                                                                                                          <  6:866447/81‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                                                                                                                                                   >  8:971679/1‑88 (MQ=255)
                                                                                                                                                                                                                                          aaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCGCTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAt                                                                                                                                                                      <  7:176260/148‑1 (MQ=255)
                                                                                                                                                                                                                                                 cTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAg                                  <  2:296084/273‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATTAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                                                                                                                                                              <  7:376775/148‑1 (MQ=255)
                                                                                                                                                                                                                                                   gACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                                                                                                                                                              <  6:638223/147‑1 (MQ=255)
                                                                                                                                                                                                                                                    aCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                                                                                                                                                                    <  7:981404/92‑1 (MQ=255)
                                                                                                                                                                                                                                                       ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                                                                                                                                                                      <  5:531420/103‑1 (MQ=255)
                                                                                                                                                                                                                                                        cATAATATTCCCCGACACCATCAGGGTGACGGGGATTa                                                                                                                                                                                                                                                                      >  5:918376/1‑38 (MQ=255)
                                                                                                                                                                                                                                                        cATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaa                                                                                                                                                                                                                              >  8:257897/1‑78 (MQ=255)
                                                                                                                                                                                                                                                         ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                                                                                                                                                                                                                                        >  8:134759/1‑67 (MQ=255)
                                                                                                                                                                                                                                                         ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                                                                                                                                                                                                                                        <  7:134759/67‑1 (MQ=255)
                                                                                                                                                                                                                                                         ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc                                                                                                                                                                                                                          <  8:808455/81‑1 (MQ=255)
                                                                                                                                                                                                                                                         ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                                                                                                                                                                     >  6:853461/1‑118 (MQ=255)
                                                                                                                                                                                                                                                         ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                                                                                                                                                                     <  5:853461/118‑1 (MQ=255)
                                                                                                                                                                                                                                                          taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                                                                                                                                                                    <  5:827488/86‑1 (MQ=255)
                                                                                                                                                                                                                                                          taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                                                                                                                                                                    >  6:827488/1‑86 (MQ=255)
                                                                                                                                                                                                                                                          taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                                                                                                                                                           >  7:886930/1‑111 (MQ=255)
                                                                                                                                                                                                                                                             taatCCCCGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                                                                                                                                                                    >  8:307170/4‑99 (MQ=37)
                                                                                                                                                                                                                                                                tCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCt                                                                                                                                   >  3:16927/1‑161 (MQ=255)
                                                                                                                                                                                                                                                                  cccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                                                                                                                                                                       <  7:769610/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                                                                                                                                                                        >  1:138384/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                   ccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagc                                                                                                                                                                                                                                                >  7:1019739/1‑49 (MQ=255)
                                                                                                                                                                                                                                                                   ccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagc                                                                                                                                                                                                                                                <  8:1019739/49‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                                                                                                                                                                        <  6:22930/87‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                                                                                                                                                                        >  5:22930/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                       cacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATNGACCGGGTTATCGCTGATGCCTTTATTAATATAGCTCTTCGTTACGGTTAAAAATCAGGCGACAGTTCGTTCTCTGTGCTNACTAACATAATGCAGCGGCCTCCTANCGAATAAGTCTTTAAAAAAAAAAATaaaa          >  3:85316/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                       cacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGAt                                                                                                                                >  4:296311/1‑157 (MQ=255)
                                                                                                                                                                                                                                                                           aTCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCg                                                                                                                            <  2:296312/157‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aTCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTacg      <  2:37169/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aTCAGAGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAAACATAATTAAACTCCTGCGAAGGCTAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTacg      <  1:124443/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                                                                                                                                                          >  6:555302/1‑78 (MQ=255)
                                                                                                                                                                                                                                                                            tCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                                                                                                                                                          <  5:555302/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAtt                                                                                                                                                                                            >  8:647692/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                             cAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAtt                                                                                                                                                                                            <  7:647692/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                  
TGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACGACGT  >  NC_000913/4336713‑4337257

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: