Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,788,336 A→T 5.5% intergenic (+58/‑99) ppsR → / → aroH PEP synthase kinase and PEP synthase pyrophosphorylase/3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,788,3360AT 5.5% 391.6 / 9.7 109intergenic (+58/‑99)ppsR/aroHPEP synthase kinase and PEP synthase pyrophosphorylase/3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base A (49/54);  new base T (3/3);  total (52/57)
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.83e-01

Only 100 of 113 total aligned reads displayed.

ATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCT  >  NC_000913/1788189‑1788478
                                                                                                                                                   |                                                                                                                                              
aTCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGac                                                                                                                                               >  3:278901/1‑149 (MQ=255)
aTCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGa                                                                                                                                                >  1:317148/1‑148 (MQ=255)
   ccGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGacaca                                                                                                                                            <  4:726194/149‑1 (MQ=255)
     gTGGATTAACTGTACCAACTATTCGGTAGAAGAGATTGCCACCAAGATCCTCCATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTCGCTTATTTTTTTGTTATCATGCTAAGCACCCGGCGAGGTGTGacacag                                                                                                                                           >  2:659025/1‑147 (MQ=255)
       ggATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCt                                                                                                                                         <  1:729317/148‑1 (MQ=255)
           tAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCa                                                                                                                                     <  1:722249/148‑1 (MQ=255)
            aaCAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGTCACACCTCGccgg                                                                                                                                   <  3:25077/149‑4 (MQ=255)
                 tACCAATTATTCGGTAGAAGAGATTGCCACCAAGGTCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGaa                                                                                                                               >  2:104047/1‑148 (MQ=255)
                 tACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGaaa                                                                                                                              <  3:466619/149‑1 (MQ=255)
                 tACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGaa                                                                                                                               >  2:419972/1‑148 (MQ=255)
                 tACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGaa                                                                                                                               <  1:277711/148‑1 (MQ=255)
                 tACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGaa                                                                                                                               >  2:736549/1‑148 (MQ=255)
                   ccAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATc                                                                                                                            >  4:383648/1‑149 (MQ=255)
                    cAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGAGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATca                                                                                                                           >  4:650803/1‑149 (MQ=255)
                     aaTTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATcag                                                                                                                          <  3:52531/149‑1 (MQ=255)
                     aaTTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATca                                                                                                                           >  2:525104/1‑148 (MQ=255)
                        tattCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATcagcag                                                                                                                       <  3:386641/149‑1 (MQ=255)
                            cGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGa                                                                                                                    <  2:317148/148‑1 (MQ=255)
                             ggTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGAtt                                                                                                                  >  3:147561/1‑149 (MQ=255)
                               tAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTgg                                                                                                                <  4:499374/149‑1 (MQ=255)
                               tAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTg                                                                                                                 <  1:736549/148‑1 (MQ=255)
                               tAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTg                                                                                                                 <  1:419972/148‑1 (MQ=255)
                                agaagaGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTgg                                                                                                                <  2:586712/148‑1 (MQ=255)
                                agaagaGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGt                                                                                                               >  3:363588/1‑149 (MQ=255)
                                agaagaGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGt                                                                                                               <  4:680725/149‑1 (MQ=255)
                                  aagaGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATca                                                                                                                           >  2:132180/1‑135 (MQ=255)
                                  aagaGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATca                                                                                                                           <  1:132180/135‑1 (MQ=255)
                                     agaTTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTAt                                                                                                           <  2:174070/148‑1 (MQ=255)
                                     agaTTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTAt                                                                                                           <  1:250721/148‑1 (MQ=255)
                                     agaTTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATc                                                                                                          <  4:390740/149‑1 (MQ=255)
                                     agaTTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCCCGAACTTGAAATCAGCAGCGATTGGTTTATc                                                                                                          <  4:555061/149‑1 (MQ=255)
                                      gaTTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCg                                                                                                         >  4:825246/1‑149 (MQ=255)
                                         tGCCACCAAGAACCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAACTCAGCAGCGATTGGTTTATCGTGa                                                                                                      >  3:595077/1‑149 (MQ=255)
                                               caAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATgcgc                                                                                                 <  1:57265/148‑1 (MQ=255)
                                                   aTCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCAc                                                                                            >  3:705083/1‑149 (MQ=255)
                                                   aTCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGCTTATCGTGATGCGCATCa                                                                                             >  2:758467/1‑148 (MQ=255)
                                                      cTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATcagcag                                                                                                                       <  3:236043/119‑1 (MQ=255)
                                                      cTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATcagcag                                                                                                                       >  4:236043/1‑119 (MQ=255)
                                                      cTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATc                                                                                                          <  3:38219/132‑1 (MQ=255)
                                                      cTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTc                                                                                         >  3:21101/1‑149 (MQ=255)
                                                      cTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTc                                                                                         >  4:26502/1‑149 (MQ=255)
                                                        cGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTcc                                                                                        >  1:638306/1‑148 (MQ=255)
                                                         gATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTccc                                                                                       <  1:727108/148‑1 (MQ=255)
                                                          atatCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCg                                                                                      <  2:638306/148‑1 (MQ=255)
                                                               aTGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGt                                                                                 >  2:798417/1‑148 (MQ=255)
                                                                 gggCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGt                                                                                 >  2:654327/1‑146 (MQ=255)
                                                                 gggCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGt                                                                                 <  1:654327/146‑1 (MQ=255)
                                                                   gCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGc                                                                            <  4:375909/149‑1 (MQ=255)
                                                                         gTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGGACTGCCgggaa                                                                      <  3:374541/149‑5 (MQ=255)
                                                                          tCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGaa                                                                     <  4:595077/149‑1 (MQ=255)
                                                                           cGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATca                                                                                                                           >  4:465718/1‑94 (MQ=255)
                                                                           cGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATca                                                                                                                           <  3:465718/94‑1 (MQ=255)
                                                                              cGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGca                                                                  <  1:32788/148‑1 (MQ=255)
                                                                              cGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGca                                                                  <  2:770732/148‑1 (MQ=255)
                                                                                       tAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTc                                                         <  1:798417/148‑1 (MQ=255)
                                                                                       tAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCt                                                        >  3:785241/1‑149 (MQ=255)
                                                                                        agagAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCt                                                        >  2:50785/1‑148 (MQ=255)
                                                                                         gagaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTga                                                      <  4:278901/149‑1 (MQ=255)
                                                                                         gagaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTga                                                      <  3:330427/149‑1 (MQ=255)
                                                                                         gagaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTga                                                      >  3:754886/1‑149 (MQ=255)
                                                                                         gagaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTg                                                       <  1:525104/148‑1 (MQ=255)
                                                                                           gaACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCGGTAGAAGCAACAAATTTCTgag                                                     >  1:753249/1‑148 (MQ=255)
                                                                                              cTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACtt                                                 <  3:632545/149‑1 (MQ=255)
                                                                                              cTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACt                                                  <  2:92870/148‑1 (MQ=255)
                                                                                               tAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTg                                                >  3:336324/1‑149 (MQ=255)
                                                                                               tAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTg                                                >  3:600468/1‑149 (MQ=255)
                                                                                                aGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTg                                                >  2:81013/1‑148 (MQ=255)
                                                                                                     aTTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGaaa                                                                                                                              <  3:142103/65‑1 (MQ=255)
                                                                                                     aTTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGaaa                                                                                                                              >  4:142103/1‑65 (MQ=255)
                                                                                                     aTTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGa                                          >  3:215097/1‑149 (MQ=255)
                                                                                                      ttAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGaa                                         <  4:754886/149‑1 (MQ=255)
                                                                                                      ttAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGaa                                         <  4:785241/149‑1 (MQ=255)
                                                                                                      ttAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGa                                          <  1:837309/148‑1 (MQ=255)
                                                                                                           ttATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACTCCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGTGACTTGTAATGAACAGa                                     <  2:621354/148‑1 (MQ=255)
                                                                                                           ttATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGaa                                         <  3:279725/144‑1 (MQ=255)
                                                                                                           ttATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGaa                                         >  4:279725/1‑144 (MQ=255)
                                                                                                           ttATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGaa                                    >  4:46706/1‑149 (MQ=255)
                                                                                                           ttATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGa                                     >  1:487507/1‑148 (MQ=255)
                                                                                                             aTTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACt                                  <  4:21101/149‑1 (MQ=255)
                                                                                                                  tttGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACg                              >  1:638608/1‑148 (MQ=255)
                                                                                                                      ttATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTga                                                      >  3:374378/1‑120 (MQ=255)
                                                                                                                      ttATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTga                                                      <  4:374378/120‑1 (MQ=255)
                                                                                                                       tATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTc                         <  1:50785/148‑1 (MQ=255)
                                                                                                                         tCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGt                      >  4:698931/1‑149 (MQ=255)
                                                                                                                                 accaccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTgcgc               >  1:914011/1‑148 (MQ=255)
                                                                                                                                  ccaccCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGcgcg              >  2:770806/1‑148 (MQ=255)
                                                                                                                                    accCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGaa                                                                     <  4:593574/91‑1 (MQ=255)
                                                                                                                                    accCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGaa                                                                     >  3:593574/1‑91 (MQ=255)
                                                                                                                                       cGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTg         <  1:104047/148‑1 (MQ=255)
                                                                                                                                       cGGCGAGGTGTGACACACCTCACACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTg         >  2:479208/1‑148 (MQ=255)
                                                                                                                                        tgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTgaga                                                    <  4:129471/103‑1 (MQ=255)
                                                                                                                                        tgCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTgaga                                                    >  3:129471/2‑104 (MQ=255)
                                                                                                                                         gCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTgag       >  2:814733/1‑148 (MQ=255)
                                                                                                                                         gCGAGGTGTGTCACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTgag       >  2:350784/1‑148 (MQ=255)
                                                                                                                                         gCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGaa                            >  3:84317/1‑127 (MQ=255)
                                                                                                                                           gAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGc    <  3:741428/149‑1 (MQ=255)
                                                                                                                                             ggTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCt  <  4:26826/149‑1 (MQ=255)
                                                                                                                                              gtgtGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCt  <  1:770806/148‑1 (MQ=255)
                                                                                                                                                gtgACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGGACTGCCggga                                                                       >  1:35775/1‑74 (MQ=255)
                                                                                                                                                gtgACACACCTCGCACTTGAAAACAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGGACTGCCggga                                                                       <  2:35775/77‑4 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCCTCGATATCATGGGCCTTAGTCGCCGAATGTACTAGAGAACTAGTGCATTAGCTTATTTTTTTGTTATCATGCTAACCACCCGGCGAGGTGTGACACACCTCGCACTTGAAATCAGCAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCT  >  NC_000913/1788189‑1788478

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