Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,586 A→C 6.9% intergenic (+46/‑35) rpmF → / → plsX 50S ribosomal subunit protein L32/putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,5860AC 6.9% 498.5 / 25.3 161intergenic (+46/‑35)rpmF/plsX50S ribosomal subunit protein L32/putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (68/82);  new base C (5/6);  total (73/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

Only 100 of 161 total aligned reads displayed.

CACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTT  >  NC_000913/1147441‑1147733
                                                                                                                                                 |                                                                                                                                                   
cACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTgg                                                                                                                                                 <  2:2597647/150‑1 (MQ=255)
   cAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGa                                                                                                                                              >  1:1143992/1‑150 (MQ=255)
   cAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGa                                                                                                                                              >  1:1586731/1‑150 (MQ=255)
   aaGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGa                                                                                                                                              >  1:1941554/2‑150 (MQ=255)
     gCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGTTGCCTGGGGaaa                                                                                                                                            <  1:45881/150‑1 (MQ=255)
     gCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGTTGCCTGGGGaaa                                                                                                                                            <  1:1782486/150‑1 (MQ=255)
     gCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGaaa                                                                                                                                            <  1:2795695/150‑1 (MQ=255)
     gCCGGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGTTGCCTGGGGaaa                                                                                                                                            <  1:1782515/150‑1 (MQ=255)
      ccTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGTTGCCTGGGGaaaa                                                                                                                                           <  1:1745462/150‑2 (MQ=255)
        tgtctgtAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGAc                                                                                                                                         <  2:318554/150‑1 (MQ=255)
        ccaccgcaGACACAACCTCTGGCGAACAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGAc                                                                                                                                         <  1:1627626/143‑1 (MQ=255)
         gtctgtAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGTGAAAGAcc                                                                                                                                        >  2:53107/1‑150 (MQ=255)
          tctgtAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGTTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCa                                                                                                                                       >  1:1070700/1‑150 (MQ=255)
            tgtAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCaaa                                                                                                                                     >  2:856653/1‑150 (MQ=255)
                 aCAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCggg                                                                                                                                <  1:2104102/150‑1 (MQ=255)
                  cAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGACAGACCAAACCGggc                                                                                                                               >  2:1559965/1‑150 (MQ=255)
                    aaaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggcgg                                                                                                                             <  2:1750857/150‑1 (MQ=255)
                            ggTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATa                                                                                                                     <  2:1593649/150‑1 (MQ=255)
                            ggTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATa                                                                                                                     <  1:864036/150‑1 (MQ=255)
                            ggTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATa                                                                                                                     <  2:2017169/150‑1 (MQ=255)
                            ggTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATa                                                                                                                     <  2:2017380/150‑1 (MQ=255)
                            ggTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGAGGAAAGACCAAACCGGGCGGCGACGATa                                                                                                                     <  1:1534723/150‑1 (MQ=255)
                                 aaaaCACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTg                                                                                                                >  1:1120192/1‑150 (MQ=255)
                                  aaaCACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGa                                                                                                               <  1:893837/150‑1 (MQ=255)
                                     cacCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGacac                                                                                                            <  2:2765095/150‑1 (MQ=255)
                                     cacCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGacac                                                                                                            <  2:2765118/150‑1 (MQ=255)
                                       ccTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGt                                                                                                          <  2:1255658/150‑1 (MQ=255)
                                       ccTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGt                                                                                                          >  2:2896073/1‑150 (MQ=255)
                                       ccTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGt                                                                                                          <  1:432766/150‑1 (MQ=255)
                                       ccTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGt                                                                                                          >  2:922640/1‑150 (MQ=255)
                                       ccTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGt                                                                                                          <  1:1502861/150‑1 (MQ=255)
                                           cGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCTACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTaa                                                                                                      <  2:1190667/150‑1 (MQ=255)
                                                  acaTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGc                                                                                               <  1:2565326/150‑1 (MQ=255)
                                                      cACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAACCCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTa                                                                                           >  2:441703/1‑150 (MQ=255)
                                                      cACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTa                                                                                           <  2:1826812/150‑1 (MQ=255)
                                                      cACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTa                                                                                           <  1:1065906/150‑1 (MQ=255)
                                                      cACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTa                                                                                           >  1:2308655/1‑150 (MQ=255)
                                                       aCTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAg                                                                                          <  2:2755882/150‑1 (MQ=255)
                                                        cTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGa                                                                                         <  2:1408637/150‑1 (MQ=255)
                                                           ccGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGt                                                                                      >  2:1756908/1‑150 (MQ=255)
                                                           ccGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGt                                                                                      >  2:1636272/1‑150 (MQ=255)
                                                           ccGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGt                                                                                      >  2:908285/1‑150 (MQ=255)
                                                           ccGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGGGCTGACGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAACGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGt                                                                                      >  2:163062/1‑150 (MQ=255)
                                                             gACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCa                                                                                    >  1:566952/1‑150 (MQ=255)
                                                             gACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCa                                                                                    <  1:1035380/150‑1 (MQ=255)
                                                                   tactacCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAg                                                                              <  2:1494865/150‑1 (MQ=255)
                                                                   tactacCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAg                                                                              <  2:2404858/150‑1 (MQ=255)
                                                                   tactacCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAg                                                                              <  2:2404865/150‑1 (MQ=255)
                                                                    actacCGCGGCCGCCAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAgg                                                                             <  2:365416/150‑1 (MQ=255)
                                                                    actacCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAgg                                                                             <  2:1927188/150‑1 (MQ=255)
                                                                     ctacCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAggg                                                                            <  2:2669560/150‑1 (MQ=255)
                                                                        ccGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGAt                                                                         >  1:864677/1‑150 (MQ=255)
                                                                           cgGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGAtttt                                                                      >  1:346468/1‑150 (MQ=255)
                                                                             gCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTgg                                                                    <  2:1120192/150‑1 (MQ=255)
                                                                              ccGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGc                                                                   <  1:951439/150‑1 (MQ=255)
                                                                              ccGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGc                                                                   <  1:2216634/150‑1 (MQ=255)
                                                                              ccGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGc                                                                   <  1:838893/150‑1 (MQ=255)
                                                                                gCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCACCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGccc                                                                 >  2:124844/1‑150 (MQ=255)
                                                                                    ggTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTcc                                                             >  1:1194349/1‑150 (MQ=255)
                                                                                       cATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTg                                                          <  1:1890206/150‑1 (MQ=255)
                                                                                        aTCGCTAAGTNATCANGCATNTGCGTGATGAAGCTNAGTNAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGa                                                         >  2:2467464/1‑150 (MQ=255)
                                                                                         tCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGAc                                                        >  2:679480/1‑150 (MQ=255)
                                                                                         tCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGAc                                                        >  2:1172256/1‑150 (MQ=255)
                                                                                          cGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACa                                                       >  2:1413874/1‑150 (MQ=255)
                                                                                           gCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCCACTGGGGAAAGACCACACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCAGGGGAGGGGATTTTGGCCCTTCCGTGACAg                                                      >  2:2518977/1‑150 (MQ=255)
                                                                                             tAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTg                                                    >  2:2054589/1‑150 (MQ=255)
                                                                                             tAAGGACTCACGCATCTGCGGGATGAAGCAGAGTGAGGATTATCCCCAGGCAACAGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGGTCCATGACCTCGGAGGGGAGTTAGGCCCTTCCGTGACAGTg                                                    >  2:1695112/1‑150 (MQ=255)
                                                                                              taGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGc                                                   <  2:2292655/149‑1 (MQ=255)
                                                                                                    tCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcag                                             >  1:1886073/1‑150 (MQ=255)
                                                                                                    tCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcag                                             >  1:1886059/1‑150 (MQ=255)
                                                                                                      aCGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCCAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcagca                                           >  1:242664/1‑150 (MQ=255)
                                                                                                      aCGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCCAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcagca                                           >  1:2571226/1‑150 (MQ=255)
                                                                                                       cGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCAt                                          <  1:344994/150‑1 (MQ=255)
                                                                                                       cGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCAt                                          <  1:690380/150‑1 (MQ=255)
                                                                                                        gCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCAGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcagcagt                                         >  2:1500034/1‑148 (MQ=255)
                                                                                                           tCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCa                                      >  2:1138083/1‑150 (MQ=255)
                                                                                                           tCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCa                                      <  2:2493927/150‑1 (MQ=255)
                                                                                                           tCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCa                                      <  2:516668/150‑1 (MQ=255)
                                                                                                           tCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCa                                      <  2:1273079/150‑1 (MQ=255)
                                                                                                            cTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAg                                     >  2:1005453/1‑150 (MQ=255)
                                                                                                              gCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGc                                   >  1:292146/1‑150 (MQ=255)
                                                                                                              gCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGc                                   >  1:455629/1‑150 (MQ=255)
                                                                                                               cGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCa                                  >  1:1028703/1‑150 (MQ=255)
                                                                                                               cGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAATGCTGCAGGCa                                  <  2:726612/150‑1 (MQ=255)
                                                                                                                gTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCAc                                 <  2:566952/150‑1 (MQ=255)
                                                                                                                    tgaAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGaa                             <  2:1744680/150‑1 (MQ=255)
                                                                                                                       aGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTc                          >  1:326900/1‑150 (MQ=255)
                                                                                                                          ttAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTATCCCTGGCGTTAGATGTCCGGGGAGGGGATTTTGGCCCTTCCGTGTCAGTGCCTGCAGCATTGCAGGCGCTGGATTCtct                       >  1:300882/1‑148 (MQ=255)
                                                                                                                               gAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGc                  >  1:1022768/1‑150 (MQ=255)
                                                                                                                               gAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGc                  >  1:2598709/1‑150 (MQ=255)
                                                                                                                               gAGGATTGTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGc                  <  1:916963/150‑1 (MQ=255)
                                                                                                                                   tctttCCCCAGTTGCCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATGTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAcct              >  2:2821314/3‑150 (MQ=255)
                                                                                                                                   tctttCCCCAGTTGCCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGGTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAact              >  2:1106980/3‑150 (MQ=255)
                                                                                                                                   aTTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAact              >  2:884554/1‑150 (MQ=255)
                                                                                                                                    ctttCCCCAGTTGCCTGGGGAAAGACCAAAACGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAactc             >  1:1415841/2‑150 (MQ=255)
                                                                                                                                     tttCCCCAGTTGCCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGCGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCCactca            >  2:2191107/1‑150 (MQ=255)
                                                                                                                                        ccccAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAactcactc         <  2:1922091/150‑1 (MQ=255)
                                                                                                                                           cAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACtcttc      <  1:713545/150‑1 (MQ=255)
                                                                                                                                              gCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCCACTCACTCTTCttt   >  1:1341594/1‑150 (MQ=255)
                                                                                                                                              gCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACtcttctgt   >  2:1010050/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTT  >  NC_000913/1147441‑1147733

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