Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,440,173 T→C E286G (GAG→GGG)  rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,440,1730TC98.9% 305.7 / NA 89E286G (GAG→GGG) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (44/44);  minor base A (1/0);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGT  >  NC_000913_3_pae_pgi/3440026‑3440319
                                                                                                                                                   |                                                                                                                                                  
aaCCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCg                                                                                                                                                  <  1:438011/150‑1 (MQ=255)
   cTGTGATCCGGTTTCTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTa                                                                                                                                               <  4:603068/150‑1 (MQ=255)
    tgtgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTAc                                                                                                                                              >  2:97596/1‑150 (MQ=255)
    tgtgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTAc                                                                                                                                              >  2:267373/1‑150 (MQ=255)
    tgtgATCCGGTTACTCGTCAACGGTGCTTGCCGATGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGGCTTTAATCTCAGTAAGACATTTTCTACCAAGCTTAGCCGTTTTAAGGAGCTCAACCCCCGTAc                                                                                                                                              >  4:147100/1‑150 (MQ=255)
         tCCGGTTACTCGTGAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGTACGTCTGTGATCTCAGTAAGAGATGTTTTGCCAGGTTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGt                                                                                                                                         <  4:71226/150‑1 (MQ=255)
          ccGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAGGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTa                                                                                                                                        >  2:7892/1‑150 (MQ=255)
           cGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTAc                                                                                                                                       <  3:365594/150‑1 (MQ=255)
           cGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTAc                                                                                                                                       <  2:32744/150‑1 (MQ=255)
             gTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCACGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCa                                                                                                                                     <  1:301327/150‑1 (MQ=255)
                  tCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCa                                                                                                                                >  1:408578/1‑150 (MQ=255)
                       aGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGat                                                                                                                           >  4:713917/1‑150 (MQ=255)
                        gCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGata                                                                                                                          <  2:1086800/150‑1 (MQ=255)
                          gATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGAtata                                                                                                                        >  4:520760/1‑150 (MQ=255)
                               ttGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGa                                                                                                                   <  2:1028572/150‑1 (MQ=255)
                               ttGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGa                                                                                                                   >  2:297478/1‑150 (MQ=255)
                                      tGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTc                                                                                                            >  1:152612/1‑150 (MQ=255)
                                         ccAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                                                                                                         >  1:646985/1‑150 (MQ=255)
                                          cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                                                                                                        >  1:396168/1‑150 (MQ=255)
                                          cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                                                                                                        >  2:429031/1‑150 (MQ=255)
                                           aGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt                                                                                                       >  3:584192/1‑150 (MQ=255)
                                             ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa                                                                                                     <  3:501778/150‑1 (MQ=255)
                                             ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa                                                                                                     >  3:93853/1‑150 (MQ=255)
                                             ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa                                                                                                     >  1:133557/1‑150 (MQ=255)
                                              tttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa                                                                                                    <  2:274176/150‑1 (MQ=255)
                                               ttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGGACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg                                                                                                   >  1:637513/1‑150 (MQ=255)
                                                tCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg                                                                                                  <  3:155436/150‑1 (MQ=255)
                                                    ggCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                                                                                              >  1:343637/1‑150 (MQ=255)
                                                    ggCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                                                                                              <  1:429031/150‑1 (MQ=255)
                                                    ggCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                                                                                              >  1:428910/1‑150 (MQ=255)
                                                     gcgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt                                                                                             <  3:641874/150‑1 (MQ=255)
                                                     gcgcATGCCCAGAAACAGTCCACGGAAAGCCAGGACGGCGTTAAGCTCAGGGAGAGATGGTTTACCAAGGGTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt                                                                                             <  3:212177/150‑1 (MQ=255)
                                                       gcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGa                                                                                                                            <  2:378919/117‑1 (MQ=255)
                                                       gcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGa                                                                                                                            >  1:378919/1‑117 (MQ=255)
                                                        cATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc                                                                                          >  2:248045/1‑150 (MQ=255)
                                                         aTGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGCGCTCAACCACGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTCAGGCAGTTAGCa                                                                                         >  3:272541/1‑150 (MQ=255)
                                                         aTGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa                                                                                         >  3:272162/1‑150 (MQ=255)
                                                            cccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTg                                                                                                                     <  2:216921/119‑1 (MQ=255)
                                                            cccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTg                                                                                                                     >  1:216921/1‑119 (MQ=255)
                                                             ccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc                                                                                     <  3:147100/150‑1 (MQ=255)
                                                              cAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc                                                                                     <  1:592695/149‑1 (MQ=255)
                                                              cAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc                                                                                     >  2:592695/1‑149 (MQ=255)
                                                              cAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCg                                                                                    <  3:179294/150‑1 (MQ=255)
                                                                     aGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAggc                                                                             >  4:539441/1‑148 (MQ=255)
                                                                       tCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTTTTTTTTTTAGACGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATTGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCac                                                                           >  2:21018/1‑149 (MQ=255)
                                                                         cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc                                                                                          <  3:548212/133‑1 (MQ=255)
                                                                         cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc                                                                                          >  4:548212/1‑133 (MQ=255)
                                                                          aCGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTTTTTTTTTTTTTCGTCTTAAGGAGCTCAACCCCGGTACCCTGTATGAGATCACCGCTATAGTGGATATCTTCTCCTTTCAGGCAGTTAGCAGAGCGGCCATGCAATTc                                                                        >  2:1122242/1‑150 (MQ=255)
                                                                            gggAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCa                                                                      <  2:152612/150‑1 (MQ=255)
                                                                               aaGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATc                                                                  <  1:599312/150‑1 (MQ=255)
                                                                                  ccAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGt                                                                <  1:97596/150‑1 (MQ=255)
                                                                                    aGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCa                                                              >  4:638020/1‑150 (MQ=255)
                                                                                    aGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCa                                                              >  2:1057703/1‑150 (MQ=255)
                                                                                     gCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCaa                                                             <  3:460705/150‑1 (MQ=255)
                                                                                          tCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAACTCCAGATCGTCAACAggg                                                        >  3:275428/1‑150 (MQ=255)
                                                                                             ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGTTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGgcgc                                                     <  1:267373/150‑1 (MQ=255)
                                                                                                 tctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCgcagca                                                 >  2:267873/1‑150 (MQ=255)
                                                                                                 tctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCgcagca                                                 <  3:671574/150‑1 (MQ=255)
                                                                                                     aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggat                                             <  1:703024/150‑1 (MQ=255)
                                                                                                      gTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatc                                            <  3:570022/150‑1 (MQ=255)
                                                                                                      gTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatc                                            <  2:977780/150‑1 (MQ=255)
                                                                                                        aaGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgg                                          <  4:93853/150‑1 (MQ=255)
                                                                                                         agagATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgga                                         <  2:995119/150‑1 (MQ=255)
                                                                                                                 tttACCAAGGTTAGGCGGTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAActc                                 >  4:199090/1‑150 (MQ=255)
                                                                                                                  ttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAActct                                <  2:428910/150‑1 (MQ=255)
                                                                                                                     ccAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCCTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGt                             <  3:243266/150‑1 (MQ=255)
                                                                                                                      cAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGtt                            >  4:102539/1‑150 (MQ=255)
                                                                                                                      cAAGGTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGtt                            >  4:102728/1‑150 (MQ=255)
                                                                                                                         ggTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCttt                                                                                                      <  1:5610/73‑1 (MQ=255)
                                                                                                                         ggTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCttt                                                                                                      >  2:5610/1‑73 (MQ=255)
                                                                                                                         ggTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtct                         >  1:534601/1‑150 (MQ=255)
                                                                                                                         ggTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtct                         <  3:713917/150‑1 (MQ=255)
                                                                                                                         ggTTAGGCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtct                         <  1:7892/150‑1 (MQ=255)
                                                                                                                               gCGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTcttctt                   >  4:604987/1‑150 (MQ=255)
                                                                                                                                cGTTTTAAGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                  <  2:565460/150‑1 (MQ=255)
                                                                                                                                      aaGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa                                                                                                >  3:714837/1‑66 (MQ=255)
                                                                                                                                      aaGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa                                                                                                <  4:714837/66‑1 (MQ=255)
                                                                                                                                      aaGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTc            <  3:459905/150‑1 (MQ=255)
                                                                                                                                       aGGAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCa           <  2:1063502/150‑1 (MQ=255)
                                                                                                                                        ggAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAg          >  1:1089116/1‑150 (MQ=255)
                                                                                                                                         gAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAgg         >  2:115085/1‑150 (MQ=255)
                                                                                                                                         gAGCTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAgg         >  2:913598/1‑150 (MQ=255)
                                                                                                                                            cTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTg      >  3:503568/1‑150 (MQ=255)
                                                                                                                                            cTCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTg      >  4:703954/1‑150 (MQ=255)
                                                                                                                                             tCAACCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGa     <  1:267873/150‑1 (MQ=255)
                                                                                                                                               aaCCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTATCATAGCGGa                                                                                  <  4:65458/71‑1 (MQ=255)
                                                                                                                                               aaCCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa                                                                                  >  3:65458/1‑71 (MQ=255)
                                                                                                                                               aaCCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACg   <  2:343637/150‑1 (MQ=255)
                                                                                                                                                aCCCCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGt  <  1:537105/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGT  >  NC_000913_3_pae_pgi/3440026‑3440319

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