Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,830 G→A 9.6% intergenic (‑263/+23) bglG ← / ← phoU transcriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,8300GA 9.6% 277.2 / 16.4 95intergenic (‑263/+23)bglG/phoUtranscriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator
Reads supporting (aligned to +/- strand):  ref base G (48/37);  new base A (7/2);  total (55/40)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

CGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACG  >  NC_000913/3906683‑3906978
                                                                                                                                                    |                                                                                                                                                    
cGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGt                                                                                                                                                     >  3:381150/1‑149 (MQ=255)
 gAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGtt                                                                                                                                                    >  1:628761/1‑149 (MQ=255)
    cATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTAtt                                                                                                                                                 <  2:541665/149‑1 (MQ=255)
     aTCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTg                                                                                                                                                >  2:246456/1‑149 (MQ=255)
       ccAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTggg                                                                                                                                              >  1:519212/1‑149 (MQ=255)
        cAGGTTAGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGAGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAAAATTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGa                                                                                                                                             <  2:300168/149‑1 (MQ=255)
         aGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGAt                                                                                                                                            <  1:484362/149‑1 (MQ=255)
         aGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGAt                                                                                                                                            >  1:49948/1‑149 (MQ=255)
           gTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGAttt                                                                                                                                          <  3:375955/149‑1 (MQ=255)
               ggAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCt                                                                                                                                      <  1:333518/149‑1 (MQ=255)
                gAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTg                                                                                                                                     >  2:73658/1‑149 (MQ=255)
                  aattaattTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGt                                                                                                                                   <  1:457811/149‑1 (MQ=255)
                  aattaattTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGt                                                                                                                                   <  4:257607/149‑1 (MQ=255)
                    ttaatttaatTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGa                                                                                                                                 >  3:242727/1‑149 (MQ=255)
                        tttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAatta                                                                                                                             <  4:170556/149‑1 (MQ=255)
                             ttGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGt                                                                                                                        <  1:697730/149‑1 (MQ=255)
                                  tGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTa                                                                                                                   <  4:537934/149‑1 (MQ=255)
                                  tGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTa                                                                                                                   >  3:374558/1‑149 (MQ=255)
                                   gCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTAt                                                                                                                  >  1:455023/1‑149 (MQ=255)
                                     agagAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCt                                                                                                                >  1:371057/1‑149 (MQ=255)
                                      gagaAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCtt                                                                                                               >  1:410339/1‑149 (MQ=255)
                                        gaAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTc                                                                                                             <  3:408695/149‑1 (MQ=255)
                                           aaTAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTcccc                                                                                                          >  3:642600/1‑149 (MQ=255)
                                                cAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcag                                                                                                     >  2:94837/1‑149 (MQ=255)
                                                  aTGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagca                                                                                                   <  1:77930/149‑1 (MQ=255)
                                                   tGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagcag                                                                                                  >  3:658923/1‑149 (MQ=255)
                                                   tGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagcag                                                                                                  >  3:27235/1‑149 (MQ=255)
                                                   tGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagcag                                                                                                  >  3:25783/1‑149 (MQ=255)
                                                     cgcTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGtt                                                                                                <  2:485478/149‑1 (MQ=255)
                                                        tATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTAt                                                                                             >  4:349924/1‑149 (MQ=255)
                                                        tATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTAt                                                                                             >  3:754198/1‑149 (MQ=255)
                                                         aTTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATc                                                                                            <  1:737194/149‑1 (MQ=255)
                                                            gATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAg                                                                                         >  1:302117/1‑149 (MQ=255)
                                                                  aaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATc                                                                                   >  2:446134/1‑149 (MQ=255)
                                                                          cATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGa                                                                           >  3:271260/1‑149 (MQ=255)
                                                                           aTGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGAc                                                                          <  2:771122/149‑1 (MQ=255)
                                                                               tCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGa                                                                      <  3:584076/149‑1 (MQ=255)
                                                                               tCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGa                                                                      >  4:402736/1‑149 (MQ=255)
                                                                               tCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGa                                                                      >  1:594801/1‑149 (MQ=255)
                                                                                 gCAGTTATTAACTTTGTGTAATTTTAGGAATTTCTAAAGCTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGCTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCGAGCTCATCGCCACCTACGTGACg                                                                    >  3:410642/1‑149 (MQ=255)
                                                                                 gCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGAc                                                                     >  4:135548/1‑148 (MQ=255)
                                                                                 gCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGAc                                                                     <  3:135548/148‑1 (MQ=255)
                                                                                 gCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACg                                                                    <  2:65792/149‑1 (MQ=255)
                                                                                    gTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGaa                                                                 <  3:277240/149‑1 (MQ=255)
                                                                                     ttattaACTTTGTGTAAGTTTAGGAATGTATAAAGTTATATATAACAAATCCCAATAATTAAGTTACTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAATTTATCCAGCTCATCGCCACCCACGTGACGGaaa                                                                <  4:640263/149‑1 (MQ=255)
                                                                                        ttaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATcc                                                             >  4:519260/1‑149 (MQ=255)
                                                                                         taACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCt                                                            <  1:712207/149‑1 (MQ=255)
                                                                                          aaCTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTg                                                           >  1:189547/1‑149 (MQ=255)
                                                                                           aCTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGc                                                          <  2:594801/149‑1 (MQ=255)
                                                                                             tttGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAAGTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTAt                                                                                                                  <  3:462271/91‑1 (MQ=255)
                                                                                             tttGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCACAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTAt                                                                                                                  >  4:462271/1‑91 (MQ=255)
                                                                                                 tgtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGACATCCTGCCCCTTc                                                    >  4:263904/1‑149 (MQ=255)
                                                                                                 tgtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTc                                                    >  4:433192/1‑149 (MQ=255)
                                                                                                   tAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCAc                                                  >  1:706967/1‑149 (MQ=255)
                                                                                                   tAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCAc                                                  >  3:648958/1‑149 (MQ=255)
                                                                                                   tAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCAc                                                  <  2:49948/149‑1 (MQ=255)
                                                                                                       tttAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACgtag                                              >  2:444519/1‑149 (MQ=255)
                                                                                                            gAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                         >  2:421560/1‑149 (MQ=255)
                                                                                                            gAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                         <  3:433192/149‑1 (MQ=255)
                                                                                                            gAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAACTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                         <  1:73658/149‑1 (MQ=255)
                                                                                                               tttATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGAt                                      >  3:114503/1‑149 (MQ=255)
                                                                                                                 tataAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATaa                                    >  1:18758/1‑149 (MQ=255)
                                                                                                                  ataAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGTAAATCCTGCCCCTTCACGTAGTAGAAGATaaa                                   >  2:326720/1‑149 (MQ=255)
                                                                                                                   taAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAAc                                  >  2:140668/1‑149 (MQ=255)
                                                                                                                     aaGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGACACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTc                                <  4:164411/149‑1 (MQ=255)
                                                                                                                      aGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTca                               <  4:658923/149‑1 (MQ=255)
                                                                                                                      aGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTca                               <  4:27235/149‑1 (MQ=255)
                                                                                                                        ttATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTAttt                                                                                                                          <  3:86495/56‑1 (MQ=255)
                                                                                                                        ttATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTAttt                                                                                                                          >  4:86495/1‑56 (MQ=255)
                                                                                                                           tatataACAAATCCCAATAATTAAGATATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAat                          >  3:470922/1‑149 (MQ=255)
                                                                                                                              ataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATAtt                       >  3:580210/1‑149 (MQ=255)
                                                                                                                              ataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATAtt                       <  4:615677/149‑1 (MQ=255)
                                                                                                                               taACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTc                      <  4:159220/149‑1 (MQ=255)
                                                                                                                               taACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTc                      <  2:455023/149‑1 (MQ=255)
                                                                                                                               taACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGACCCTTCACGTAGTAGAAGATAAACTCACAAATATTc                      >  4:513142/1‑149 (MQ=255)
                                                                                                                               agaCAAATCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGAc                                                                          >  3:677350/3‑97 (MQ=37)
                                                                                                                               agaCAAATCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGAc                                                                          <  4:677350/95‑1 (MQ=38)
                                                                                                                                aaCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCt                     >  2:353675/1‑149 (MQ=255)
                                                                                                                                aaCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCt                     <  3:349924/149‑1 (MQ=255)
                                                                                                                                   aaaTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGc                                                                                        >  4:321065/1‑79 (MQ=255)
                                                                                                                                   aaaTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGc                                                                                        <  3:321065/79‑1 (MQ=255)
                                                                                                                                     aTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAg                >  4:244913/1‑149 (MQ=255)
                                                                                                                                     aTCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAg                >  2:135031/1‑149 (MQ=255)
                                                                                                                                      tCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGc               >  2:92734/1‑149 (MQ=255)
                                                                                                                                      tCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCACCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATGAACTCACAAATATTCTGGCAGc               >  1:663737/1‑149 (MQ=255)
                                                                                                                                       cccAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGGAGATAAACTCACAAATATTCTGGCAGCg              >  2:623155/1‑149 (MQ=255)
                                                                                                                                       cccAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCg              >  3:442260/1‑149 (MQ=255)
                                                                                                                                        ccAATAATTAAGTTATGGGGATTTGGCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGGTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCgg             >  4:404121/1‑149 (MQ=255)
                                                                                                                                        ccAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                         <  2:314881/121‑1 (MQ=255)
                                                                                                                                        ccAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                         >  1:314881/1‑121 (MQ=255)
                                                                                                                                          aataatTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTc           >  3:56330/1‑149 (MQ=255)
                                                                                                                                          aataatTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTc           <  2:628761/149‑1 (MQ=255)
                                                                                                                                            taatTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGc         <  3:513142/149‑1 (MQ=255)
                                                                                                                                              atTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCa       >  3:397922/1‑149 (MQ=255)
                                                                                                                                                   aGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATAc   <  4:374558/149‑1 (MQ=255)
                                                                                                                                                    gTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACg  >  4:4568/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACG  >  NC_000913/3906683‑3906978

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