Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,970,392 C→T P87L (CCA→CTA)  wecB → UDP‑N‑acetyl glucosamine‑2‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,3920CT99.0% 354.2 / ‑4.2 95P87L (CCA→CTA) wecBUDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (58/36);  total (59/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

AGCATCGGGAGATGCTCGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTAC  >  NC_000913/3970245‑3970539
                                                                                                                                                   |                                                                                                                                                   
agcaTCGGGAGATGCTCGATCAGGTGCTGTAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGa                                                                                                                                                  <  1:967098/151‑1 (MQ=255)
  caTCGGGAGATGCTCGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGAcg                                                                                                                                                <  2:87709/151‑1 (MQ=255)
     cGGGAGATGCTCGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGAcgtc                                                                                                                                              <  1:1344486/150‑1 (MQ=255)
       ggAGATGCTCGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTg                                                                                                                                           <  1:829268/151‑1 (MQ=255)
          gaTGCTCGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTg                                                                                                                                        <  2:455966/151‑1 (MQ=255)
             gCTCGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGt                                                                                                                                      <  1:708788/150‑1 (MQ=255)
             gCTCGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACCACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGtt                                                                                                                                     >  2:4321/1‑151 (MQ=255)
                cGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCAc                                                                                                                                  <  1:604771/151‑1 (MQ=255)
                cGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCAc                                                                                                                                  >  1:300818/1‑151 (MQ=255)
                 gATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACg                                                                                                                                 >  1:1292278/1‑151 (MQ=255)
                  aTCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGGGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACgg                                                                                                                                >  2:1403572/1‑151 (MQ=255)
                  aTCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACgg                                                                                                                                >  1:492123/1‑151 (MQ=255)
                  aTCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACgg                                                                                                                                >  2:553257/1‑151 (MQ=255)
                     aGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCg                                                                                                                              >  2:1255817/1‑150 (MQ=255)
                         gCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATacg                                                                                                                         >  1:887869/1‑151 (MQ=255)
                         gCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATacg                                                                                                                         >  2:722814/1‑151 (MQ=255)
                          cTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATacg                                                                                                                         >  2:1117042/1‑150 (MQ=255)
                            gAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATacgacg                                                                                                                      <  1:1104445/151‑1 (MQ=255)
                             aaaCTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTTACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATacgacg                                                                                                                      >  1:748785/1‑150 (MQ=255)
                             aaaCTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATacgacg                                                                                                                      >  1:1115512/1‑150 (MQ=255)
                                  ctTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCATGGCGATacgacgacgacg                                                                                                                >  2:652820/1‑151 (MQ=255)
                                   ttttCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATacgacgacgacg                                                                                                                >  1:691197/1‑150 (MQ=255)
                                      tCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTgg                                                                                                            >  1:1335143/1‑151 (MQ=255)
                                      tCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTg                                                                                                             <  1:805217/150‑1 (MQ=255)
                                        cATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGc                                                                                                           >  2:241387/1‑150 (MQ=255)
                                         aTTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCaa                                                                                                         <  2:342992/151‑1 (MQ=255)
                                            gTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCa                                                                                                      <  1:1280419/151‑1 (MQ=255)
                                               ccTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGcc                                                                                                   >  2:943897/1‑151 (MQ=255)
                                               ccTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGcc                                                                                                   <  1:934395/151‑1 (MQ=255)
                                                   aCTACGATCTCAACCAAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGACGATACGACGACGACGCTGGCAACCAGCCTggc                                                                                               <  1:223186/151‑1 (MQ=255)
                                                    cTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTggcg                                                                                              >  1:663081/1‑151 (MQ=255)
                                                      aCGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTggcggc                                                                                            >  1:693946/1‑151 (MQ=255)
                                                       cGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTggcggc                                                                                            >  1:1122699/1‑150 (MQ=255)
                                                        gATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGt                                                                                          >  2:834257/1‑151 (MQ=255)
                                                        gATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGt                                                                                          >  1:1316460/1‑151 (MQ=255)
                                                          tctcAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGttt                                                                                        >  1:495070/1‑151 (MQ=255)
                                                          tctcAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGttt                                                                                        >  2:941613/1‑151 (MQ=255)
                                                             cAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCACCCAGCCTGGCGGCGTTTTAt                                                                                     <  2:462754/151‑1 (MQ=255)
                                                             cAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGtttt                                                                                       <  1:389420/149‑1 (MQ=255)
                                                             cAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGt                                                                                          <  2:311594/146‑1 (MQ=255)
                                                             cAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTAt                                                                                     >  1:427862/1‑151 (MQ=255)
                                                              caCATAATGCAGCCAGGACAGGGTCTGACAGAGATAACATGTCGGGTTCTGGAAGTGCTAAAACGTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGAAACGACGACGACGCTCGCAACCAGCCTGGCGGCGTTTTAt                                                                                     >  1:1167670/2‑150 (MQ=255)
                                                                cATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAg                                                                                  <  1:444924/151‑1 (MQ=255)
                                                                 ttAATGCAACCAGGACAGGGCCTGACAGAGATAACGTGCATGTTTCTGGAAGTACTAAAATCTATTCTTGCCGAGTTCAAACTAGACGGCGTGGTGGTTCACGGCGATACGACGACGACGCTGGCAACTAGCCTGGCGGCGGTATATCAGc                                                                                 >  1:154153/2‑151 (MQ=255)
                                                                 aTAATGCAGCCAGGACAGGGGCTGACAAAGATAACCTGTCGGATTCTGGCAGGGCTAAAACCTAGTCTTGCCGAGTGCAAACTAGACGTCGTGCTGGTTGACGGCGAGACGACGACGACGCTGGCAACCAGCCTTGCGGCGTTTTATAAg                                                                                  >  1:615080/1‑150 (MQ=255)
                                                                   aaTGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGt                                                                               >  2:586612/1‑151 (MQ=255)
                                                                    aTGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGt                                                                               >  2:1180679/1‑150 (MQ=255)
                                                                     tGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTCTCAGCGTAt                                                                             >  1:587393/1‑151 (MQ=255)
                                                                       cAGCAAGGACAGGGCCTGAAAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGCCAACAAGCCTGGCGGCGTTTTATCAGCGTATTc                                                                           >  1:107833/1‑151 (MQ=255)
                                                                        aGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTcc                                                                          >  2:1172621/1‑151 (MQ=255)
                                                                         gCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCt                                                                         <  2:902809/151‑1 (MQ=255)
                                                                         gCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCt                                                                         <  1:1412654/151‑1 (MQ=255)
                                                                            aGGACAGTGCCAGACAGAGATAACCAGGCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGt                                                                       >  1:712079/1‑150 (MQ=255)
                                                                                 aGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCa                                                                 <  2:494431/151‑1 (MQ=255)
                                                                                 aGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCa                                                                 <  2:450578/151‑1 (MQ=255)
                                                                                 aGGGCATGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTGTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGGTACGTCGACGACGTCGGCAACCAGTCAGGCGGCGGTTTATCAGGGTATTCCCGGTgg                                                                    >  1:1277792/1‑148 (MQ=255)
                                                                                  gggCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCAc                                                                <  2:878919/151‑1 (MQ=255)
                                                                                  gggCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCAc                                                                >  1:643274/1‑151 (MQ=255)
                                                                                     ccTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGtt                                                             <  2:422174/151‑1 (MQ=255)
                                                                                        gACAGAGATACCCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGTCGATACGACGACGGCGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAg                                                          <  2:41130/151‑1 (MQ=255)
                                                                                          cAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGc                                                        >  1:880241/1‑151 (MQ=255)
                                                                                           agagATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGc                                                        >  1:226369/1‑150 (MQ=255)
                                                                                           agagATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGc                                                        >  1:1416844/1‑150 (MQ=255)
                                                                                            gagaTAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTg                                                      >  2:1306549/1‑151 (MQ=255)
                                                                                               aTAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTc                                                   >  1:187716/1‑151 (MQ=255)
                                                                                                tAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTc                                                   >  1:1384015/1‑150 (MQ=255)
                                                                                                tAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCt                                                  >  1:115223/1‑151 (MQ=255)
                                                                                                tAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCt                                                  >  1:1387627/1‑151 (MQ=255)
                                                                                                 aaCCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCt                                                  <  1:775233/150‑1 (MQ=255)
                                                                                                      gTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTgcgc                                              <  1:962967/149‑1 (MQ=255)
                                                                                                      gTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCAc                                            <  1:684681/151‑1 (MQ=255)
                                                                                                        cGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACg                                           <  2:576096/150‑1 (MQ=255)
                                                                                                                  aaGGGCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATctct                                 >  1:903793/1‑150 (MQ=255)
                                                                                                                    gggCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTAtt                              >  2:3262/1‑151 (MQ=255)
                                                                                                                     ggCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTc                             <  2:550703/151‑1 (MQ=255)
                                                                                                                     ggCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTc                             >  2:752838/1‑151 (MQ=255)
                                                                                                                      gCTAAAACCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCg                            <  2:148219/151‑1 (MQ=255)
                                                                                                                      gCTAAAACCTATTCTAGCCGAGTTCAAACTAGACGTCGAGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCg                            >  1:1375385/1‑151 (MQ=255)
                                                                                                                         aaaaCCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCg                         >  2:29284/1‑151 (MQ=255)
                                                                                                                          aaaCCTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGt                        >  2:327128/1‑151 (MQ=255)
                                                                                                                             ccTATTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGc                     >  2:1227008/1‑151 (MQ=255)
                                                                                                                                aTTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCgg                  >  1:836052/1‑151 (MQ=255)
                                                                                                                                aTTCTTGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCgg                  >  2:793971/1‑151 (MQ=255)
                                                                                                                                     tGCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAag               >  2:356080/1‑149 (MQ=255)
                                                                                                                                      gCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAagag             <  1:139589/150‑1 (MQ=255)
                                                                                                                                      gCCGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAgg            >  2:1007868/1‑151 (MQ=255)
                                                                                                                                       ccGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGc           <  1:948493/151‑1 (MQ=255)
                                                                                                                                        cGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCt          <  1:1045549/151‑1 (MQ=255)
                                                                                                                                        cGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCt          >  2:109218/1‑151 (MQ=255)
                                                                                                                                        cGAGTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCt          <  2:691229/151‑1 (MQ=255)
                                                                                                                                           gTTCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAAc       >  1:856071/1‑151 (MQ=255)
                                                                                                                                             tCAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCg     >  2:82864/1‑151 (MQ=255)
                                                                                                                                              cAAACTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGt    >  2:755553/1‑151 (MQ=255)
                                                                                                                                               aaaCTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTa   <  1:264330/151‑1 (MQ=255)
                                                                                                                                                aaCTAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTAc  <  1:525435/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCATCGGGAGATGCTCGATCAGGTGCTGAAACTCTTTTCCATTGTACCTGACTACGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTAC  >  NC_000913/3970245‑3970539

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