Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,788,336 A→T 8.8% intergenic (+58/‑99) ppsR → / → aroH posphoenolpyruvate synthetase regulatory protein/3‑deoxy‑7‑phosphoheptulonate synthase, Trp‑sensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,788,3360AT 8.8% 380.8 / 18.4 135intergenic (+58/‑99)ppsR/aroHposphoenolpyruvate synthetase regulatory protein/3‑deoxy‑7‑phosphoheptulonate synthase, Trp‑sensitive
Reads supporting (aligned to +/- strand):  ref base A (62/59);  new base T (6/6);  total (70/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 142 total aligned reads displayed.

GATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACG  >  NC_000913/1788188‑1788485
                                                                                                                                                    |                                                                                                                                                     
gATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                                                                      <  1:283377/150‑1 (MQ=255)
 aTCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGtca                                                                                                                                                     <  2:930744/150‑1 (MQ=255)
  tCCCGTGGATTAACAGTTCCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGacac                                                                                                                                                    <  1:428975/150‑1 (MQ=255)
        ggATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGTACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCg                                                                                                                                              >  2:1981703/1‑150 (MQ=255)
         gATTAACAGTACCAATTATTCTGTAGAAGAGATTTCCACCAAGATCCTCGATATCATGTGCCTTATTCGCCGATTGTACTAGAGAAGTTGTGCATTATATTATTTTTTTGTTATCATTCTAACCACACGGCGAGGTGTGTCACACCTCGc                                                                                                                                             <  2:1772849/150‑1 (MQ=255)
         gATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGCCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGc                                                                                                                                             <  1:874937/150‑1 (MQ=255)
           ttAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTCACCACCCGGCGAGGTGTGTCACACCTCGccg                                                                                                                                           <  2:294469/150‑3 (MQ=255)
           ttAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccg                                                                                                                                           <  2:884064/150‑3 (MQ=255)
           ttAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCAc                                                                                                                                           >  2:1387525/1‑150 (MQ=255)
            tAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccgg                                                                                                                                          <  1:1913094/150‑4 (MQ=255)
            tAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACt                                                                                                                                          >  1:1235299/1‑150 (MQ=255)
             aaCAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACtt                                                                                                                                         >  2:866753/1‑150 (MQ=255)
             aaCAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACtt                                                                                                                                         >  1:1950298/1‑150 (MQ=255)
             aaCAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACtt                                                                                                                                         >  1:1915338/1‑150 (MQ=255)
             aaCAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACtt                                                                                                                                         >  1:1155378/1‑150 (MQ=255)
               cAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGa                                                                                                                                       >  1:2076037/1‑150 (MQ=255)
                aGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccgggtgg                                                                                                                                      <  1:820957/150‑8 (MQ=255)
                 gTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAAt                                                                                                                                    <  2:516098/150‑1 (MQ=255)
                  tACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAAt                                                                                                                                    <  1:358079/150‑1 (MQ=255)
                      aaTTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATaagc                                                                                                                                >  2:1525268/1‑150 (MQ=255)
                        ttattCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATcagcag                                                                                                                              <  2:1567435/150‑1 (MQ=255)
                         tattCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCagcagc                                                                                                                             >  2:467743/1‑150 (MQ=255)
                               gTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTATAGAACTAGTTCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGGGTGACACACCTCGCACTTGAAATCAGCAGCGATTgg                                                                                                                       <  2:1498720/150‑1 (MQ=255)
                                 agaagaGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGtt                                                                                                                     >  2:1705921/1‑150 (MQ=255)
                                 agaagaGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACCAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGtt                                                                                                                     >  2:2107283/1‑150 (MQ=255)
                                     gagaTTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATc                                                                                                                 >  2:1857673/1‑150 (MQ=255)
                                      agaTTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATc                                                                                                                                   >  2:1982534/1‑131 (MQ=255)
                                          tGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGAt                                                                                                            >  1:2115040/1‑150 (MQ=255)
                                           gCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATg                                                                                                           >  2:2177963/1‑150 (MQ=255)
                                           gCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATg                                                                                                           >  2:2098989/1‑150 (MQ=255)
                                           gCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATg                                                                                                           <  2:78036/150‑1 (MQ=255)
                                            ccaccaAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATgc                                                                                                          <  1:1105743/150‑1 (MQ=255)
                                               ccaAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCa                                                                                                       <  2:465533/150‑1 (MQ=255)
                                               ccaAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCa                                                                                                       <  1:1464916/150‑1 (MQ=255)
                                               ccaAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCa                                                                                                       <  1:1769598/150‑1 (MQ=255)
                                               ccaAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCa                                                                                                       <  1:1778512/150‑1 (MQ=255)
                                                   gATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCAc                                                                                                   >  1:1909445/1‑150 (MQ=255)
                                                      ccTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCAcgat                                                                                                <  2:1155378/150‑4 (MQ=255)
                                                        tCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATACTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTccc                                                                                              <  2:28954/150‑1 (MQ=255)
                                                             atCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAg                                                                                         <  2:1145471/150‑1 (MQ=255)
                                                             atCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAg                                                                                         <  1:189633/150‑1 (MQ=255)
                                                              tCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGt                                                                                        <  1:810241/150‑1 (MQ=255)
                                                                aTGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTGACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTcc                                                                                      <  2:1303614/150‑1 (MQ=255)
                                                                    gCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGcc                                                                                  <  2:901731/150‑1 (MQ=255)
                                                                      cTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGt                                                                                >  1:1338978/1‑150 (MQ=255)
                                                                         aGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTCATGCGCATCACTTCCCGGCAGTCCTGCCGTAGa                                                                             >  2:2059712/1‑150 (MQ=255)
                                                                          gTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCTCGGCAGGACTGCCgggaag                                                                            <  1:2076539/150‑6 (MQ=255)
                                                                          gTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGaa                                                                            <  1:2141317/150‑1 (MQ=255)
                                                                          gTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGGACTGCCgggaag                                                                            <  1:714474/150‑6 (MQ=255)
                                                                           gcGCCGAATGTACTCGAGAACTCGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAg                                                                           <  2:1908993/149‑1 (MQ=255)
                                                                            cGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGc                                                                          >  1:519997/1‑150 (MQ=255)
                                                                             gCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGca                                                                         <  2:517925/150‑1 (MQ=255)
                                                                               cGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGcaac                                                                       <  1:693464/150‑1 (MQ=255)
                                                                                  aTGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACaaa                                                                    <  1:445717/150‑1 (MQ=255)
                                                                                  aTGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACaaa                                                                    <  1:240298/150‑1 (MQ=255)
                                                                                          gagaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTgag                                                            <  1:467743/150‑1 (MQ=255)
                                                                                            taaCTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGAc                                                          <  1:1090503/149‑1 (MQ=255)
                                                                                            gaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGAc                                                          >  1:1296416/1‑150 (MQ=255)
                                                                                            gaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGAc                                                          <  1:1884765/150‑1 (MQ=255)
                                                                                            gaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGAc                                                          >  1:1287950/1‑150 (MQ=255)
                                                                                            gaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGAc                                                          >  1:255165/1‑150 (MQ=255)
                                                                                            gaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGAc                                                          <  1:866753/150‑1 (MQ=255)
                                                                                               cTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTg                                                       >  2:827002/1‑150 (MQ=255)
                                                                                               cTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACGTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTg                                                       <  2:393735/150‑1 (MQ=255)
                                                                                                 aGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTa                                                     >  1:1618934/1‑150 (MQ=255)
                                                                                                 aGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTa                                                     >  2:2444674/1‑150 (MQ=255)
                                                                                                 aGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTa                                                     <  2:2318796/150‑1 (MQ=255)
                                                                                                   tGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAAt                                                   >  1:1248018/1‑150 (MQ=255)
                                                                                                    gCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCGCTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATg                                                  <  2:1015064/150‑1 (MQ=255)
                                                                                                          gCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGa                                            >  1:2308521/1‑150 (MQ=255)
                                                                                                            ttATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAAc                                          >  2:2469342/1‑150 (MQ=255)
                                                                                                             tATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACt                                         >  1:2172877/1‑150 (MQ=255)
                                                                                                              aTTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTg                                        <  1:2424069/150‑1 (MQ=255)
                                                                                                               tttttttGTTATCATGCTAACCACCCGGCGCGGTGTGACACCCCTCGCACTTGAAATCAGCAGCGCTTGGTTTATCGTGATGCGCATCCCTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGCACAGAACTgc                                       >  2:1029188/1‑149 (MQ=255)
                                                                                                               tttttttGTTATCATGCTAACCACCCGGCGAGGTGTGCCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGa                                       >  2:2417595/1‑150 (MQ=255)
                                                                                                                ttttttGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGAc                                      <  2:806846/150‑1 (MQ=255)
                                                                                                                         aTCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGACCTCCGt                             >  2:1581239/1‑150 (MQ=255)
                                                                                                                         aTCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGt                             >  2:1291061/1‑150 (MQ=255)
                                                                                                                           cATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTAc                           >  1:1366233/1‑150 (MQ=255)
                                                                                                                           cATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTAc                           <  2:2343557/150‑1 (MQ=255)
                                                                                                                                    caccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTAt                  >  1:39400/1‑150 (MQ=255)
                                                                                                                                    caccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTAt                  >  1:655746/1‑150 (MQ=255)
                                                                                                                                     caggtgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTAtt                 >  1:1028188/6‑150 (MQ=255)
                                                                                                                                     accCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTAtt                 >  2:1463026/1‑150 (MQ=255)
                                                                                                                                       ccGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTga               <  2:1239027/150‑1 (MQ=255)
                                                                                                                                        gtgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTgag              >  2:137769/3‑150 (MQ=255)
                                                                                                                                        gtgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTgag              >  1:1970607/3‑150 (MQ=255)
                                                                                                                                         tgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTgaga             >  1:2287988/2‑150 (MQ=255)
                                                                                                                                         ggCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCTCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTgaga             <  1:611552/150‑1 (MQ=255)
                                                                                                                                         ggCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTgaga             >  2:601152/1‑150 (MQ=255)
                                                                                                                                            gAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGcc          <  1:1463026/150‑1 (MQ=255)
                                                                                                                                            gAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGcc          <  1:997782/150‑1 (MQ=255)
                                                                                                                                             aGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCt         <  2:655746/150‑1 (MQ=255)
                                                                                                                                              ggTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTg        >  2:1363557/1‑150 (MQ=255)
                                                                                                                                               gtgtGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTgg       <  1:601152/150‑1 (MQ=255)
                                                                                                                                                tgtgACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGCACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGt      >  2:1683701/1‑150 (MQ=255)
                                                                                                                                                 gtgAAACACCTCGCACTTGAAATCAGGAGCGATTGGTTTATAGTGATGCGAATCACTTCCCGGCAGTCCTTCCGTAGAAGCAACAAATTTCTGATACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTa     <  2:1338978/150‑1 (MQ=255)
                                                                                                                                                  tgACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTaa    <  1:409238/150‑1 (MQ=255)
                                                                                                                                                   tacacacCTCGCACTTGAAATCAGAAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAAc   <  1:587240/149‑1 (MQ=255)
                                                                                                                                                   gACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAAc   >  2:2239404/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACG  >  NC_000913/1788188‑1788485

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