Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 93,800 T→A 5.2% L212H (CTT→CAT)  murE → UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091393,8000TA 5.2% 705.4 / 18.8 191L212H (CTT→CAT) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base T (87/94);  new base A (5/5);  total (92/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 202 total aligned reads displayed.

GTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATG  >  NC_000913/93652‑93948
                                                                                                                                                    |                                                                                                                                                    
gTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCt                                                                                                                                                      >  1:1112655/1‑149 (MQ=255)
 tCGATTTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATAACCtt                                                                                                                                                     <  4:414649/149‑1 (MQ=255)
  cGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTg                                                                                                                                                    <  3:869131/149‑1 (MQ=255)
   gATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGa                                                                                                                                                   <  4:81158/149‑1 (MQ=255)
   gATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGa                                                                                                                                                   <  3:46510/149‑1 (MQ=255)
    aTGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGAt                                                                                                                                                  <  2:1260342/149‑1 (MQ=255)
      gTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTa                                                                                                                                                >  4:583022/1‑149 (MQ=255)
        tCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCAGCTTGATTATc                                                                                                                                              <  4:235067/149‑1 (MQ=255)
        tCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATc                                                                                                                                              <  4:1678307/149‑1 (MQ=255)
         cAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATc                                                                                                                                              <  1:242950/148‑1 (MQ=255)
         cAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATc                                                                                                                                              >  2:242950/1‑148 (MQ=255)
          aGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCAt                                                                                                                                            <  2:315345/149‑1 (MQ=255)
            cATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATgg                                                                                                                                          <  4:1141488/149‑1 (MQ=255)
             aTGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGt                                                                                                                                         <  1:357542/149‑1 (MQ=255)
             aTGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGt                                                                                                                                         <  3:1447549/149‑1 (MQ=255)
                aGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGat                                                                                                                                      >  3:643905/1‑149 (MQ=255)
                  cTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGatat                                                                                                                                    >  2:722324/1‑149 (MQ=255)
                           cTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTa                                                                                                                           <  1:722324/149‑1 (MQ=255)
                           cTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTa                                                                                                                           <  1:1530341/149‑1 (MQ=255)
                            tggtggATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTAc                                                                                                                          <  2:545881/149‑1 (MQ=255)
                              gtggATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGa                                                                                                                        <  3:1848771/149‑1 (MQ=255)
                                ggATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGacg                                                                                                                      <  4:1338038/149‑3 (MQ=255)
                                 gATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCt                                                                                                                                                      >  1:344548/1‑116 (MQ=255)
                                 gATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAaca                                                                                                                              >  4:339920/1‑140 (MQ=255)
                                  aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATc                                                                                                                                              <  2:1734672/123‑1 (MQ=255)
                                  aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATc                                                                                                                                              >  1:1734672/1‑123 (MQ=255)
                                  aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGa                                                                                                                                       <  3:328783/130‑1 (MQ=255)
                                  aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGa                                                                                                                                       >  4:328783/1‑130 (MQ=255)
                                  aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTAc                                                                                                                          >  3:715068/1‑143 (MQ=255)
                                  aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTAc                                                                                                                          >  3:611877/1‑143 (MQ=255)
                                     aGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGaa                                                                                                                       >  2:230536/1‑143 (MQ=255)
                                     aGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGaa                                                                                                                       <  1:230536/143‑1 (MQ=255)
                                     aGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGaa                                                                                                                       >  4:666006/1‑143 (MQ=255)
                                     aGGGCGCGACGTTTTGCCCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCcgcg                                                                                                                 >  2:1809931/1‑149 (MQ=255)
                                        gcgcGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGaaa                                                                                                              <  1:521359/149‑1 (MQ=255)
                                            gACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACt                                                                                                                            >  1:834449/1‑131 (MQ=255)
                                            gACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgc                                                                                                          >  3:1350203/1‑149 (MQ=255)
                                             aCGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATgg                                                                                                                                  <  3:1130650/124‑1 (MQ=255)
                                              cGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgctg                                                                                                        <  2:1112655/149‑1 (MQ=255)
                                                ttttGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTAt                                                                                                                                               <  4:6292/108‑1 (MQ=255)
                                                ttttGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTAt                                                                                                                                               >  1:572825/1‑108 (MQ=255)
                                                ttttGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgctgct                                                                                                      <  4:1316577/149‑1 (MQ=255)
                                                    gcgcAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgctg                                                                                                        >  1:561015/1‑143 (MQ=255)
                                                     cgcAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTAtt                                                                                                 <  1:1041385/149‑1 (MQ=255)
                                                         aTGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGa                                                                                             <  1:681388/149‑1 (MQ=255)
                                                          tGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAg                                                                                                                      >  3:11124/1‑123 (MQ=255)
                                                          tGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAg                                                                                            <  4:1350203/149‑1 (MQ=255)
                                                            gAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGa                                                                                             >  4:1623337/1‑146 (MQ=255)
                                                            gAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACCAAGCCGCGAAATGGCTGCTTTATTCTGa                                                                                             <  3:1623337/146‑1 (MQ=255)
                                                                  tcctccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTg                                                                                    >  2:1396846/1‑149 (MQ=255)
                                                                   cctccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAGCACTACGAAGCCGCGAAAGGGCTGCTTTATTCTGAGCATCATTGc                                                                                   <  3:790758/149‑1 (MQ=255)
                                                                    ctccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAg                                                                                            <  2:57085/139‑1 (MQ=255)
                                                                    ctccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCg                                                                                  >  2:1718701/1‑149 (MQ=255)
                                                                     tccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCgg                                                                                 >  3:50977/1‑149 (MQ=255)
                                                                       ccACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgctg                                                                                                        <  2:683870/124‑1 (MQ=255)
                                                                        cACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCa                                                                              <  1:1809931/149‑1 (MQ=255)
                                                                        cACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCa                                                                              <  1:1385296/149‑1 (MQ=255)
                                                                        cACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCa                                                                              <  2:1341830/149‑1 (MQ=255)
                                                                         aCGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAg                                                                             <  4:101412/149‑1 (MQ=255)
                                                                              cTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                        <  2:654008/149‑1 (MQ=255)
                                                                              cTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                        >  1:94852/1‑149 (MQ=255)
                                                                                 gTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGattat                                                                     >  1:71656/1‑149 (MQ=255)
                                                                                 gTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGattat                                                                     <  4:1848659/149‑1 (MQ=255)
                                                                                          cGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGatt                                                                       <  4:513533/138‑1 (MQ=255)
                                                                                                  ggcATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAg                                                    <  4:731591/149‑1 (MQ=255)
                                                                                                   gcATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCAtt                                                                                     >  3:1416626/1‑115 (MQ=255)
                                                                                                   gcATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCAtt                                                                                     <  4:1416626/115‑1 (MQ=255)
                                                                                                      ttGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGaa                                                                                                                       >  1:1455172/1‑78 (MQ=255)
                                                                                                      ttGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGa                                                                                             <  1:1273724/104‑1 (MQ=255)
                                                                                                      ttGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGatt                                                                       <  1:1885773/126‑1 (MQ=255)
                                                                                                      ttGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGatt                                                                       >  2:1885773/1‑126 (MQ=255)
                                                                                                         aaaTTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                        >  2:1559709/1‑122 (MQ=255)
                                                                                                         aaaTTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                        <  1:1559709/122‑1 (MQ=255)
                                                                                                                 ggcgTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCgctggc                                     <  2:909560/149‑1 (MQ=255)
                                                                                                                  gcgTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                        >  1:291205/1‑113 (MQ=255)
                                                                                                                        gTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGt                                                                                                                                         <  3:1012214/42‑1 (MQ=255)
                                                                                                                        gTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTc                                                                                                >  2:472826/1‑83 (MQ=255)
                                                                                                                        gTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTc                                                                                                <  1:472826/83‑1 (MQ=255)
                                                                                                                        gTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCaaa                              >  4:48894/1‑149 (MQ=255)
                                                                                                                         tCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGctggctgg                                  <  2:1054645/144‑1 (MQ=255)
                                                                                                                           tttACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACt                           >  2:1486650/1‑149 (MQ=255)
                                                                                                                           tttACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACt                           >  4:1566472/1‑149 (MQ=255)
                                                                                                                                cAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCgg                      >  4:1793156/1‑149 (MQ=255)
                                                                                                                                cAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCgg                      >  3:1553162/1‑149 (MQ=255)
                                                                                                                                 aaCTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGa                     >  1:1053607/1‑149 (MQ=255)
                                                                                                                                  aCTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGt                                                   >  2:356175/1‑118 (MQ=255)
                                                                                                                                        gCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCt                                                                                               <  2:999272/68‑1 (MQ=255)
                                                                                                                                          ccataTCACCATGATAATCAAGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGcc                        <  2:746082/133‑1 (MQ=255)
                                                                                                                                             gATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgctgc                                                                                                       <  3:1316709/55‑1 (MQ=255)
                                                                                                                                                cACCATGATAATCAAGGTGATATGGAACACTACGAAGCCGCGAAATGgctgc                                                                                                       <  3:994033/52‑1 (MQ=25)
                                                                                                                                                cACCATGATAATCAAGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATc      >  2:1465525/1‑149 (MQ=255)
                                                                                                                                                  ccTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGa                                                                                             >  4:1124240/1‑60 (MQ=255)
                                                                                                                                                  ccTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCaa    <  1:1718701/149‑1 (MQ=255)
                                                                                                                                                   cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGattatt                                                                    >  3:558982/1‑84 (MQ=255)
                                                                                                                                                   cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGAcgc                  <  2:1186047/134‑1 (MQ=255)
                                                                                                                                                   cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGAcgc                  >  1:1186047/1‑134 (MQ=255)
                                                                                                                                                   cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAAt   >  1:1349005/1‑149 (MQ=255)
                                                                                                                                                   cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAAt   <  1:1065045/149‑1 (MQ=255)
                                                                                                                                                   cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGATTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGGAAAACTGCCGGACGCGGTTGCGGTATCAAt   >  1:1219166/1‑149 (MQ=255)
                                                                                                                                                    ttGATTATCATGGTGATATGGAACACTACGAAGTCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATg  >  3:1197435/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATG  >  NC_000913/93652‑93948

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