Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 1,623,718 C→A D45Y (GAT→TAT)  dgcZ ← diguanylate cyclase, zinc‑sensing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,623,7180CA95.5% 476.3 / 10.3 134D45Y (GAT→TAT) dgcZdiguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (62/66);  total (65/69)
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.66e-01

Only 100 of 134 total aligned reads displayed.

CCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCACTCCTTTT  >  NC_000913_3_pae_pgi/1623577‑1623863
                                                                                                                                             |                                                                                                                                                 
 cAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTTATTTc                                                                                                                                          >  4:937986/1‑150 (MQ=255)
  aGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCt                                                                                                                                         <  2:93749/150‑1 (MQ=255)
    cATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTgg                                                                                                                                       <  3:1053841/150‑1 (MQ=255)
    cATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTgg                                                                                                                                       <  4:136253/150‑1 (MQ=255)
    cATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTgg                                                                                                                                       <  2:745702/150‑1 (MQ=255)
     aTTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGc                                                                                                                                      <  2:447605/150‑1 (MQ=255)
       tAATTCCCGACCACAGTTATGCATATGTTGATGGGCGGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCtt                                                                                                                                    >  1:1020504/1‑150 (MQ=255)
          ttCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACt                                                                                                                                 <  3:762998/150‑1 (MQ=255)
           tCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTg                                                                                                                                >  1:1367101/1‑150 (MQ=255)
                aCCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATCCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTTATTTCTGGCTTACTGGCAtt                                                                                                                           >  3:863290/1‑149 (MQ=255)
                     aGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTcg                                                                                                                      <  3:408658/150‑1 (MQ=255)
                      gTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTcgc                                                                                                                     <  1:1072908/150‑1 (MQ=255)
                      gTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTcgc                                                                                                                     <  1:868877/150‑1 (MQ=255)
                       ttATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTcgcg                                                                                                                    <  2:496159/150‑1 (MQ=255)
                         aTGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGAc                                                                                                                  <  1:1026226/150‑1 (MQ=255)
                             atatGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACg                                                                                                              >  3:683544/1‑150 (MQ=255)
                               atGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAGCTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCt                                                                                                            <  4:583203/150‑1 (MQ=255)
                                tGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtg                                                                                                           <  4:354458/150‑1 (MQ=255)
                                 gTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtgt                                                                                                          >  4:846145/1‑150 (MQ=255)
                                  ttGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtgtg                                                                                                         >  2:1026395/1‑150 (MQ=255)
                                   tGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGc                                                                                                                              <  3:552312/128‑1 (MQ=255)
                                   tGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGc                                                                                                                              >  4:552312/1‑128 (MQ=255)
                                   tGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGa                                                                                                        <  1:965281/150‑1 (MQ=255)
                                    gATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGaa                                                                                                       <  1:754371/150‑1 (MQ=255)
                                     aTGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGaa                                                                                                       <  4:389527/149‑1 (MQ=255)
                                     aTGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGaa                                                                                                       >  3:389527/1‑149 (MQ=255)
                                         gCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATa                                                                                                  >  2:1506395/1‑150 (MQ=255)
                                         gCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATa                                                                                                  <  1:1633742/150‑1 (MQ=255)
                                         gCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATa                                                                                                  >  1:331424/1‑150 (MQ=255)
                                          cAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATAc                                                                                                 >  1:1054309/1‑150 (MQ=255)
                                          cAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATAc                                                                                                 >  1:1054153/1‑150 (MQ=255)
                                          cAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATAc                                                                                                 <  4:17452/150‑1 (MQ=255)
                                                ccATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCa                                                                                           <  3:752661/150‑1 (MQ=255)
                                                 cATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAg                                                                                          >  3:259930/1‑150 (MQ=255)
                                                 cATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAg                                                                                          >  2:554243/1‑150 (MQ=255)
                                                  aTTAGCCGAACATAAGGTAATTCATTGTTATTGAGTAGCCTCAGATGATCATTTCACCGACCAAATTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTTGCGACCACGCTGTGAAACATACTCACCatc                                                                                         >  2:1120575/1‑148 (MQ=255)
                                                      gCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACt                                                                                     >  4:354929/1‑150 (MQ=255)
                                                      gCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACt                                                                                     <  3:278473/150‑1 (MQ=255)
                                                       ccGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTg                                                                                    <  2:929315/150‑1 (MQ=255)
                                                       ccGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGAcc                                                                                                                 <  2:800349/121‑1 (MQ=255)
                                                       ccGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGAcc                                                                                                                 >  1:800349/1‑121 (MQ=255)
                                                        cGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTgg                                                                                   >  1:328950/1‑150 (MQ=255)
                                                         gAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGt                                                                                  >  4:17142/1‑150 (MQ=255)
                                                               aaGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtgt                                                                                                          <  2:993722/120‑1 (MQ=255)
                                                               aaGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGACACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtgt                                                                                                          >  1:993722/1‑120 (MQ=255)
                                                                aGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGc                                                                           <  2:1367101/150‑1 (MQ=255)
                                                                 ggTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACt                                                                                                                                 >  1:157515/1‑95 (MQ=255)
                                                                 ggTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACt                                                                                                                                 <  2:157515/95‑1 (MQ=255)
                                                                  gTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCAt                                                                         >  1:457605/1‑150 (MQ=255)
                                                                  gTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCAt                                                                         <  4:683544/150‑1 (MQ=255)
                                                                  gTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCAt                                                                         >  2:313440/1‑150 (MQ=255)
                                                                  gTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCAt                                                                         >  3:560785/1‑150 (MQ=255)
                                                                    aaTTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCg                                                                       <  3:447010/150‑1 (MQ=255)
                                                                     aTTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGa                                                                      <  2:1178352/150‑1 (MQ=255)
                                                                      tacATCGTTATCGAGTGGCCCCAGATGATCCATCCACCGACCCAACTGGCACAGTCCATAAGAATGGTTCTACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGAt                                                                     <  2:839093/148‑1 (MQ=255)
                                                                       tCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATa                                                                    >  4:615658/1‑150 (MQ=255)
                                                                         aTCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGc                                                                  >  2:942846/1‑150 (MQ=255)
                                                                         aTCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGc                                                                  >  1:566261/1‑150 (MQ=255)
                                                                            gTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCtt                                                               >  3:201155/1‑150 (MQ=255)
                                                                             ttATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACg                                                                                                                                                 <  2:75283/67‑1 (MQ=255)
                                                                                     tGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACAt                                                                                                   >  4:289127/1‑105 (MQ=255)
                                                                                     tGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGAtt                                                      <  1:658189/150‑1 (MQ=255)
                                                                                       gCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTa                                                    >  1:56337/1‑150 (MQ=255)
                                                                                        ccccAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATATGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTAGGCATCGATAGCCTTATTGAGATTTaa                                                   >  1:392823/1‑150 (MQ=255)
                                                                                           cAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACaa                                                >  4:302376/1‑150 (MQ=255)
                                                                                             gatgatCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGa                                              <  2:1020504/150‑1 (MQ=255)
                                                                                              atgatCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGAt                                             >  4:829233/1‑150 (MQ=255)
                                                                                              atgatCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGAt                                             <  2:1566860/150‑1 (MQ=255)
                                                                                               tgttCAAGCCACCGACCATACTGGAACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGAATCGCTCGCGACCACGCTGTGAAACATACTCACAAGCGACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTGACAAGATg                                            <  4:713657/147‑1 (MQ=255)
                                                                                                     aTCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCt                                                                                                                                         <  3:296617/51‑1 (MQ=255)
                                                                                                     aTCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCt                                                                                                                                         >  4:296617/1‑51 (MQ=255)
                                                                                                     aTCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGCGATTTAACAAGATGGCATCa                                      >  4:1054868/1‑150 (MQ=255)
                                                                                                     aTCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCa                                      >  1:964723/1‑150 (MQ=255)
                                                                                                         aCCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTga                                                          >  3:46077/1‑126 (MQ=255)
                                                                                                         aCCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTga                                                          <  4:46077/126‑1 (MQ=255)
                                                                                                         aCCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAAttt                                  <  4:682295/150‑1 (MQ=255)
                                                                                                         aCCGACCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAAttt                                  >  2:701714/1‑150 (MQ=255)
                                                                                                           cGACCAAGCTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTcc                                <  2:1266133/150‑1 (MQ=255)
                                                                                                             aCCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATgg                                           >  2:425671/1‑137 (MQ=255)
                                                                                                             aCCAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGt                              >  3:294645/1‑150 (MQ=255)
                                                                                                              tcAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGtt                             <  4:1034660/149‑1 (MQ=255)
                                                                                                              ccAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGtt                             <  1:550910/150‑1 (MQ=255)
                                                                                                              ccAAACTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGtt                             >  1:620811/1‑150 (MQ=255)
                                                                                                                aaaCTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGt                           >  1:1664393/1‑150 (MQ=255)
                                                                                                                 aaCTGGCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGtc                          <  2:510448/150‑1 (MQ=255)
                                                                                                                     ggCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGtcttct                      <  1:1418873/150‑1 (MQ=255)
                                                                                                                     ggCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGtcttct                      >  2:77691/1‑150 (MQ=255)
                                                                                                                     ggCACAGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGtcttct                      <  1:1005443/150‑1 (MQ=255)
                                                                                                                       cacaGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTg                    <  2:331424/150‑1 (MQ=255)
                                                                                                                        acaGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAg                                                                                >  1:779463/1‑89 (MQ=255)
                                                                                                                        acaGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAg                                                                                <  2:779463/89‑1 (MQ=255)
                                                                                                                         caGTCCATAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGAt                  <  2:635511/150‑1 (MQ=255)
                                                                                                                               aTAAGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAAc                                                  >  1:941413/1‑112 (MQ=255)
                                                                                                                                 aaGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCa          <  2:52810/150‑1 (MQ=255)
                                                                                                                                  aGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTga                                                          >  3:947351/1‑101 (MQ=255)
                                                                                                                                  aGAATGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTga                                                          <  4:947351/101‑1 (MQ=255)
                                                                                                                                     aTGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCACTcc      <  2:432965/150‑1 (MQ=255)
                                                                                                                                     aTGGTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCACTcc      >  1:432965/1‑150 (MQ=255)
                                                                                                                                        gTTATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCACTCCttt   <  2:278125/150‑1 (MQ=255)
                                                                                                                                         ttATACGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCACTCCtttt  <  1:1468882/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CCAGCATTAATTCCCGACCACAGTTATGCATATGTTGATGGGCAGAATCCATTAGCCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCACTCCTTTT  >  NC_000913_3_pae_pgi/1623577‑1623863

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