Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,336,979 A→C 8.0% M18R (ATG→AGG)  adiC ← arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,336,9790AC 8.0% 425.7 / 25.7 133M18R (ATG→AGG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (61/60);  new base C (8/4);  total (69/64)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

Only 100 of 138 total aligned reads displayed.

TACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACG  >  NC_000913/4336705‑4337253
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                  
tACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCa                                                                                                                                                                                                                                                                                    >  2:26073/1‑275 (MQ=255)
        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)
                                                                       cGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCNGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCAACAGGGTGACGGGGANAAAGCCCACTTTGTGAGCATCAACANCCGAAGAAAAAATAAAACAACTGCGaaaa                                                                                                                                                                                                             >  3:93342/1‑273 (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)
                                                                                                          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)
                                                                                                                              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)
                                                                                                                                                          aTCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACtt                                                                                                                                                                                                                                                         >  6:208470/1‑148 (MQ=255)
                                                                                                                                                            tCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                                                                                                                                                                                                                                      <  7:182556/134‑1 (MQ=255)
                                                                                                                                                            tCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                                                                                                                                                                                                                                      >  8:182556/1‑134 (MQ=255)
                                                                                                                                                                 agCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                                                                                                                                                      >  4:458846/1‑209 (MQ=255)
                                                                                                                                                                 agCCCAAGTGCACCGGTAATCGTCACCAACCATCCATAAAGAGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGtta                                                                                                                     <  2:512432/274‑1 (MQ=255)
                                                                                                                                                                  gCCCAAGTGCACCGATAATCGTCACCAACCCGCCAGAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat                                                                                                                    <  4:70433/274‑1 (MQ=255)
                                                                                                                                                                    cccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                                                                                                                                                                <  8:947755/116‑1 (MQ=255)
                                                                                                                                                                     ccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                                                                                                                                                                  <  2:458847/209‑1 (MQ=255)
                                                                                                                                                                                  aaccgtgcccaaccaTCCTAAACTGCCACTCCCGCCAGTAGAGGCCAGGTTTGCATGTACCAGAACAACGCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGATCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAAGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGtt                                                                                                   <  4:444624/267‑1 (MQ=255)
                                                                                                                                                                                  aaccGTCCCCACCCCTCCACAAAAGCCAACCCCGCCAGTAGAGGCCCGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGtt                                                                                                   <  4:23983/272‑1 (MQ=255)
                                                                                                                                                                                     cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgggg                                                                                                                                                                                                                                                                      >  6:377542/1‑109 (MQ=255)
                                                                                                                                                                                     cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                                                                                                                                               <  5:410024/148‑1 (MQ=255)
                                                                                                                                                                                      gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCg                                                                                                                                                                                                                                                                         <  6:785728/105‑1 (MQ=255)
                                                                                                                                                                                       tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCg                                                                                                                                                                                                                                                                         >  7:485879/1‑104 (MQ=255)
                                                                                                                                                                                       tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgg                                                                                                                                                                                                                                                                        <  7:657052/105‑1 (MQ=255)
                                                                                                                                                                                       tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAg                                                                                                                                                                                                                                                                >  8:800626/1‑113 (MQ=255)
                                                                                                                                                                                       tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGCGATTAAg                                                                                                                                                                                                                                                                <  7:800626/113‑1 (MQ=255)
                                                                                                                                                                                             accaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACg                                                                                                                                                                                                                                                                         <  6:51273/98‑1 (MQ=255)
                                                                                                                                                                                             accaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACg                                                                                                                                                                                                                                                                         >  5:51273/1‑98 (MQ=255)
                                                                                                                                                                                                 tccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGGCGGGGAt                                                                                                                                                                                                                                                                    <  7:1049377/99‑1 (MQ=255)
                                                                                                                                                                                                 tccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                                                                                                                                                                                                                                    >  8:1049377/1‑99 (MQ=255)
                                                                                                                                                                                                  ccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAg                                                                                    >  2:191698/1‑274 (MQ=255)
                                                                                                                                                                                                           gCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                                                                                                                                                                   >  3:134737/1‑154 (MQ=255)
                                                                                                                                                                                                            cAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTg                                                                           >  2:357303/1‑273 (MQ=255)
                                                                                                                                                                                                             aaTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCa                                                                                                                                                                                                                                                                                 <  5:418609/74‑1 (MQ=255)
                                                                                                                                                                                                             aaTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCa                                                                                                                                                                                                                                                                                 >  6:418609/1‑74 (MQ=255)
                                                                                                                                                                                                              aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                                                                                                                                                                <  6:247601/74‑1 (MQ=255)
                                                                                                                                                                                                              aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                                                                                                                                                                >  5:247601/1‑74 (MQ=255)
                                                                                                                                                                                                              aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                                                                                                                                                                                                                                    <  7:708752/86‑1 (MQ=255)
                                                                                                                                                                                                              aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                                                                                                                                                                                                                                    >  8:708752/1‑86 (MQ=255)
                                                                                                                                                                                                               tCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                                                                                                                                                               <  1:134737/154‑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                                                                                                                                                                                                                                                     <  8:587351/93‑1 (MQ=255)
                                                                                                                                                                                                                      aGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                                                                                                                                                                                                                                     >  7:587351/1‑93 (MQ=255)
                                                                                                                                                                                                                          gagGCCAGGTTTTCAGGTAACAGAAAAACACCTTACCCCACAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGAc                                                           <  4:26072/275‑1 (MQ=255)
                                                                                                                                                                                                                           agGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACGCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGAc                                                           <  2:144805/274‑1 (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)
                                                                                                                                                                                                                                          gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGa                                                                                                                                                                                                          >  6:462479/1‑116 (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)
                                                                                                                                                                                                                                                 aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                                                                                                                                               <  7:971679/88‑1 (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)
                                                                                                                                                                                                                                                           gACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                                                                                                                                                          >  5:638223/1‑147 (MQ=255)
                                                                                                                                                                                                                                                            aCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                                                                                                                                                                <  7:981404/92‑1 (MQ=255)
                                                                                                                                                                                                                                                               ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                                                                                                                                                                  >  6:531420/1‑103 (MQ=255)
                                                                                                                                                                                                                                                               ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                                                                                                                                                                  <  5:531420/103‑1 (MQ=255)
                                                                                                                                                                                                                                                                cATAATATTCCCCGACACCATCAGGGTGACGGGGATTa                                                                                                                                                                                                                                                                  >  5:918376/1‑38 (MQ=255)
                                                                                                                                                                                                                                                                cATAATATTCCCCGACACCATCAGGGTGACGGGGATTa                                                                                                                                                                                                                                                                  <  6:918376/38‑1 (MQ=255)
                                                                                                                                                                                                                                                                cATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaa                                                                                                                                                                                                                          <  7:257897/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                cATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaa                                                                                                                                                                                                                          >  8:257897/1‑78 (MQ=255)
                                                                                                                                                                                                                                                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                                                                                                                                                                                                                                    >  8:134759/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACAt                                                                                                                                                                                                                                <  7:1055777/71‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACAt                                                                                                                                                                                                                                >  8:1055777/1‑71 (MQ=255)
                                                                                                                                                                                                                                                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc                                                                                                                                                                                                                      <  8:808455/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc                                                                                                                                                                                                                      >  7:808455/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                                                                                                                                                                 >  6:853461/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                                                                                                                                                                 <  5:853461/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                  taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGa                                                                                                                                                                                                                <  5:827488/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                  taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                                                                                                                                                       >  7:886930/1‑111 (MQ=255)
                                                                                                                                                                                                                                                                        tCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCt                                                                                                                               >  3:16927/1‑161 (MQ=255)
                                                                                                                                                                                                                                                                          cccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                                                                                                                                                                    >  1:138384/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                          cccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATCCCGTAGAAGTCTTTAAAGACAGaa            >  2:290966/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                             gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                                                                                                                                                                    >  5:22930/1‑87 (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)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                  
TACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCTTCGTTAAGGTTAAAAATCAGGCGACAGTGCGTTCTCTGTGCTGACTAACATAATGCAGCGGCGTCATACCGTAGAAGTCTTTAAAGACAGAAATAAAGTACG  >  NC_000913/4336705‑4337253

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: