Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,440,176 G→A T285I (ACC→ATC)  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,1760GA99.5% 736.2 / NA 203T285I (ACC→ATC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (92/109);  minor base C (1/0);  total (93/110)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

TGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAC  >  NC_000913_3_pae_pgi/3440030‑3440321
                                                                                                                                                  |                                                                                                                                                 
   gATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCt                                                                                                                                             <  4:762345/150‑1 (MQ=255)
     tCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGt                                                                                                                                           <  4:1362090/150‑1 (MQ=255)
     tCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGt                                                                                                                                           <  3:1583/150‑1 (MQ=255)
       cGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTAc                                                                                                                                         <  4:378371/150‑1 (MQ=255)
        ggTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCCGAGCCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCCCGCTTAGGCGATCTCAGGCGCTCAACCTCGCTACGCCGTAcc                                                                                                                                        >  2:1838565/1‑150 (MQ=255)
          ttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAg                                                                                                                                      >  2:1505702/1‑150 (MQ=255)
          ttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAg                                                                                                                                      <  2:1027764/150‑1 (MQ=255)
            aCTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGAt                                                                                                                                    >  3:1093906/1‑150 (MQ=255)
                gTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCt                                                                                                                                             >  4:52550/1‑137 (MQ=255)
                gTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCt                                                                                                                                             <  3:52550/137‑1 (MQ=255)
                gTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCAcc                                                                                                                                <  3:712435/150‑1 (MQ=255)
                gTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCAACGGGCAGCCAGCCCGTATTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGAGTTTTAAGGAGCTCAACCTCGATACGCTGTACAAGATCAcc                                                                                                                                >  2:1906216/1‑150 (MQ=255)
                 tCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACTTCTTTAATCTCAGTAAGAGATGTTTTACCAAGGGTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCg                                                                                                                               >  3:665596/1‑150 (MQ=255)
                 tCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCg                                                                                                                               <  3:967414/150‑1 (MQ=255)
                 tCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAGTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATAGGCTGTACCAGATCACCg                                                                                                                               >  3:491633/1‑150 (MQ=255)
                   aGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGt                                                                                                                                           >  1:1028725/1‑136 (MQ=255)
                       aTGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCGGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAg                                                                                                                         <  3:68757/150‑1 (MQ=255)
                         gCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTg                                                                                                                       <  2:734087/150‑1 (MQ=255)
                         gCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTg                                                                                                                       >  4:680907/1‑150 (MQ=255)
                              ccGGTGGCCATTTTCCCAGTTGTATGCCCAGATACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCTATGTTAGGCTTTTTAAGTAGCTCAACCTCTTTACGCTGTACCAGATCACCTATATAGTGGATAg                                                                                                                  <  4:1106942/150‑1 (MQ=255)
                               cGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGc                                                                                                                 >  2:2217819/1‑150 (MQ=255)
                                 gTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCtt                                                                                                               >  3:375947/1‑150 (MQ=255)
                                 gTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCtt                                                                                                               >  4:500085/1‑150 (MQ=255)
                                      cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGa                                                                                                                                     <  1:128012/123‑1 (MQ=255)
                                      cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGa                                                                                                                                     >  2:128012/1‑123 (MQ=255)
                                      cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                                                                                                          >  1:780299/1‑150 (MQ=255)
                                          tttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa                                                                                                      <  2:691133/150‑1 (MQ=255)
                                           ttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg                                                                                                     <  1:1505702/150‑1 (MQ=255)
                                           ttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg                                                                                                     <  1:1009729/150‑1 (MQ=255)
                                              cAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa                                                                                                  <  2:1593667/150‑1 (MQ=255)
                                              cAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa                                                                                                  >  1:1579844/1‑150 (MQ=255)
                                               aGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg                                                                                                     >  2:293131/1‑146 (MQ=255)
                                               aGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg                                                                                                 >  4:933446/1‑150 (MQ=255)
                                                 gcgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt                                                                                               <  1:585398/150‑1 (MQ=255)
                                                     aTGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCGAGGTTAGGCGTTTTAAGGAGCTCAGCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa                                                                                           <  1:1838565/150‑1 (MQ=255)
                                                        cccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag                                                                                        >  1:309682/1‑150 (MQ=255)
                                                        cccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGGAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGATCAGCCTCGATACGCTGTACCAGATCGCCGCTATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag                                                                                        >  2:1733410/1‑150 (MQ=255)
                                                             agaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGAc                                                                                   <  1:2217819/150‑1 (MQ=255)
                                                                cAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGt                                                                                <  1:1538349/150‑1 (MQ=255)
                                                                 aGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTc                                                                               >  3:243180/1‑150 (MQ=255)
                                                                   tCCATGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCaa                                                                             >  3:1357254/1‑150 (MQ=255)
                                                                   tCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGt                                                                                <  2:1583127/147‑1 (MQ=255)
                                                                   tCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCaa                                                                             >  2:1918297/1‑150 (MQ=255)
                                                                     cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCATCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCttt                                                                                                        <  2:1396265/121‑1 (MQ=255)
                                                                     cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCttt                                                                                                        >  1:1396265/1‑121 (MQ=255)
                                                                     cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa                                                                                                       <  2:1995670/122‑1 (MQ=255)
                                                                     cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAAtt                                                                           >  2:516769/1‑150 (MQ=255)
                                                                     cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGAGACGCTGTACCAGATCACCGATATAGTGTATAGCTTCTGCTTTAAGGCAGTTGGCAGAGCGGACAGTCAAtt                                                                           <  2:1147758/150‑1 (MQ=255)
                                                                     cACGGGAAGCCAGCACGACTTTACTCTCAGTAAGAGATTTTTTACCAAGCATAGGCGTTCTAAGGAACTCAACCTCGACACGCTGTACCAGATCACCGATATAGTGCATACCTTCTCCTCTAAGGCACTTAGCAGAGCGGACAGTCAAtc                                                                           >  2:848810/1‑149 (MQ=255)
                                                                      aCGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTc                                                                          >  1:712856/1‑150 (MQ=255)
                                                                        gggAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCa                                                                        <  3:500085/150‑1 (MQ=255)
                                                                           aaGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGAt                                                                     <  1:2100127/150‑1 (MQ=255)
                                                                                aGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCg                                                                   <  1:385310/147‑1 (MQ=255)
                                                                                 gCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTc                                                                                                              >  1:2126594/1‑103 (MQ=255)
                                                                                  cACGTCTTTAATCTCAGTATGAGATTGTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCTCCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGTACAGTCAATTCCAGATCGTCAAc                                                              <  2:581215/150‑1 (MQ=255)
                                                                                   aCGTCTTTAATCTCAGTCAGAGATTTTTTTGCATGGGTAGGCGTTTTAAGGAGCTCAGCCTCGAAACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAAct                                                             <  2:2127433/150‑2 (MQ=255)
                                                                                   aCGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACa                                                             <  1:1918297/150‑1 (MQ=255)
                                                                                     gTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGa                                                                                                                                     <  4:261268/76‑1 (MQ=255)
                                                                                     gTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg                                                                                                     >  1:1798563/1‑108 (MQ=255)
                                                                                     gTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg                                                                                                     <  2:1798563/108‑1 (MQ=255)
                                                                                     gTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAgg                                                           >  1:2297647/1‑150 (MQ=255)
                                                                                        tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAgg                                                                                                    >  2:2040969/1‑106 (MQ=255)
                                                                                            aTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCgcagc                                                    >  1:941847/1‑150 (MQ=255)
                                                                                              ctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                                                                                                           >  2:2104861/1‑93 (MQ=255)
                                                                                                  gTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatc                                              >  1:2030744/1‑150 (MQ=255)
                                                                                                     agagATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcgga                                           <  3:84667/150‑1 (MQ=255)
                                                                                                       agaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                         >  1:2012620/1‑150 (MQ=255)
                                                                                                        gaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTc                                                                                                              <  3:440836/80‑1 (MQ=255)
                                                                                                        gaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTc                                                                                                              >  4:440836/1‑80 (MQ=255)
                                                                                                         aTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGa                                       >  2:1855674/1‑150 (MQ=255)
                                                                                                                 ccAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc                                                                                       >  1:825195/1‑94 (MQ=255)
                                                                                                                  cAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa                                                                                           >  1:1612048/1‑89 (MQ=255)
                                                                                                                  cAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGtt                              >  1:1585991/1‑150 (MQ=255)
                                                                                                                      gTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtctc                          <  2:2389039/150‑1 (MQ=255)
                                                                                                                        tAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtctt                        <  1:2108731/150‑1 (MQ=255)
                                                                                                                          ggCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                      >  2:2123228/1‑150 (MQ=255)
                                                                                                                           gCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAg                                                                                                                         >  1:142764/1‑50 (MQ=255)
                                                                                                                           gCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa                                                                                                      >  2:2053164/1‑69 (MQ=255)
                                                                                                                           gCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa                                                                                                      <  1:2053164/69‑1 (MQ=255)
                                                                                                                           gCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                  <  3:327608/89‑1 (MQ=255)
                                                                                                                             gTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTc                   <  1:1961289/150‑1 (MQ=255)
                                                                                                                              ttttAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                  >  3:48654/1‑86 (MQ=255)
                                                                                                                              ttttAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                  <  4:48654/86‑1 (MQ=255)
                                                                                                                              ttttAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCa                  <  2:802673/150‑1 (MQ=255)
                                                                                                                              ttttAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCGGGTTTCTCTTCTTTCa                  >  3:1307991/1‑150 (MQ=255)
                                                                                                                                ttAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACt                <  2:2297647/150‑1 (MQ=255)
                                                                                                                                   aGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatc                                              <  2:1343245/117‑1 (MQ=255)
                                                                                                                                   aGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatc                                              >  1:1343245/1‑117 (MQ=255)
                                                                                                                                       gCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCg                                                                   <  2:146648/92‑1 (MQ=255)
                                                                                                                                        cgcAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGTCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTg        <  1:1906216/148‑1 (MQ=255)
                                                                                                                                        cTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                    <  1:1161514/138‑1 (MQ=255)
                                                                                                                                          cAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTa                                                                                              <  1:1563278/62‑1 (MQ=255)
                                                                                                                                          cAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTa                                                                                              >  2:1563278/1‑62 (MQ=255)
                                                                                                                                          cAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttc                       <  4:1202206/133‑1 (MQ=255)
                                                                                                                                          cAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttc                       >  3:1202206/1‑133 (MQ=255)
                                                                                                                                           aaCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTa                                                                                              <  1:2173640/61‑1 (MQ=255)
                                                                                                                                           aaCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTa                                                                                              >  2:2173640/1‑61 (MQ=255)
                                                                                                                                           aaCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACg     >  4:153220/1‑150 (MQ=255)
                                                                                                                                           aaCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACg     >  2:1212663/1‑150 (MQ=255)
                                                                                                                                              cTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAc  <  1:1212663/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAC  >  NC_000913_3_pae_pgi/3440030‑3440321

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: