Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,161,241 C→A 9.2% intergenic (‑95/+16) ygiQ ← / ← ftsP radical SAM superfamily protein YgiQ/cell division protein required during stress conditions

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,161,2410CA 9.2% 303.4 / 14.5 98intergenic (‑95/+16)ygiQ/ftsPradical SAM superfamily protein YgiQ/cell division protein required during stress conditions
Reads supporting (aligned to +/- strand):  ref base C (48/40);  new base A (3/7);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTC  >  NC_000913/3161095‑3161387
                                                                                                                                                  |                                                                                                                                                  
gCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCggg                                                                                                                                                 >  2:983089/1‑150 (MQ=255)
  ggCCAGGAGAACAGGGCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCg                                                                                                                                               <  1:1040268/150‑1 (MQ=255)
    ccAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAg                                                                                                                                             <  1:1432207/150‑1 (MQ=255)
    ccAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAg                                                                                                                                             <  1:119943/150‑1 (MQ=255)
    ccAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAg                                                                                                                                             <  1:495263/150‑1 (MQ=255)
     cAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGt                                                                                                                                            <  1:1159207/150‑1 (MQ=255)
      aGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTa                                                                                                                                           >  1:1510117/1‑150 (MQ=255)
      aGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTa                                                                                                                                           <  1:1398746/150‑1 (MQ=255)
       ggAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTTCATACTCGCCCGAGGAGGGCGAGTAt                                                                                                                                          <  2:1116702/150‑1 (MQ=255)
       ggAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCGAGGAGGGCGAGTAt                                                                                                                                          <  2:840120/150‑1 (MQ=255)
       ggAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCGAGGAGGGCGAGTAt                                                                                                                                          <  1:1192791/150‑1 (MQ=255)
        gagaACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATg                                                                                                                                         <  2:1188903/150‑1 (MQ=255)
        gagaACAGGGCGCGATCCGGTTGGATCAGGGATATAGAGCTCATTTTGCTTCCACAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCGAGGAGGGCGAGTATg                                                                                                                                         <  2:464904/150‑1 (MQ=255)
         agaACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTTCATACTCGCCCGAGGAGGGCGAGTATgg                                                                                                                                        <  2:930501/150‑2 (MQ=255)
          gaACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGaa                                                                                                                                       <  2:574704/150‑1 (MQ=255)
          gaACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCGAGGAGGGCGAGTATgga                                                                                                                                       <  1:921131/150‑3 (MQ=255)
           aaCAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAg                                                                                                                                      <  2:222519/150‑1 (MQ=255)
             cAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCGAGGAGGGCGAGTATggaata                                                                                                                                    <  1:1004770/150‑6 (MQ=255)
               ggTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTa                                                                                                                                  >  2:490912/1‑150 (MQ=255)
                gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTAc                                                                                                                                 >  1:148625/1‑150 (MQ=255)
                 tCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACg                                                                                                                                <  2:880544/150‑1 (MQ=255)
                  cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACgg                                                                                                                               <  1:1065680/150‑1 (MQ=255)
                         cGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGa                                                                                                                        <  1:983089/150‑1 (MQ=255)
                          ggTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGAt                                                                                                                       >  2:263324/1‑150 (MQ=255)
                          ggTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGAt                                                                                                                       <  2:362362/150‑1 (MQ=255)
                            ttGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTg                                                                                                                     <  2:1015850/150‑1 (MQ=255)
                             tGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGa                                                                                                                    >  1:574961/1‑150 (MQ=255)
                               gATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGAcc                                                                                                                  >  1:343936/1‑150 (MQ=255)
                                aTCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCa                                                                                                                 >  1:1242365/1‑150 (MQ=255)
                                  cAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAAc                                                                                                               >  2:1010241/1‑150 (MQ=255)
                                      gagaTAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGtt                                                                                                           >  1:209242/1‑150 (MQ=255)
                                         aTAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGcc                                                                                                        <  2:392169/150‑1 (MQ=255)
                                         aTAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGcc                                                                                                        <  2:178451/150‑1 (MQ=255)
                                         aTAGAGCTCATTTTGCTTCCAGAAATGAAAAAAAAATAATCAAAGGCCGCGGATTATAAGCCGGAACGAAAGACAAATCGAAAGGTATTACATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGcc                                                                                                        >  2:1111659/1‑150 (MQ=255)
                                          tAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGccc                                                                                                       >  2:699166/1‑150 (MQ=255)
                                          tAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGccc                                                                                                       <  2:1510117/150‑1 (MQ=255)
                                              gCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATc                                                                                                   >  2:1019703/1‑150 (MQ=255)
                                                 cATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAg                                                                                                <  2:129387/150‑1 (MQ=255)
                                                 cATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAg                                                                                                <  2:201242/150‑1 (MQ=255)
                                                  aTTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACAAACAGTTGCCCAATCGAGc                                                                                               >  1:377658/1‑150 (MQ=255)
                                                     ttGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCAc                                                                                            >  2:808065/1‑150 (MQ=255)
                                                       gCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAAGCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACgg                                                                                          <  2:910592/150‑1 (MQ=255)
                                                          tCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTcc                                                                                       >  2:1258511/1‑150 (MQ=255)
                                                           ccAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCg                                                                                      >  1:1378884/1‑150 (MQ=255)
                                                           ccAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCg                                                                                      >  1:1357845/1‑150 (MQ=255)
                                                           ccAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCACAATCGAGCCACGGTccc                                                                                      >  2:1295359/1‑149 (MQ=255)
                                                            cAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGTCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGc                                                                                     <  2:625178/150‑1 (MQ=255)
                                                            cAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTAGCACTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAAACGAGCCACGGTcccc                                                                                     >  1:1132455/1‑148 (MQ=255)
                                                             aGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGcc                                                                                    <  2:202192/150‑1 (MQ=255)
                                                             aGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGcc                                                                                    >  1:1246930/1‑150 (MQ=255)
                                                             aGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGcc                                                                                    <  2:1378884/150‑1 (MQ=255)
                                                              gAATTGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCa                                                                                   <  1:1362057/150‑1 (MQ=255)
                                                              gAAATGATAAAAAAATAATCAAAGGCCGGGGATTATTAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCa                                                                                   <  1:1362055/150‑1 (MQ=255)
                                                              gAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCa                                                                                   >  2:292119/1‑150 (MQ=255)
                                                                aaTGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCAtt                                                                                 >  1:367530/1‑150 (MQ=255)
                                                                   gATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGGGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTcc                                                                              >  1:1115951/1‑150 (MQ=255)
                                                                   gATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTcc                                                                              >  1:919661/1‑150 (MQ=255)
                                                                     tAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAg                                                                            >  1:1376972/1‑150 (MQ=255)
                                                                       aaaaaaTAATCAAAGCCCCAGGATTATAATCCGGAACCAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTAGGCATATTAAGAAGATTAAGGTACCGCATTTACAAACAGTTGCACAATCTAGCCACGGTACGCCATTTCCAGCg                                                                          >  2:76563/1‑150 (MQ=255)
                                                                       aaaaaaTAATAAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTAAGGTACCGGATTGACCAACAGTTGCCAAATCGAGCCACGGTCCGCCATTTCCAGCg                                                                          >  2:1130900/1‑150 (MQ=255)
                                                                         aaaaTAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTc                                                                        <  1:1019703/150‑1 (MQ=255)
                                                                           aataatCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTg                                                                      <  2:148625/150‑1 (MQ=255)
                                                                               atCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTg                                                                  >  1:741498/1‑150 (MQ=255)
                                                                                  aaaGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTg                                                               <  2:963471/150‑1 (MQ=255)
                                                                                  aaaGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTg                                                               <  2:1353452/150‑1 (MQ=255)
                                                                                  aaaGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTAAGTTAACGTATTGACCAACAGTTGCCCAATCTAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTg                                                               <  2:968519/150‑1 (MQ=255)
                                                                                  aaaGGCCGGGGATTATAAGCCGGAACGAAAGAGAAAGCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACTGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGGCCGCCATTTCCCGCTTCTGACTGTTg                                                               <  2:192489/150‑1 (MQ=255)
                                                                                     ggCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAg                                                            >  1:820223/1‑150 (MQ=255)
                                                                                      gCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGt                                                           <  1:16002/150‑1 (MQ=255)
                                                                                        cGGGGATTATAAGCCGGAACGAAAGAGAAATCGCAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGCTTTCGGTACCGGATGGACCGACAGTTGCCCGATCGAGCCACGGTCCGCCATTTCCAGCGTGTGACTGTTGAAGTAg                                                         >  2:1255769/1‑150 (MQ=255)
                                                                                          gggATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGaa                                                       >  1:381126/1‑150 (MQ=255)
                                                                                                ataAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGaa                                                 <  1:76563/150‑1 (MQ=255)
                                                                                                ataAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGaa                                                 <  1:541939/150‑1 (MQ=255)
                                                                                                 taAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAg                                                >  1:920097/1‑150 (MQ=255)
                                                                                                 taAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAg                                                >  1:1138815/1‑150 (MQ=255)
                                                                                                 taAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAg                                                >  1:545623/1‑150 (MQ=255)
                                                                                                      cGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTgcgc                                           <  2:390286/150‑1 (MQ=255)
                                                                                                       ggAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAACTAGAACGGGAAGTGCGcc                                          >  2:272260/1‑150 (MQ=255)
                                                                                                         aaCGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCa                                        <  1:1258511/150‑1 (MQ=255)
                                                                                                           cGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCCTTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAgg                                      >  2:988340/1‑150 (MQ=255)
                                                                                                            gAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGa                                     >  2:864508/1‑150 (MQ=255)
                                                                                                             aaaGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGaa                                    >  2:806665/1‑150 (MQ=255)
                                                                                                                gagaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGc                                 >  2:925452/1‑150 (MQ=255)
                                                                                                                 agaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCt                                >  1:1021586/1‑150 (MQ=255)
                                                                                                                 agaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCt                                <  1:292119/150‑1 (MQ=255)
                                                                                                                  gaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTg                               >  1:1113734/1‑150 (MQ=255)
                                                                                                                   aaaTCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGATTATGAAGATTACGGTACCGGATTGACCAAAAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGa                              >  1:1286034/1‑150 (MQ=255)
                                                                                                                        gAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCgca                         >  1:1520590/1‑148 (MQ=255)
                                                                                                                        gAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGaa                         >  1:599752/1‑150 (MQ=255)
                                                                                                                         aaaGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGaaa                        <  2:1232496/150‑1 (MQ=255)
                                                                                                                               aTTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACa                  >  1:353482/1‑150 (MQ=255)
                                                                                                                                ttCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACaa                 >  2:1077640/1‑150 (MQ=255)
                                                                                                                                  tcATACTCGCCCGAGGAGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAagc               >  1:1397454/2‑150 (MQ=255)
                                                                                                                                   cATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAagca              <  1:490912/150‑1 (MQ=255)
                                                                                                                                         cGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCa        >  1:1143812/1‑150 (MQ=255)
                                                                                                                                         cGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCa        >  1:1419282/1‑150 (MQ=255)
                                                                                                                                           cccTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACt      <  2:1113734/150‑1 (MQ=255)
                                                                                                                                            cctcctAGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCAAATCGAGAAACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACtt     >  2:707015/1‑150 (MQ=255)
                                                                                                                                             ctcctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTg    <  2:1242365/150‑1 (MQ=255)
                                                                                                                                               cctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTc  >  2:249997/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTC  >  NC_000913/3161095‑3161387

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