Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,928,802 G→T 7.2% T13K (ACG→AAG)  ptrB ← oligopeptidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,928,8020GT 7.2% 274.0 / 7.5 83T13K (ACG→AAG) ptrBoligopeptidase B
Reads supporting (aligned to +/- strand):  ref base G (43/34);  new base T (4/2);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

AAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATC  >  NC_000913/1928655‑1928951
                                                                                                                                                   |                                                                                                                                                     
aaGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTc                                                                                                                                                     >  1:1957480/1‑150 (MQ=255)
 aGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCa                                                                                                                                                    <  2:1087432/150‑1 (MQ=255)
    tGCGATCCTGCAAGGCTTGTTTTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATgg                                                                                                                                                 <  2:1537482/150‑1 (MQ=255)
        aTCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTg                                                                                                                                             <  1:476914/150‑1 (MQ=255)
        aTCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTg                                                                                                                                             <  1:309675/150‑1 (MQ=255)
           cTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTgggg                                                                                                                                          <  1:227023/150‑1 (MQ=255)
            tGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCgatacgc                                                                                                                                         <  1:741128/150‑7 (MQ=255)
            tGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGa                                                                                                                                         >  2:1747370/1‑150 (MQ=255)
            tGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGa                                                                                                                                         >  1:922944/1‑150 (MQ=255)
                  gCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcgg                                                                                                                                   >  2:1702576/1‑150 (MQ=255)
                      gttgtGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcgggcgg                                                                                                                               >  1:937229/1‑150 (MQ=255)
                       ttgtGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGggcggc                                                                                                                              <  2:1957480/150‑1 (MQ=255)
                       ttgtGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGggcggc                                                                                                                              >  2:784985/1‑150 (MQ=255)
                         gtgAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtt                                                                                                                            >  2:1787167/1‑150 (MQ=255)
                           gAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtttt                                                                                                                          <  1:59092/150‑1 (MQ=255)
                               ccATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTa                                                                                                                      >  2:814005/1‑150 (MQ=255)
                                cATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATTGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAg                                                                                                                     >  2:325371/1‑150 (MQ=255)
                                 aTCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGc                                                                                                                    <  2:878567/150‑1 (MQ=255)
                                    aCCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGTCGGCTTTTGGTATCATg                                                                                                                 <  2:344479/150‑1 (MQ=255)
                                    aCCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGt                                                                                                                       >  2:579559/1‑144 (MQ=255)
                                    aCCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGt                                                                                                                       <  1:579559/144‑1 (MQ=255)
                                     cccGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGt                                                                                                                <  1:1735333/150‑1 (MQ=255)
                                       cGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTa                                                                                                              <  2:1529894/150‑1 (MQ=255)
                                             ccGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGttc                                                                                                        >  2:420638/1‑150 (MQ=255)
                                              cGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGttct                                                                                                       >  1:1460283/1‑150 (MQ=255)
                                               gTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGttctt                                                                                                      >  1:1645429/1‑150 (MQ=255)
                                                   cTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTg                                                                                                  >  2:1468115/1‑150 (MQ=255)
                                                    tATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTgg                                                                                                 <  1:1230907/150‑1 (MQ=255)
                                                      ttttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGtt                                                                                               >  2:1774511/1‑150 (MQ=255)
                                                        ttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTgt                                                                                             <  2:1237264/150‑2 (MQ=255)
                                                        ttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGa                                                                                             >  2:943545/1‑150 (MQ=255)
                                                          gttgttgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCTCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGaaa                                                                                           <  1:1676514/149‑1 (MQ=255)
                                                             gttgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACAt                                                                                        >  1:1240055/1‑150 (MQ=255)
                                                               tgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGTTTGAAACATCt                                                                                      <  2:87610/150‑1 (MQ=255)
                                                               tgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCt                                                                                      >  1:387189/1‑150 (MQ=255)
                                                               tgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCt                                                                                      <  1:1787167/150‑1 (MQ=255)
                                                               tgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCt                                                                                      >  2:406659/1‑150 (MQ=255)
                                                                gCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCtt                                                                                     <  2:1645429/150‑1 (MQ=255)
                                                                   ggTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTtata                                                                                  >  1:1342012/1‑150 (MQ=255)
                                                                   ggTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTtata                                                                                  <  1:784985/150‑1 (MQ=255)
                                                                   ggTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTtata                                                                                  >  1:791847/1‑150 (MQ=255)
                                                                       gTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAggg                                                                              >  1:1272406/1‑150 (MQ=255)
                                                                       gTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCATATAAggg                                                                              >  1:1227670/1‑150 (MQ=255)
                                                                          cAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTgg                                                                           >  1:1669891/1‑150 (MQ=255)
                                                                           aGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGc                                                                          <  2:1272406/150‑1 (MQ=255)
                                                                               cTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATAGAAGCGCGTATCGCCATGAAGCGTCATGGCCTGGGGAATGCGGGCGGCTTTTGGTAGCATCTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCaaaa                                                                      >  2:312324/1‑150 (MQ=255)
                                                                                ttCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAAc                                                                     <  1:1665292/150‑1 (MQ=255)
                                                                                 tCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACt                                                                    >  1:1868644/1‑150 (MQ=255)
                                                                                  cTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGAAGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTc                                                                   >  1:673213/1‑150 (MQ=255)
                                                                                        gagaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCggg                                                             >  2:667470/1‑150 (MQ=255)
                                                                                         agaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCgggg                                                            <  2:937229/150‑1 (MQ=255)
                                                                                         agaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCgggg                                                            >  2:1112924/1‑150 (MQ=255)
                                                                                          gaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGa                                                           >  1:1595650/1‑150 (MQ=255)
                                                                                                gTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGGGGCAAAACTCACCGGGGATGCAAg                                                     <  2:1532354/150‑1 (MQ=255)
                                                                                                 tATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGc                                                    >  1:1104285/1‑150 (MQ=255)
                                                                                                 tATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGc                                                    <  1:1393985/150‑1 (MQ=255)
                                                                                                  aTCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCg                                                   <  1:1949840/150‑1 (MQ=255)
                                                                                                   tCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGa                                                  <  1:1277955/150‑1 (MQ=255)
                                                                                                     gTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGaaa                                                >  1:22385/1‑150 (MQ=255)
                                                                                                      tCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAAc                                               <  2:1868644/150‑1 (MQ=255)
                                                                                                           cAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAgggg                                          <  2:387189/150‑1 (MQ=255)
                                                                                                             gCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAg                                        >  2:17338/1‑150 (MQ=255)
                                                                                                               cAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTc                                      >  1:1089161/1‑150 (MQ=255)
                                                                                                               cAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTc                                      >  1:1246077/1‑150 (MQ=255)
                                                                                                               cAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGGTATTTCTGGGTGAAGAATCTTATTAGGGTCGCAAAAGTCACGGGGGATGGGAGGGAAACAGGGGTGTc                                      >  2:867515/1‑150 (MQ=255)
                                                                                                                agtagtAATTATCGAAGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCa                                     <  2:673213/150‑1 (MQ=255)
                                                                                                                      aaTTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCtt                               >  1:362438/1‑150 (MQ=255)
                                                                                                                       aTTCTCGATGCGCTTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTa                              <  1:867515/150‑1 (MQ=255)
                                                                                                                       aTTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTa                              <  1:814005/150‑1 (MQ=255)
                                                                                                                         tATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAga                            <  1:469867/150‑1 (MQ=255)
                                                                                                                          aTCGATGCGCGTATCGCAATGAAGCGACATGGCGTGGGGAATGCGGGCGGCATTTGGTAGCATGTTATTGTTCTATCAGGTTGAAACATCTTATAAGGGTGGCAAAACTCAACGCAGATGCAAGCGAAACAGGCGAGTCATTGATTacaa                           >  2:244898/1‑147 (MQ=255)
                                                                                                                              aTGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAGTGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGatgatg                       <  2:922944/150‑1 (MQ=255)
                                                                                                                              aTGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGatgatg                       <  1:1747370/150‑1 (MQ=255)
                                                                                                                                     tATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTa                >  2:333196/1‑150 (MQ=255)
                                                                                                                                     ccacGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTa                >  2:1744459/4‑150 (MQ=255)
                                                                                                                                     ccacGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTa                >  1:442986/4‑150 (MQ=255)
                                                                                                                                       tCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAAt              >  1:1970035/1‑150 (MQ=255)
                                                                                                                                       acGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAAt              >  1:604798/2‑150 (MQ=255)
                                                                                                                                          ccATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGgc           >  1:474137/1‑150 (MQ=255)
                                                                                                                                            aTGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGgcgc         >  2:1843937/1‑150 (MQ=255)
                                                                                                                                             tGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGcgcg        <  1:943545/150‑1 (MQ=255)
                                                                                                                                                aGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGat     >  1:1710567/1‑150 (MQ=255)
                                                                                                                                                 gCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGata    >  1:1428189/1‑150 (MQ=255)
                                                                                                                                                   gTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATc  <  1:17338/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATC  >  NC_000913/1928655‑1928951

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: