Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,709,561 (C)6→7 coding (915/1692 nt) eptB ← KDO phosphoethanolamine transferase, Ca(2+)‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,709,5551.C97.9% 340.3 / ‑1.0 97coding (921/1692 nt)eptBKDO phosphoethanolamine transferase, Ca(2+)‑inducible
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (49/46);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-01

Only 100 of 142 total aligned reads displayed.

AGCAATGTTGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAATACCCAT  >  NC_000913/3709417‑3709689
                                                                                                                                             |                                                                                                                                        
 gCAATGTTGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑gcg‑‑‑c                                                                                                                                       >  2:486170/1‑139 (MQ=255)
      gttgtCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cc                                                                                                                                      >  2:1398770/1‑135 (MQ=255)
      gttgtCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cc                                                                                                                                      <  1:1398770/135‑1 (MQ=255)
       ttgtCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑ccccccc                                                                                                                                 <  2:663900/139‑2 (MQ=255)
       ttgtCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑ccccccc                                                                                                                                 <  2:430548/139‑2 (MQ=255)
         gtCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑ccccccctg                                                                                                                               <  1:1517190/139‑4 (MQ=255)
          tCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑ccccccctga                                                                                                                              >  1:1283213/1‑135 (MQ=255)
          tCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctg                                                                                                                              >  1:770868/1‑135 (MQ=255)
          tCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctg                                                                                                                              >  2:1369302/1‑135 (MQ=255)
          tCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctg                                                                                                                              >  2:1509527/1‑135 (MQ=255)
           ccGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCc‑‑gcg                                                                                                                                           <  2:417445/128‑1 (MQ=255)
            cGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctc                                                                                                                              >  2:1098672/1‑133 (MQ=255)
            cGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑ccccccc                                                                                                                                 >  1:895618/1‑133 (MQ=255)
            cGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑ccccccc                                                                                                                                 <  2:895618/134‑2 (MQ=255)
             gCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctgacg                                                                                                                           <  2:770868/139‑8 (MQ=255)
               cATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctgacgta                                                                                                                         >  1:1030766/1‑130 (MQ=255)
                         cTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctgacg                                                                                                                           >  1:880303/1‑120 (MQ=255)
                         cTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctgacg                                                                                                                           <  2:880303/127‑8 (MQ=255)
                           gTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAAc                                                                                                               <  1:72728/139‑1 (MQ=255)
                            tAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACg                                                                                                              <  2:467687/139‑1 (MQ=255)
                            tAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACg                                                                                                              >  1:16392/1‑139 (MQ=255)
                                 ccACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTga                                                                                                        >  2:434119/1‑139 (MQ=255)
                                         cGCTCTGCATAGCGTAGAGGTCAGAACTAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                               >  1:107005/1‑139 (MQ=255)
                                           ctctGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccct                                                                                                                               >  1:1650337/1‑102 (MQ=255)
                                             ctGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTg                                                                                             <  1:1369302/139‑1 (MQ=255)
                                              tGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                            <  1:1509527/139‑1 (MQ=255)
                                              tGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                            <  2:16392/139‑1 (MQ=255)
                                                cATAGCGTAGAGGTCAGAACTGAATCCTAACGGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                          <  1:286597/139‑1 (MQ=255)
                                                     cGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACaa                                                                                                                        >  1:1634778/1‑104 (MQ=255)
                                                     cGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                                     >  1:606299/1‑139 (MQ=255)
                                                     cGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATcac                                                                                    <  1:550037/139‑1 (MQ=255)
                                                        agagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAg                                                                                  >  2:1650468/1‑139 (MQ=255)
                                                        agagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAgg                                                                                 >  2:1038829/1‑139 (MQ=255)
                                                         gagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAgg                                                                                 <  1:592816/139‑1 (MQ=255)
                                                          agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGa                                                                                >  1:416876/1‑139 (MQ=255)
                                                          agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAg                                                                               >  2:367142/1‑139 (MQ=255)
                                                          agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccct                                                                                                                               <  1:1130492/90‑4 (MQ=255)
                                                          agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cc                                                                                                                                      <  2:856325/83‑1 (MQ=255)
                                                          agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cc                                                                                                                                      >  1:856325/1‑83 (MQ=255)
                                                           ggTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGt                                                                              <  2:1636821/139‑1 (MQ=255)
                                                            gTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGt                                                                              <  2:1484739/139‑1 (MQ=255)
                                                               aGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGagag                                                                                                     >  1:417821/1‑113 (MQ=255)
                                                               aGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccct                                                                                                                               >  2:1315541/1‑82 (MQ=255)
                                                                     gAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGaa                                                                    <  1:434119/139‑1 (MQ=255)
                                                                      aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑ccccccc                                                                                                                                 >  1:490201/1‑75 (MQ=255)
                                                                           tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                          <  1:1160924/112‑1 (MQ=255)
                                                                           tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                          >  2:1160924/1‑112 (MQ=255)
                                                                           tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccg                                                               <  2:103159/139‑1 (MQ=255)
                                                                           tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccg                                                               >  1:103159/1‑139 (MQ=255)
                                                                           tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccgc                                                              >  1:74319/1‑139 (MQ=255)
                                                                            aaCTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                               <  1:940342/105‑1 (MQ=255)
                                                                            aaCTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                               >  2:940342/1‑105 (MQ=255)
                                                                             aCTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                          >  1:980770/1‑109 (MQ=255)
                                                                             aCTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                          <  2:980770/109‑1 (MQ=255)
                                                                              cTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccc                                                                                                                                  >  2:163221/1‑67 (MQ=255)
                                                                              cTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccc                                                                                                                                  <  1:163221/67‑1 (MQ=255)
                                                                                gCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGAt                                                         >  2:68225/1‑139 (MQ=255)
                                                                                gCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctgacg                                                                                                                           >  2:649841/1‑65 (MQ=255)
                                                                                gCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑cccccccctgacg                                                                                                                           <  1:649841/72‑8 (MQ=255)
                                                                                  ttCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGAtt                                                        >  2:1251198/1‑139 (MQ=255)
                                                                                         ccGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATcac                                                                                    <  2:515457/104‑1 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccgc                                                              >  2:1363992/1‑123 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTgg                                              >  2:109311/1‑139 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTgg                                              >  2:55318/1‑139 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGt                                                                                                          <  1:1287127/78‑1 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGt                                                                                                          >  2:1287127/1‑78 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGaa                                                                    <  1:658755/116‑1 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGGACGCCGCCAGATTTTTCTCCTggg                                             >  2:1484267/1‑139 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTggg                                             >  2:1299574/1‑139 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTggg                                             >  1:1330/1‑139 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑ccccccc                                                                                                                                 <  1:604193/54‑2 (MQ=255)
                                                                                            cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑ccccccc                                                                                                                                 >  2:604193/1‑53 (MQ=255)
                                                                                                     tctgttctTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCg                                    >  2:636100/1‑139 (MQ=255)
                                                                                                        gttctTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCg                                 <  2:1065797/139‑1 (MQ=255)
                                                                                                         ttctTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACaa                                                                                                                        >  1:501245/1‑52 (MQ=255)
                                                                                                         ttctTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACaa                                                                                                                        <  2:501245/52‑1 (MQ=255)
                                                                                                         ttctTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCg                                 <  1:752247/139‑1 (MQ=255)
                                                                                                          tctTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTg                               <  1:1023354/139‑1 (MQ=255)
                                                                                                               tAATGTGCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGACCGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTAtt                          <  1:1080238/139‑1 (MQ=255)
                                                                                                                   gtgCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                               <  2:857489/67‑1 (MQ=255)
                                                                                                                   gtgCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                               >  1:857489/1‑67 (MQ=255)
                                                                                                                   gtgCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACg                       >  1:1560259/1‑139 (MQ=255)
                                                                                                                   gtgCGCTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACg                       >  2:1264492/1‑139 (MQ=255)
                                                                                                                   gtgCGCTGCGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGc                      >  1:280735/1‑139 (MQ=255)
                                                                                                                      cgcTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                                 <  2:1616319/110‑1 (MQ=255)
                                                                                                                      cgcTGCGGATTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTc                    >  1:132864/1‑139 (MQ=255)
                                                                                                                       gcTGCGGATTATCTTCCGTGCC‑‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCa                   <  1:1406336/139‑1 (MQ=255)
                                                                                                                           cGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTAtt                          <  1:843851/127‑1 (MQ=255)
                                                                                                                           cGGATTATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTAtt                          >  2:843851/1‑127 (MQ=255)
                                                                                                                              aTTATCTTCC‑‑GCGCCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                                 >  2:146682/1‑103 (MQ=255)
                                                                                                                              aTTATCTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGa            <  1:55318/139‑1 (MQ=255)
                                                                                                                                tATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCg                                 >  2:593719/1‑115 (MQ=255)
                                                                                                                                tATCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGaaaa         >  2:1452025/1‑139 (MQ=255)
                                                                                                                                  tCTTCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAATa       <  1:1484267/139‑1 (MQ=255)
                                                                                                                                   cTTCC‑‑GCGCCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAATa       <  1:1251198/139‑1 (MQ=255)
                                                                                                                                     tCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAATa       >  2:811719/1‑136 (MQ=255)
                                                                                                                                     tCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAATa       <  1:811719/136‑1 (MQ=255)
                                                                                                                                     tCC‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAATAccc    >  2:653253/1‑139 (MQ=255)
                                                                                                                                     cgc‑‑gcc‑‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAATACCCa   >  2:1284984/3‑139 (MQ=255)
                                                                                                                                      cc‑‑GCGC‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAATACCCa   <  1:147532/139‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                        
AGCAATGTTGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCC‑‑GCG‑‑‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAATACCCAT  >  NC_000913/3709417‑3709689

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