Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,928,799 A→C 5.9% L14R (CTT→CGT)  ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,928,7990AC 5.9% 304.6 / 9.2 91L14R (CTT→CGT) ptrBprotease II
Reads supporting (aligned to +/- strand):  ref base A (37/48);  new base C (4/2);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

ATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCA  >  NC_000913/1928530‑1929058
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                   
aTGCCGATTGAGGCGGGTAAATGGCATATTCAAAGCCGGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTc                                                                                                                                                                                                                                                                <  4:395638/275‑1 (MQ=255)
aTGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTc                                                                                                                                                                                                                                                                <  2:412996/275‑1 (MQ=255)
           cGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAg                                                                                                                                                                                                                                                                    >  2:395643/1‑260 (MQ=255)
                   gATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcggg                                                                                                                                                                                                                                             <  1:140746/275‑1 (MQ=255)
                                  gCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCg                                                                                                                                                                                                                                                          >  4:236042/1‑247 (MQ=255)
                                   ccTGGTTCATAAATATGGCGGTAGCGGTAGCCATTTTTGAGTGAGGGCGCAGAACCTGTTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGTTCCTGCAAGGCTTGTGGTGATGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTAGCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCa                                                                                                                                                                                                                              <  4:16173/274‑1 (MQ=255)
                                      ggTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTggg                                                                                                                                                                                                                                                      <  2:236046/247‑1 (MQ=255)
                                       gTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTTGCAGGttc                                                                                                                                                                                                                         >  2:67684/1‑274 (MQ=255)
                                             aaaTATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTCACTATCGATGCGCGTACGGTCTCGAAGCTTCATGGCGCGGGGATTGCGGGCGGCTTTTGGTGGCATGTTATTGttc                                                                                                                                                                                                                   >  1:74094/1‑275 (MQ=255)
                                                  tGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGttctttct                                                                                                                                                                                                              >  2:513614/1‑275 (MQ=255)
                                                          gaggTAGCCATTTTTGAGGTAGGTCGCAGAAACTTCCCGTTCGGGAATGCGGTCGATGATTTCCTTTAGGTTTCGACCCTTCAAGGCTTGTTGCGGGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGGAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGaaa                                                                                                                                                                                                      <  4:67683/273‑1 (MQ=255)
                                                          gCGGTAGCCAGTTTTGATGTAGGGCGCAGAAACTGCTCGTGGCGGAGGGCGGTGGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGCCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACACGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGaaa                                                                                                                                                                                                      <  2:422893/275‑1 (MQ=255)
                                                           cGGTAGCCTGTTTTGATGTAGGGCGAGTAAACTTCTCTTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTGGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAAc                                                                                                                                                                                                     <  4:493306/275‑1 (MQ=255)
                                                            ggTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACa                                                                                                                                                                                                    <  2:69853/275‑1 (MQ=255)
                                                                            tgaGGGCGCAGAAACTTCTCGTGGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGGCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGc                                                                                                                                                                                     <  1:55953/272‑1 (MQ=255)
                                                                                                         tgcggtggtggttttCCTTTAAGATGCGTTCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACGATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAgg                                                                                                                                                       <  4:5907/263‑1 (MQ=255)
                                                                                                         tGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGaaa                                                                                                                                                           >  2:5907/1‑271 (MQ=255)
                                                                                                                          tttAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAg                                                                                                                                                                                                                                                                    >  5:391427/1‑149 (MQ=255)
                                                                                                                             aaGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGt                                                                                                                                                                                                                                                                 <  5:175111/149‑1 (MQ=255)
                                                                                                                                   cGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCGTCCTGGCGTGGGGAATGCGGGCTGCTTTTGGTAGCATGTTCTTGTTCTTTCTGGTTGCAACCTCTTATAAGGGTGGCAAAACTTCCCGTGGCTGCAAGCGAACCAGGGGAATCCTTGCCTAATTTGTGACAgg                                                                                                                             >  4:470175/1‑275 (MQ=255)
                                                                                                                                           cAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGaa                                                                                                                                                                                                                                                   <  7:72426/149‑1 (MQ=255)
                                                                                                                                             aGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATg                                                                                                                                                                                                                                                 <  6:116305/149‑1 (MQ=255)
                                                                                                                                                 ttgttgTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcggg                                                                                                                                                                                                                                             >  7:249968/1‑149 (MQ=255)
                                                                                                                                                  tgttgtGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCaa                                                                                                                                                                 >  4:343450/1‑224 (MQ=255)
                                                                                                                                                      gtgAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGa                                                                                                                                                             <  2:343454/224‑1 (MQ=255)
                                                                                                                                                       tgAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtt                                                                                                                                                                                                                                       >  5:372842/1‑149 (MQ=255)
                                                                                                                                                         aGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtttt                                                                                                                                                                                                                                     <  8:349924/149‑1 (MQ=255)
                                                                                                                                                         aGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtttt                                                                                                                                                                                                                                     <  8:13042/149‑1 (MQ=255)
                                                                                                                                                         aGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtttt                                                                                                                                                                                                                                     <  8:57100/149‑1 (MQ=255)
                                                                                                                                                           gCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTgg                                                                                                                                                                                                                                   >  5:103059/1‑149 (MQ=255)
                                                                                                                                                              aTCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGc                                                                                                  >  2:52341/1‑275 (MQ=255)
                                                                                                                                                              aTCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGGTATGc                                                                                                  >  2:101722/1‑275 (MQ=255)
                                                                                                                                                                cACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGc                                                                                                                                                                                                                                                           <  7:286416/120‑1 (MQ=255)
                                                                                                                                                                cACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGc                                                                                                                                                                                                                                                           >  8:286416/1‑120 (MQ=255)
                                                                                                                                                                      aTGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGa                                                                                                                                                                                                        >  2:498040/1‑165 (MQ=255)
                                                                                                                                                                      aTGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATTATGACAGGTAATGGCCCGGATATCGAATGTTATGCAGACagag                                                                                          >  2:31676/1‑275 (MQ=255)
                                                                                                                                                                      aTGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGcagacaga                                                                                           <  2:470180/274‑1 (MQ=255)
                                                                                                                                                                      aTGACCGTAACTAGTTGTTGGTGGCAGGTAGTCCAGTACGTCTGGCTGAGAAAGTGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCAAGGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGAAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGcagacaga                                                                                           <  2:280290/274‑1 (MQ=255)
                                                                                                                                                                       tGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTAtt                                                                                                                                                                                                                       <  6:267595/149‑1 (MQ=255)
                                                                                                                                                                         aCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGt                                                                                                                                                                                                                     >  6:249253/1‑149 (MQ=255)
                                                                                                                                                                          ccGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACa                                                                                                                                                                                                    <  4:498035/165‑1 (MQ=255)
                                                                                                                                                                              aaCTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGttcttt                                                                                                                                                                                                                <  5:190307/149‑1 (MQ=255)
                                                                                                                                                                               aCTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGttctttc                                                                                                                                                                                                               <  6:465438/149‑1 (MQ=255)
                                                                                                                                                                               aCTATTTTCTTGTTGCAGGTAGTCCAGGACGTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGc                                                                                 <  2:214450/275‑1 (MQ=255)
                                                                                                                                                                               aCTATGTTCTTGTTGACGGGATGCCATGTCTTCTGGTTGAAAACCCTGTTCGACCCGCGGCCAGTAGTAATTAGCGATGCGCGTGTTGCCATGAAGGGTCATGGTGTGGGGAATGCGGGCTGCTTTTTGTATCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTTTGCAGACAGAGAAAATCAGc                                                                                 <  4:101721/275‑1 (MQ=255)
                                                                                                                                                                                 tATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTg                                                                                                                                                                                                             >  5:131335/1‑149 (MQ=255)
                                                                                                                                                                                  aTTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTgg                                                                                                                                                                                                            <  7:278626/149‑1 (MQ=255)
                                                                                                                                                                                    tttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGtt                                                                                                                                                                                                          >  8:213922/1‑149 (MQ=255)
                                                                                                                                                                                            tgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATc                                                                                                                                                                                                  <  6:409256/149‑1 (MQ=255)
                                                                                                                                                                                               aGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTta                                                                                                                                                                                               >  8:100851/1‑149 (MQ=255)
                                                                                                                                                                                                ggTAGTCCGTGATGTCTGGGTTGGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATAGATGGGCGGATCGCCATGAAGCGTCATGGCGTGGTGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACAGCTTATAAGGGTGGCAAAACGCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCtgtg                                                                <  4:52340/275‑1 (MQ=255)
                                                                                                                                                                                                      ccAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGt                                                                                                                                                                                        <  5:249253/149‑1 (MQ=255)
                                                                                                                                                                                                       cAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTCTTGTTCTTTCTGGTTGAAACATCTTATAAGGGTg                                                                                                                                                                                       >  7:389040/1‑149 (MQ=255)
                                                                                                                                                                                                        aGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTgg                                                                                                                                                                                      >  6:21351/1‑149 (MQ=255)
                                                                                                                                                                                                             ttCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACAt                                                                                                                                                                                                   <  8:140520/131‑1 (MQ=255)
                                                                                                                                                                                                             ttCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACAt                                                                                                                                                                                                   >  7:140520/1‑131 (MQ=255)
                                                                                                                                                                                                             ttCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCaaaa                                                                                                                                                                                 <  5:162866/149‑1 (MQ=255)
                                                                                                                                                                                                               cTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTCTAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATTGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTggtgatcc                                                 >  2:137687/1‑269 (MQ=255)
                                                                                                                                                                                                                tGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCaaaa                                                                                                                                                                                 >  8:27299/1‑146 (MQ=255)
                                                                                                                                                                                                                tGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCaaaa                                                                                                                                                                                 <  7:27299/146‑1 (MQ=255)
                                                                                                                                                                                                                tGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTc                                                                                                                                                                              >  6:464145/1‑149 (MQ=255)
                                                                                                                                                                                                                    tGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTg                                                                                                                                                                                                                                    >  5:309752/1‑91 (MQ=255)
                                                                                                                                                                                                                    tGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTg                                                                                                                                                                                                                                    <  6:309752/91‑1 (MQ=255)
                                                                                                                                                                                                                       gaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCgggg                                                                                                                                                                       >  6:436092/1‑149 (MQ=255)
                                                                                                                                                                                                                                cagTCCCGCAGCCAGTAGTAATTATCGATGCGCGTAACGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTc                                 <  2:482385/272‑1 (MQ=255)
                                                                                                                                                                                                                                   tCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCtt                                                                                                                                                                                                <  8:448201/112‑1 (MQ=255)
                                                                                                                                                                                                                                   tCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCtt                                                                                                                                                                                                >  7:448201/1‑112 (MQ=255)
                                                                                                                                                                                                                                   tCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGaaa                                                                                                                                                           <  6:372842/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGa                                                                                                                                                                                                        >  3:38738/1‑100 (MQ=255)
                                                                                                                                                                                                                                       gCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAgg                                                                                                                                                       <  7:213922/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAgg                                                                                                                                                       <  7:100851/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTTAACTTCAGGCTATc                         >  2:424916/1‑275 (MQ=255)
                                                                                                                                                                                                                                           ccAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACa                                                                                                                                                                                                    <  1:38738/100‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaTTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCt                                                                                                                                           <  5:216785/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCt                                                                                                                                                                                                 <  6:111863/93‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCt                                                                                                                                                                                                 >  5:111863/1‑93 (MQ=255)
                                                                                                                                                                                                                                                       aTCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCAt                                                                                                                                               >  7:151436/1‑141 (MQ=255)
                                                                                                                                                                                                                                                       aTCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCAt                                                                                                                                               <  8:151436/141‑1 (MQ=255)
                                                                                                                                                                                                                                                       aTCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAga                                                                                                                                       <  8:249968/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               tccccacGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACa                                                                                                                               >  5:111318/7‑149 (MQ=255)
                                                                                                                                                                                                                                                               gcgTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGCTGGCAAAACTCACTGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGACATCGACTGTTATGTAGACACAGAAAATCAGCCTGTTCAAGTGGCTGTGCGCTTCTGGATAGCCCGAACGAGTAAACTTCACGCTACCCACAGAGCGGACATATTCCCCCa  >  4:431412/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               gcgTATCGCCATGAAGCGTCATGGCGTGGGGAACGCGGGCGGCTTTTGGTAGCATGTTATTGttctttct                                                                                                                                                                                                              <  7:443779/70‑1 (MQ=255)
                                                                                                                                                                                                                                                               gcgTATCGCCATGAAGCGTCATGGCGTGGGGAACGCGGGCGGCTTTTGGTAGCATGTTATTGttctttct                                                                                                                                                                                                              >  8:443779/1‑70 (MQ=255)
                                                                                                                                                                                                                                                                  tATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAggg                                                                                                                                                      >  2:541301/1‑123 (MQ=255)
                                                                                                                                                                                                                                                                   cacGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCaa                                                                                                                                                                 <  7:262256/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cacGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCaa                                                                                                                                                                 >  8:262256/3‑111 (MQ=255)
                                                                                                                                                                                                                                                                      gCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGt                                                                                                                                                  <  4:541296/123‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATg                                                                                                                        <  6:103059/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ccATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATgg                                                                                                                       >  8:65454/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         aTGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGaa                                                                                                                                                                                                       <  5:426449/67‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aTGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGaa                                                                                                                                                                                                       >  6:426449/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                   
ATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCA  >  NC_000913/1928530‑1929058

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: