Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,721,507 A→C 6.9% I82L (ATC→CTC)  ydhJ → putative membrane fusion protein (MFP) of YdhJK efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,721,5070AC 6.9% 304.8 / 13.9 87I82L (ATC→CTC) ydhJputative membrane fusion protein (MFP) of YdhJK efflux pump
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base C (3/3);  total (43/44)
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.63e-01

GGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAATTTTATTTCTGCCGAAGAACTCGATACCGCTAACCTC  >  NC_000913/1721364‑1721641
                                                                                                                                               |                                                                                                                                      
ggACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGc                                                                                                                                    <  2:981810/148‑1 (MQ=255)
     cgcgACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGCTCGGCGatc                                                                                                                               >  2:128776/1‑146 (MQ=255)
          cGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGCTCGGCGATCTGaaa                                                                                                                          <  1:128776/148‑2 (MQ=255)
            ggAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGc                                                                                                                        >  2:215583/1‑148 (MQ=255)
               aaaTACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGcaca                                                                                                                     <  1:215583/148‑1 (MQ=255)
                    cgcgCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTc                                                                                                                <  1:657889/148‑1 (MQ=255)
                                  ttCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGCTCGGCGATCTGaaag                                                                                                                         >  2:415835/1‑123 (MQ=255)
                                  ttCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGCTCGGCGATCTGaaag                                                                                                                         <  1:415835/125‑3 (MQ=255)
                                       tCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTg                                                                                                                             >  1:129041/1‑116 (MQ=255)
                                       tCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTg                                                                                                                             <  2:129041/116‑1 (MQ=255)
                                               caGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGCTCGGCGATCTGaa                                                                                                                           <  1:286572/110‑1 (MQ=255)
                                               caGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGCTCGGCGATCTGaa                                                                                                                           >  2:286572/1‑110 (MQ=255)
                                                   tgtCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAgcg                                                                                                  <  4:457527/131‑1 (MQ=255)
                                                   tgtCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAgcg                                                                                                  >  3:457527/1‑131 (MQ=255)
                                                        ggACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTgg                                                                                                          <  1:857540/118‑1 (MQ=255)
                                                        ggACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTgg                                                                                                          >  2:857540/1‑118 (MQ=255)
                                                          aCGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGcac                                                                                                                      <  4:544133/104‑1 (MQ=255)
                                                          aCGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGcac                                                                                                                      >  3:544133/1‑104 (MQ=255)
                                                              aTTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGa                                                                                                                                        >  1:407390/1‑82 (MQ=255)
                                                              aTTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGa                                                                                                                                        <  2:407390/82‑1 (MQ=255)
                                                                             ataAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAg                                                                                                                                >  4:907639/1‑75 (MQ=255)
                                                                             ataAAAGACAACCAGCTGGTTAACGCGGGGGATATTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAg                                                                                                                                <  3:907639/75‑1 (MQ=255)
                                                                               aaaaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgtcg                                                         >  3:155369/1‑144 (MQ=255)
                                                                               aaaaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgtcg                                                         <  4:155369/144‑1 (MQ=255)
                                                                               aaaaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgt                                                           <  1:1011471/142‑1 (MQ=255)
                                                                               aaaaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgt                                                           >  2:1011471/1‑142 (MQ=255)
                                                                                aaaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTgg                                                                                                          >  1:833772/1‑94 (MQ=255)
                                                                                aaaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTgg                                                                                                          <  2:833772/94‑1 (MQ=255)
                                                                                aaaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCa                                                                            <  3:459380/124‑1 (MQ=255)
                                                                                aaaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCa                                                                            >  4:459380/1‑124 (MQ=255)
                                                                                 aaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCt                                                                                                                 <  2:12647/86‑1 (MQ=255)
                                                                                 aaGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCt                                                                                                                 >  1:12647/1‑86 (MQ=255)
                                                                                  aGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCg                                                                                                                                  >  3:310676/1‑68 (MQ=255)
                                                                                  aGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCg                                                                                                                                  >  1:677571/1‑68 (MQ=255)
                                                                                  aGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCg                                                                                                                                  <  2:677571/68‑1 (MQ=255)
                                                                                  aGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCg                                                                                                                                  <  4:310676/68‑1 (MQ=255)
                                                                                  aGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGca                                                                                                                       <  2:865997/79‑1 (MQ=255)
                                                                                  aGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGca                                                                                                                       >  1:865997/1‑79 (MQ=255)
                                                                                     cAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAATGCTCAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCaaa                                                                                 >  4:737481/1‑114 (MQ=255)
                                                                                     cAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCa                                                                                                <  1:726611/99‑1 (MQ=255)
                                                                                     cAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCa                                                                                                >  2:726611/1‑99 (MQ=255)
                                                                                     cAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCaaa                                                                                 <  3:737481/114‑1 (MQ=255)
                                                                                         cAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGaa                                                                                                                           <  4:431596/68‑1 (MQ=255)
                                                                                         cAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGaa                                                                                                                           >  3:431596/1‑68 (MQ=255)
                                                                                              ggTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCAtc                                                 >  4:482956/1‑137 (MQ=255)
                                                                                              ggTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCAtc                                                 <  3:482956/137‑1 (MQ=255)
                                                                                                ttAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCa                                                                            >  1:24920/1‑108 (MQ=255)
                                                                                                ttAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCa                                                                            <  2:24920/108‑1 (MQ=255)
                                                                                                ttAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGa                                                                     >  4:622003/1‑115 (MQ=255)
                                                                                                ttAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGa                                                                     <  3:622003/115‑1 (MQ=255)
                                                                                                 tAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCg                                                                                                                                  <  4:684443/53‑1 (MQ=255)
                                                                                                 tAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCg                                                                                                                                  >  3:684443/1‑53 (MQ=255)
                                                                                                 tAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCa                                                                                                <  2:43714/87‑1 (MQ=255)
                                                                                                 tAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCa                                                                                                >  1:43714/1‑87 (MQ=255)
                                                                                                 tAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCaaa                                                                                 <  1:1043924/102‑1 (MQ=255)
                                                                                                 tAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCaaa                                                                                 >  2:1043924/1‑102 (MQ=255)
                                                                                                 tAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAg                                                                               >  4:545672/1‑104 (MQ=255)
                                                                                                 tAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAg                                                                               <  3:545672/104‑1 (MQ=255)
                                                                                                      cgGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGaa                                                                                                      <  1:1115149/76‑1 (MQ=255)
                                                                                                      cgGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGaa                                                                                                      >  2:1115149/1‑76 (MQ=255)
                                                                                                         gggATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGa                                                                >  3:476984/1‑111 (MQ=255)
                                                                                                         gggATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAAACTGACCTCGCAAAAGCCAATAACGAAGCGa                                                                <  4:476984/111‑1 (MQ=255)
                                                                                                             tcttctGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgt                                                           >  2:61583/1‑112 (MQ=255)
                                                                                                             tcttctGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgt                                                           <  1:61583/112‑1 (MQ=255)
                                                                                                                tctGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAAttt                                    >  1:1105807/1‑132 (MQ=255)
                                                                                                                tctGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAAttt                                    <  2:1105807/132‑1 (MQ=255)
                                                                                                                            cgatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCaa                                                                           <  4:451317/81‑1 (MQ=255)
                                                                                                                            cgatAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCaa                                                                           >  3:451317/1‑81 (MQ=255)
                                                                                                                              atAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAATTTTATTTCTGCCgaaga                     <  1:978993/133‑1 (MQ=255)
                                                                                                                              atAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAATTTTATTTCTGCCgaaga                     >  2:978993/1‑133 (MQ=255)
                                                                                                                                aaaaCGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATcgtcgtcgt                                                     >  4:894142/1‑99 (MQ=255)
                                                                                                                                aaaaCGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATcgtcgtcgt                                                     <  3:894142/99‑1 (MQ=255)
                                                                                                                                 aaaCGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAg                                                                                                    <  4:264404/51‑1 (MQ=255)
                                                                                                                                 aaaCGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAg                                                                                                    >  3:264404/1‑51 (MQ=255)
                                                                                                                                  aaCGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAATTTTATTTCTGCCGAAGAACTCGATACCGCTAACCTc  >  3:536131/1‑148 (MQ=255)
                                                                                                                                    cGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgt                                                           >  1:567511/1‑89 (MQ=255)
                                                                                                                                    cGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgt                                                           <  2:567511/89‑1 (MQ=255)
                                                                                                                                      cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgt                                                           >  2:99521/1‑87 (MQ=255)
                                                                                                                                      cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAAtcgt                                                           <  1:99521/87‑1 (MQ=255)
                                                                                                                                      cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAATTTTAttt                               >  2:154875/1‑115 (MQ=255)
                                                                                                                                      cccTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAATTTTAttt                               <  1:154875/115‑1 (MQ=255)
                                                                                                                                        cTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCa                                                                            >  1:1093585/1‑68 (MQ=255)
                                                                                                                                        cTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCa                                                                            <  2:1093585/68‑1 (MQ=255)
                                                                                                                                          ttCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCAtctctc                                             <  3:490026/97‑1 (MQ=255)
                                                                                                                                          ttCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCAtctctc                                             >  4:490026/1‑97 (MQ=255)
                                                                                                                                          ttCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCAtctct                                              >  2:994469/1‑96 (MQ=255)
                                                                                                                                          ttCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCAtctct                                              <  1:378626/96‑1 (MQ=255)
                                                                                                                                          ttCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCAtctct                                              <  1:994469/96‑1 (MQ=255)
                                                                                                                                          ttCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCAtctct                                              >  2:378626/1‑96 (MQ=255)
                                                                                                                                               aTCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAATTTTATTTCTGCCGAAGAACTCGATACCGCTAAcc    >  2:1083548/1‑133 (MQ=255)
                                                                                                                                               aTCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAATTTTATTTCTGCCGAAGAACTCGATACCGCTAAcc    <  1:1083548/133‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                      
GGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGGATCTTCTGCTCACGATCGATAAAACGCCCTTTCAGATCGCCGAGCTGAACGCACAGGCTCAGCTGGCGAAAGCGCAATCTGACCTCGCAAAAGCCAATAACGAAGCGAATCGTCGTCGTCATCTCTCACAAAATTTTATTTCTGCCGAAGAACTCGATACCGCTAACCTC  >  NC_000913/1721364‑1721641

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