Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,970,817:1 +G 7.4% coding (685/1131 nt) wecB → UDP‑N‑acetylglucosamine 2‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,8171.G 7.4% 373.4 / 26.3 95coding (685/1131 nt)wecBUDP‑N‑acetylglucosamine 2‑epimerase
Reads supporting (aligned to +/- strand):  ref base . (42/46);  new base G (2/5);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

GGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGG  >  NC_000913/3970668‑3970966
                                                                                                                                                      |                                                                                                                                                     
gggTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACgcgc                                                                                                                                                        >  2:1255684/1‑150 (MQ=255)
  gTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTg                                                                                                                                                     <  2:694708/150‑1 (MQ=255)
   tGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTgg                                                                                                                                                    <  2:1033539/150‑1 (MQ=255)
   tGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCAGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTgg                                                                                                                                                    <  1:387944/150‑1 (MQ=255)
    gCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGc                                                                                                                                                   >  1:1166285/1‑150 (MQ=255)
    gCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGc                                                                                                                                                   <  2:172680/150‑1 (MQ=255)
     cGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCa                                                                                                                                                  >  1:1446765/1‑150 (MQ=255)
      gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTAGCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCTCAGGCGTGAGAGTTTCGTTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAg                                                                                                                                                 <  2:376175/150‑1 (MQ=255)
      gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAg                                                                                                                                                 >  1:58684/1‑150 (MQ=255)
         aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGatt                                                                                                                                              <  1:1023976/150‑3 (MQ=255)
         aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGatt                                                                                                                                              <  2:241027/150‑3 (MQ=255)
         aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGGTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTCGCTTTGAAGAAATCTGCCACGCGCTGGCAGACa                                                                                                                                              <  1:526347/150‑1 (MQ=255)
          ccAGGTGATGAGCCGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGattt                                                                                                                                             <  2:1503056/150‑4 (MQ=255)
           cAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATc                                                                                                                                            >  1:797878/1‑150 (MQ=255)
              gTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGcc                                                                                                                                         <  1:71759/150‑1 (MQ=255)
              gTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGatttcttc                                                                                                                                         <  1:901688/150‑8 (MQ=255)
              gTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGGGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGatttcttc                                                                                                                                         <  1:1556798/150‑8 (MQ=255)
               tgatgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGcca                                                                                                                                        <  2:166262/150‑1 (MQ=255)
                gatgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGCCCCCGATAAAAAGATGATTCTGGTGACCGGGCACAGGCGTGAGAGTTTCGGTCGGGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccac                                                                                                                                       <  1:352979/150‑1 (MQ=255)
                   gaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCAGGCAGACATCGCcaccac                                                                                                                                    >  1:486977/1‑150 (MQ=255)
                     gcagcGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGc                                                                                                                                  <  1:738654/150‑1 (MQ=255)
                     gcagcGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGGCACAGGCGCGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGc                                                                                                                                  <  2:1439949/150‑1 (MQ=255)
                       agcGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGACTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCAc                                                                                                                                <  2:1144220/150‑1 (MQ=255)
                           aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCAGCACGCACCAgg                                                                                                                            <  2:256932/150‑1 (MQ=255)
                           aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAgg                                                                                                                            <  1:190850/150‑1 (MQ=255)
                                       cAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                                <  2:1004091/150‑1 (MQ=255)
                                         gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGTTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTc                                                                                                              <  2:1379964/150‑1 (MQ=255)
                                               gcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTc                                                                                                              >  2:1572950/1‑144 (MQ=255)
                                               gcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTc                                                                                                              <  1:1572950/144‑1 (MQ=255)
                                               gcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCg                                                                                                        <  2:1150110/150‑1 (MQ=255)
                                                cggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCgg                                                                                                       <  2:432247/150‑1 (MQ=255)
                                                cggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCgg                                                                                                       >  1:1267908/1‑150 (MQ=255)
                                                 ggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGt                                                                                                      <  1:1369649/150‑1 (MQ=255)
                                                 ggcAAATTACCCGTTTATCGACCCCAATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGt                                                                                                      >  2:195964/1‑150 (MQ=255)
                                                   cAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGCCATCCAGATTGTCTATCCGGTGc                                                                                                    >  1:496152/1‑150 (MQ=255)
                                                    aaaTTACCCGTTTATCGACCCAGATAAAAAGATGATTCTGGTGACCTGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTCCa                                                                                                   <  1:719735/150‑1 (MQ=255)
                                                         aCCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGATTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCa                                                                                              <  2:1548418/150‑1 (MQ=255)
                                                          cccGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                             >  2:1099742/1‑150 (MQ=255)
                                                           ccGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                            <  2:58684/150‑1 (MQ=255)
                                                            cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACa                                                                                                                         <  1:840707/120‑1 (MQ=255)
                                                            cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACa                                                                                                                         >  2:840707/1‑120 (MQ=255)
                                                            cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAcc                                                                                           >  2:928187/1‑150 (MQ=255)
                                                            cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAcc                                                                                           <  1:675213/150‑1 (MQ=255)
                                                            cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAcc                                                                                           >  1:149528/1‑150 (MQ=255)
                                                                tATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGaa                                                                                       >  2:1094119/1‑150 (MQ=255)
                                                                    gACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTc                                                                                   <  2:24791/150‑1 (MQ=255)
                                                                     aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGCCATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCa                                                                                  >  2:1082832/1‑150 (MQ=255)
                                                                     aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCa                                                                                  <  1:1255684/150‑1 (MQ=255)
                                                                       cccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCaga                                                                                >  2:203179/1‑150 (MQ=255)
                                                                        ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTc                                                                                                              >  1:650797/1‑119 (MQ=255)
                                                                        ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTc                                                                                                              <  2:650797/119‑1 (MQ=255)
                                                                        ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCagag                                                                               >  1:210956/1‑150 (MQ=255)
                                                                        ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCagag                                                                               <  2:612976/150‑1 (MQ=255)
                                                                         cGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAgaga                                                                              >  1:131987/1‑150 (MQ=255)
                                                                            tAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAAcc                                                                           <  1:1436711/150‑1 (MQ=255)
                                                                             aaaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCg                                                                          >  1:1473593/1‑150 (MQ=255)
                                                                              aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCgg                                                                         >  2:190549/1‑150 (MQ=255)
                                                                              aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCgg                                                                         >  2:966184/1‑150 (MQ=255)
                                                                                  gatgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                                     >  1:516752/1‑150 (MQ=255)
                                                                                   atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAAt                                                                    <  1:1162205/150‑1 (MQ=255)
                                                                                   atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAAt                                                                    <  2:749576/150‑1 (MQ=255)
                                                                                       ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTGagag                                                                               >  1:1303005/1‑135 (MQ=255)
                                                                                       ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCagag                                                                               <  2:1303005/135‑1 (MQ=255)
                                                                                       ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCa                                                                >  2:1587581/1‑150 (MQ=255)
                                                                                           ggTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCt                                                            >  1:984923/1‑150 (MQ=255)
                                                                                              gACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTggg                                                         <  1:59537/150‑1 (MQ=255)
                                                                                                ccGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGc                                                       >  1:290311/1‑150 (MQ=255)
                                                                                                  ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAGCCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                                     >  1:1236470/1‑150 (MQ=255)
                                                                                                  ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGa                                                                                        <  2:962340/115‑1 (MQ=255)
                                                                                                  ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGa                                                                                        >  1:962340/1‑115 (MQ=255)
                                                                                                  ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACg                                                                                     <  1:1593783/118‑1 (MQ=255)
                                                                                                  ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACg                                                                                     >  2:1593783/1‑118 (MQ=255)
                                                                                                     cacaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgtg                                                  <  2:1166285/150‑1 (MQ=255)
                                                                                                       caGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaa                                                >  1:760468/1‑150 (MQ=255)
                                                                                                             tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGt                                          >  1:1392931/1‑150 (MQ=255)
                                                                                                               agagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCa                                        <  1:1587581/150‑1 (MQ=255)
                                                                                                                gagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACTCGCTTGCAGACATCGCCACGACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCTCATTCTGGGGCATGTGAAAAATGTCAt                                       <  2:1353056/150‑1 (MQ=255)
                                                                                                                     tCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTga                                  >  2:1203817/1‑150 (MQ=255)
                                                                                                                       ggTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatc                                >  1:1053757/1‑150 (MQ=255)
                                                                                                                       ggTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatc                                <  1:203179/150‑1 (MQ=255)
                                                                                                                        gTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcg                               >  1:626773/1‑150 (MQ=255)
                                                                                                                        gTCGTGGCGTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCATATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTTAAAAATGTCATTCTgatcg                               <  2:486977/150‑1 (MQ=255)
                                                                                                                           gTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcgatc                            >  1:1058696/1‑150 (MQ=255)
                                                                                                                             ggCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATccc                          >  2:1482743/1‑150 (MQ=255)
                                                                                                                             ggCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATccc                          <  2:210956/150‑1 (MQ=255)
                                                                                                                                tttGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAg                       <  1:289206/150‑1 (MQ=255)
                                                                                                                                 ttGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAgg                      <  2:619382/150‑1 (MQ=255)
                                                                                                                                  tGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                     >  2:1414090/1‑150 (MQ=255)
                                                                                                                                  tGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCAGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                     >  2:1441646/1‑150 (MQ=255)
                                                                                                                                   gaagaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg                    >  1:185261/1‑150 (MQ=255)
                                                                                                                                     cgatgtCTGCCAGCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa                  >  1:1064666/6‑150 (MQ=255)
                                                                                                                                     cgatgtCTGCCAGCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa                  >  1:585934/6‑150 (MQ=255)
                                                                                                                                         aTCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTa              <  1:1011001/150‑1 (MQ=255)
                                                                                                                                            tGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTAcag           <  1:1262726/150‑3 (MQ=255)
                                                                                                                                                aCGCGCTGGCAGCCATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGCACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTg       >  2:1495004/1‑150 (MQ=255)
                                                                                                                                                aCGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTg       >  2:980452/1‑150 (MQ=255)
                                                                                                                                                  gcgcTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCTTAATTTAATTTAACGGGTAAATCGGATTGTTTGTCTTGTGAAAAATGTTTTTCTTATCGATGCTCATTTGTTTTTACGGtttttg     <  2:149528/150‑4 (MQ=255)
                                                                                                                                                   cgcTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCt    <  2:131987/150‑1 (MQ=255)
                                                                                                                                                     cTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTgg  >  2:191382/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGG  >  NC_000913/3970668‑3970966

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