Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,830 G→A 7.4% 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 7.4% 505.4 / 18.4 164intergenic (‑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 (69/81);  new base A (4/8);  total (73/91)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

CGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATA  >  NC_000913/3906683‑3906976
                                                                                                                                                   |                                                                                                                                                  
cGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGt                                                                                                                                                   >  1:267069/1‑149 (MQ=255)
 gAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGtt                                                                                                                                                  >  2:1440935/1‑149 (MQ=255)
    cATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTAtt                                                                                                                                               <  1:1319789/149‑1 (MQ=255)
         aGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGAt                                                                                                                                          <  2:831280/149‑1 (MQ=255)
         aGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGAt                                                                                                                                          <  2:93717/149‑1 (MQ=255)
           gTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGAttt                                                                                                                                        <  1:538839/149‑1 (MQ=255)
               ggAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCt                                                                                                                                    >  1:246232/1‑149 (MQ=255)
               ggAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCt                                                                                                                                    <  2:801189/149‑1 (MQ=255)
               ggAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGAGAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGTATTTGTCt                                                                                                                                    >  1:210755/1‑149 (MQ=255)
                gAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTg                                                                                                                                   <  2:1470433/149‑1 (MQ=255)
                 aaaTTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACAATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTgt                                                                                                                                  >  1:324963/1‑148 (MQ=255)
                   attaattTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTg                                                                                                                                >  2:423630/1‑149 (MQ=255)
                        tttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAatta                                                                                                                           <  1:164860/149‑1 (MQ=255)
                        tttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGATCGCATTTATTAAATTTGTGTAATTTTAGTAATTTATAAAGTTATATATAACAAATCCCAATAATTAAATTATTTTTATTTGTCTTTTTAatta                                                                                                                           >  1:99254/1‑149 (MQ=255)
                         ttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAattat                                                                                                                          >  1:556203/1‑149 (MQ=255)
                         ttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAattat                                                                                                                          <  2:777273/149‑1 (MQ=255)
                          taattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAattatt                                                                                                                         <  1:761766/149‑1 (MQ=255)
                           aattGCGTGCAGAGAAAATAGCAATGCTCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTAttt                                                                                                                        <  2:68616/149‑1 (MQ=255)
                            attGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTg                                                                                                                       <  1:821192/149‑1 (MQ=255)
                             tttcGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATCCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGt                                                                                                                      <  2:275380/146‑1 (MQ=255)
                             ttGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGt                                                                                                                      <  1:1265842/149‑1 (MQ=255)
                                  tGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTa                                                                                                                 >  2:668252/1‑149 (MQ=255)
                                     agagAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCt                                                                                                              >  1:776515/1‑149 (MQ=255)
                                       agaAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCttt                                                                                                            >  1:642195/1‑149 (MQ=255)
                                        gaAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTc                                                                                                           >  2:1599085/1‑149 (MQ=255)
                                        gaAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTc                                                                                                           >  2:1290322/1‑149 (MQ=255)
                                         aaaaTAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTcc                                                                                                          >  2:307259/1‑149 (MQ=255)
                                         aaaaTAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTcc                                                                                                          <  2:1333285/149‑1 (MQ=255)
                                             tAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGc                                                                                                      <  2:1047936/149‑1 (MQ=255)
                                                cAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcag                                                                                                   <  1:1080116/149‑1 (MQ=255)
                                                 aaTGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGTATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagc                                                                                                  >  2:588655/1‑149 (MQ=255)
                                                   tGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagcag                                                                                                <  2:1351638/149‑1 (MQ=255)
                                                   tGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagcag                                                                                                >  1:1584139/1‑149 (MQ=255)
                                                     cgcTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGtt                                                                                              <  2:652125/149‑1 (MQ=255)
                                                       cTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTa                                                                                            <  2:45144/149‑1 (MQ=255)
                                                           tGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCa                                                                                        <  2:594738/149‑1 (MQ=255)
                                                             aTAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAAAAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGATTTATTTTTCGATATATTTCCCCGCAAGCAGTTTATCCAGc                                                                                      >  2:1279296/1‑149 (MQ=255)
                                                               aaaaaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTc                                                                                    >  2:955522/1‑149 (MQ=255)
                                                                aaaaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCa                                                                                   >  1:1185598/1‑149 (MQ=255)
                                                                 aaaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCAt                                                                                  >  2:358147/1‑149 (MQ=255)
                                                                  aaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATc                                                                                 <  1:1589951/149‑1 (MQ=255)
                                                                  aaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATc                                                                                 <  2:3564/149‑1 (MQ=255)
                                                                   atatGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCg                                                                                >  2:20910/1‑149 (MQ=255)
                                                                        aCCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCAcc                                                                           <  1:456921/149‑1 (MQ=255)
                                                                          cATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGa                                                                         >  1:1436243/1‑149 (MQ=255)
                                                                               tCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGa                                                                    >  2:179758/1‑149 (MQ=255)
                                                                               tCGCAGTTATTAACTTTGTGTAATTTTAGGAATATATAAAGTTATATATAACAAATCCCAAAAAATAACTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGa                                                                    <  1:1309996/149‑1 (MQ=255)
                                                                                 gCAGTTATTAACTTTGTGTAATTTTAGGAATTTGTAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACg                                                                  >  2:372633/1‑149 (MQ=255)
                                                                                 gCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCAAATAATTAAGTTATGGGGATTTGTCTGGTGAATTATTTTTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCCTCTCCACCGACGTGACg                                                                  >  2:264300/1‑149 (MQ=255)
                                                                                     ttattaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGaaa                                                              <  1:459501/149‑1 (MQ=255)
                                                                                      tattaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAAt                                                             <  1:1199298/149‑1 (MQ=255)
                                                                                       attaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATc                                                            <  1:1440935/149‑1 (MQ=255)
                                                                                       attaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATc                                                            >  1:1355321/1‑149 (MQ=255)
                                                                                         taACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCt                                                          <  1:842468/149‑1 (MQ=255)
                                                                                           aCTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGc                                                        <  2:1513892/149‑1 (MQ=255)
                                                                                                gtgtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCtt                                                   <  2:635853/149‑1 (MQ=255)
                                                                                                  gtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCa                                                 <  1:668252/149‑1 (MQ=255)
                                                                                                  gtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCa                                                 <  1:307259/149‑1 (MQ=255)
                                                                                                      ttttAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACgta                                             <  1:1381077/149‑1 (MQ=255)
                                                                                                        ttAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACgtagt                                           <  2:1014252/149‑1 (MQ=255)
                                                                                                            gAATTTATAAAGTTATATATAACAAATCCCAATTTTTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCGCCGCCAGCAGTTTATCAAGCTCATCGCCACAGACGTGACGGAAATCCTGCCCCTTCACGTAGTagat                                       >  1:34964/1‑148 (MQ=255)
                                                                                                            gAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGCCGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                       >  2:234978/1‑149 (MQ=255)
                                                                                                            gAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                       <  2:776515/149‑1 (MQ=255)
                                                                                                            gAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                       >  1:671557/1‑149 (MQ=255)
                                                                                                            gAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCAAGCTCATCGCAACGGTCGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                       >  1:29775/1‑149 (MQ=255)
                                                                                                             aaTTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaag                                      <  1:914628/149‑1 (MQ=255)
                                                                                                              aTTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCCGCTCATCGCCACCGACGTGACGCAAATCCTGCCCCTTCACGTAGTagaaga                                     >  2:118474/1‑149 (MQ=255)
                                                                                                               tttATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGAt                                    <  2:1406443/149‑1 (MQ=255)
                                                                                                               tttATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGAt                                    <  1:595402/149‑1 (MQ=255)
                                                                                                                  ataAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATaaa                                 <  1:714934/149‑1 (MQ=255)
                                                                                                                   taAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAAc                                >  1:525692/1‑149 (MQ=255)
                                                                                                                   taAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAAc                                >  1:1579025/1‑149 (MQ=255)
                                                                                                                   taAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACAGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAAc                                <  1:1033079/149‑1 (MQ=255)
                                                                                                                   taAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGCTGAATTATTTTTCGATATCTTACCCAGCCAGCAGTTTATCCAGATCATAGACACAGACGTGACTGAAATCCTGACCCTTCACGTAATataaaataaaa                                >  2:218764/1‑148 (MQ=255)
                                                                                                                    aaaGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACt                               >  2:659069/1‑149 (MQ=255)
                                                                                                                          atatataACAAATCCCAATAATTAAGTTTTTGGGTTTTGTCTGGTGAATTATTTGTCGATATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACaaa                         <  2:975998/149‑1 (MQ=255)
                                                                                                                           tatataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAat                        <  2:890786/149‑1 (MQ=255)
                                                                                                                            atataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAata                       >  2:1022285/1‑149 (MQ=255)
                                                                                                                              ataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATAtt                     <  2:324963/149‑1 (MQ=255)
                                                                                                                              ataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATAtt                     <  2:1054873/149‑1 (MQ=255)
                                                                                                                              ataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCACGCCAGCAGTTTATCCAGCTCATCGCCACAGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATAtt                     >  1:608900/1‑149 (MQ=255)
                                                                                                                               taACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTc                    <  1:1526578/149‑1 (MQ=255)
                                                                                                                               taACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTc                    >  1:544864/1‑149 (MQ=255)
                                                                                                                                aaCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCt                   <  1:1466165/149‑1 (MQ=255)
                                                                                                                                aaCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCt                   <  2:1355321/149‑1 (MQ=255)
                                                                                                                                 aCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTg                  <  1:358147/149‑1 (MQ=255)
                                                                                                                                 aCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTg                  <  2:1381021/149‑1 (MQ=255)
                                                                                                                                  cAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTgg                 <  1:1290322/149‑1 (MQ=255)
                                                                                                                                  cAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTgg                 >  2:990663/1‑149 (MQ=255)
                                                                                                                                     aTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCATTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAATCTCACAAATATTCTGGCAg              >  1:79844/1‑149 (MQ=255)
                                                                                                                                      tctcAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGc             >  2:533468/4‑149 (MQ=255)
                                                                                                                                       cccAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCg            >  2:684181/1‑149 (MQ=255)
                                                                                                                                        ccAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCgg           <  1:600576/149‑1 (MQ=255)
                                                                                                                                        ccAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCgg           >  2:1493177/1‑149 (MQ=255)
                                                                                                                                          aataatTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTTTGGCAGCGGTc         >  2:680056/1‑149 (MQ=255)
                                                                                                                                          aataatTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTc         >  2:702722/1‑149 (MQ=255)
                                                                                                                                          aataatTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTc         <  1:1507659/149‑1 (MQ=255)
                                                                                                                                            taatTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTAATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGc       <  1:643116/149‑1 (MQ=255)
                                                                                                                                            taatTAAGTTATTGGGATTTGGCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGc       <  1:972452/149‑1 (MQ=255)
                                                                                                                                               ttAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCaa    >  1:838943/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATA  >  NC_000913/3906683‑3906976

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