Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 93,810 T→A 5.3% H215Q (CAT→CAA)  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,8100TA 5.3% 691.7 / 21.3 189H215Q (CAT→CAA) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base T (88/91);  new base A (5/5);  total (93/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 191 total aligned reads displayed.

AGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATA  >  NC_000913/93662‑93958
                                                                                                                                                    |                                                                                                                                                    
aGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCAt                                                                                                                                                      <  2:315345/149‑1 (MQ=255)
 gCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATg                                                                                                                                                     >  1:545881/1‑149 (MQ=255)
  cATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATgg                                                                                                                                                    <  3:583022/149‑1 (MQ=255)
  cATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATgg                                                                                                                                                    <  4:193004/149‑1 (MQ=255)
      aGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGat                                                                                                                                                >  3:643905/1‑149 (MQ=255)
      aGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCATGATAATCAAGGTGat                                                                                                                                                <  4:1759861/149‑1 (MQ=255)
        cTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGatat                                                                                                                                              >  2:722324/1‑149 (MQ=255)
         tGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATg                                                                                                                                             >  2:189874/1‑149 (MQ=255)
          ggCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATgg                                                                                                                                            <  1:1109189/149‑1 (MQ=255)
                 cTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTa                                                                                                                                     <  1:722324/149‑1 (MQ=255)
                  tggtggATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTAc                                                                                                                                    <  2:545881/149‑1 (MQ=255)
                    gtggATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGa                                                                                                                                  <  3:1848771/149‑1 (MQ=255)
                      ggATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAg                                                                                                                                <  3:1529617/149‑1 (MQ=255)
                        aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGa                                                                                                                                                 >  4:328783/1‑130 (MQ=255)
                        aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGa                                                                                                                                                 <  3:328783/130‑1 (MQ=255)
                        aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTAc                                                                                                                                    <  4:611877/143‑1 (MQ=255)
                        aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTAc                                                                                                                                    >  3:715068/1‑143 (MQ=255)
                        aTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGcc                                                                                                                              >  3:1316577/1‑149 (MQ=255)
                           aGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGaa                                                                                                                                 <  1:230536/143‑1 (MQ=255)
                           aGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGaa                                                                                                                                 >  4:666006/1‑143 (MQ=255)
                           aGGGCGCGACGTTTTGCCCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCcgcg                                                                                                                           >  2:1809931/1‑149 (MQ=255)
                                  gACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACt                                                                                                                                      >  1:834449/1‑131 (MQ=255)
                                  gACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACt                                                                                                                                      <  2:834449/131‑1 (MQ=255)
                                  gACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgc                                                                                                                    >  3:1350203/1‑149 (MQ=255)
                                   aCGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATgg                                                                                                                                            >  4:1130650/1‑124 (MQ=255)
                                          gcgcAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgctg                                                                                                                  >  1:561015/1‑143 (MQ=255)
                                           cgcAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTAtt                                                                                                           <  1:1041385/149‑1 (MQ=255)
                                               aTGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGa                                                                                                       <  1:681388/149‑1 (MQ=255)
                                                tGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAg                                                                                                                                <  4:11124/123‑1 (MQ=255)
                                                tGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAg                                                                                                      >  3:101412/1‑149 (MQ=255)
                                                tGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAg                                                                                                      <  4:1350203/149‑1 (MQ=255)
                                                 ggAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgct                                                                                                                   >  3:792310/1‑135 (MQ=255)
                                                      tttCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATg                                                                                                                                             >  2:1200192/1‑104 (MQ=255)
                                                        tcctccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAaca                                                                                                                                        >  2:231772/1‑107 (MQ=255)
                                                        tcctccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTg                                                                                              >  2:1396846/1‑149 (MQ=255)
                                                         cctccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAGCACTACGAAGCCGCGAAAGGGCTGCTTTATTCTGAGCATCATTGc                                                                                             <  3:790758/149‑1 (MQ=255)
                                                          ctccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAg                                                                                                      <  2:57085/139‑1 (MQ=255)
                                                          ctccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCg                                                                                            <  4:347435/149‑1 (MQ=255)
                                                           tccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCAtt                                                                                               <  3:1161538/145‑1 (MQ=255)
                                                           tccCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCAtt                                                                                               >  4:1161538/1‑145 (MQ=255)
                                                             ccACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTc                                                                                         >  1:1341830/1‑149 (MQ=255)
                                                              cACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCa                                                                                        <  1:1809931/149‑1 (MQ=255)
                                                                    cTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                                  <  2:654008/149‑1 (MQ=255)
                                                                    cTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                                  <  1:104154/149‑1 (MQ=255)
                                                                    cTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                                  >  1:94852/1‑149 (MQ=255)
                                                                       gTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGa                                                                                                                                           >  2:1722901/1‑89 (MQ=255)
                                                                       gTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGa                                                                                                                                           <  1:1722901/89‑1 (MQ=255)
                                                                       gTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGattat                                                                               >  1:71656/1‑149 (MQ=255)
                                                                          cAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTaa                                                                            <  4:183609/149‑1 (MQ=255)
                                                                              aCCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGcc                                                                        >  3:1008850/1‑149 (MQ=255)
                                                                                cGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGatt                                                                                 >  3:513533/1‑138 (MQ=255)
                                                                                cGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGatt                                                                                 <  4:513533/138‑1 (MQ=255)
                                                                                       cggcATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACTCCGACGATGaa                                                               <  1:1396846/149‑1 (MQ=255)
                                                                                        ggcATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAg                                                              <  4:731591/149‑1 (MQ=255)
                                                                                         gcATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCAtt                                                                                               <  4:1416626/115‑1 (MQ=255)
                                                                                         gcATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCAtt                                                                                               >  3:1416626/1‑115 (MQ=255)
                                                                                          cATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTg                                                            >  1:276663/1‑149 (MQ=255)
                                                                                           aTTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTgg                                                           >  1:1357433/1‑149 (MQ=255)
                                                                                           aTTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTgg                                                           >  4:1657154/1‑149 (MQ=255)
                                                                                            ttGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGa                                                                                                       >  2:1273724/1‑104 (MQ=255)
                                                                                            ttGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGatt                                                                                 <  1:1885773/126‑1 (MQ=255)
                                                                                            ttGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGatt                                                                                 >  2:1885773/1‑126 (MQ=255)
                                                                                               aaaTTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                                  <  1:1559709/122‑1 (MQ=255)
                                                                                               aaaTTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                                  >  2:1559709/1‑122 (MQ=255)
                                                                                                              gTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGt                                                                                                                                                   >  4:1012214/1‑42 (MQ=255)
                                                                                                              gTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTc                                                                                                          >  2:472826/1‑83 (MQ=255)
                                                                                                              gTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCaaa                                        >  4:48894/1‑149 (MQ=255)
                                                                                                               tCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGctggctgg                                            >  1:1054645/1‑144 (MQ=255)
                                                                                                                 tttACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACt                                     >  2:1486650/1‑149 (MQ=255)
                                                                                                                 tttACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACt                                     >  4:1566472/1‑149 (MQ=255)
                                                                                                                    aCCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGcc                                  >  2:1065045/1‑149 (MQ=255)
                                                                                                                    aCCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGcc                                  <  2:94852/149‑1 (MQ=255)
                                                                                                                     ccAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGat                                                                                  <  1:1481180/100‑1 (MQ=255)
                                                                                                                      cAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCgg                                >  3:1553162/1‑149 (MQ=255)
                                                                                                                       aaCTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGa                               >  2:676062/1‑149 (MQ=255)
                                                                                                                        aCTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGATGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGAc                              >  1:419788/1‑149 (MQ=255)
                                                                                                                        aCTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGt                                                             <  1:356175/118‑1 (MQ=255)
                                                                                                                               tccataTCACCATGATAATCAAGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCgctgg                                                >  4:90863/6‑124 (MQ=255)
                                                                                                                                  cgATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCgg                    <  2:419788/149‑1 (MQ=255)
                                                                                                                                   gATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGgctgc                                                                                                                 >  4:1316709/1‑55 (MQ=255)
                                                                                                                                    aTCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTa                  >  2:728296/1‑149 (MQ=255)
                                                                                                                                      cACCATGATAATCAAGGTGATATGGAACACTACGAAGCCGCGAAATGgctgc                                                                                                                 <  3:994033/52‑1 (MQ=25)
                                                                                                                                      cACCATGATAATCAAGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATc                >  2:1465525/1‑149 (MQ=255)
                                                                                                                                       aCCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCa               >  1:1536598/1‑149 (MQ=255)
                                                                                                                                        ccTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGa                                                                                                       <  3:1124240/60‑1 (MQ=255)
                                                                                                                                         cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGAcgc                            <  2:1186047/134‑1 (MQ=255)
                                                                                                                                         cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGAcgc                            >  1:1186047/1‑134 (MQ=255)
                                                                                                                                         cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAAt             >  1:485354/1‑149 (MQ=255)
                                                                                                                                         cTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGATTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGGAAAACTGCCGGACGCGGTTGCGGTATCAAt             >  1:1219166/1‑149 (MQ=255)
                                                                                                                                           tGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATgg           <  1:626125/149‑1 (MQ=255)
                                                                                                                                            gATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGc                                   <  1:1229700/124‑1 (MQ=255)
                                                                                                                                            gATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGa          >  3:426478/1‑149 (MQ=255)
                                                                                                                                              ttATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAAt             >  1:1710353/1‑144 (MQ=255)
                                                                                                                                              ttATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAAt             <  2:1710353/144‑1 (MQ=255)
                                                                                                                                               tATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGa       <  3:48894/149‑1 (MQ=255)
                                                                                                                                               tATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGa       <  2:1349005/149‑1 (MQ=255)
                                                                                                                                                 tCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATc     <  2:1053607/149‑1 (MQ=255)
                                                                                                                                                   aTGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCAtt                                                                                               <  1:67546/57‑1 (MQ=255)
                                                                                                                                                   aTGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCat   >  2:706154/1‑149 (MQ=255)
                                                                                                                                                    tGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCata  <  3:1793156/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGATATGGAACACTACGAAGCCGCGAAATGGCTGCTTTATTCTGAGCATCATTGCGGTCAGGCGATTATTAACGCCGACGATGAAGTGGGCCGCCGCTGGCTGGCAAAACTGCCGGACGCGGTTGCGGTATCAATGGAAGATCATA  >  NC_000913/93662‑93958

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: ‑