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.4% 594.1 / NA 168T285I (ACC→ATC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (81/86);  total (81/87)
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

Only 100 of 168 total aligned reads displayed.

TGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGT  >  NC_000913_3_pae_pgi/3440030‑3440319
                                                                                                                                                  |                                                                                                                                               
tgtgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATAc                                                                                                                                              >  1:809344/1‑150 (MQ=255)
    aTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCACTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTg                                                                                                                                          <  3:851749/150‑1 (MQ=255)
       cGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTAc                                                                                                                                       <  1:911311/150‑1 (MQ=255)
         gTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCa                                                                                                                                     <  2:1542396/150‑1 (MQ=255)
                   aGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGat                                                                                                                           >  3:790500/1‑150 (MQ=255)
                    gCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGata                                                                                                                          <  1:1555342/150‑1 (MQ=255)
                         gCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTg                                                                                                                     <  3:404972/150‑1 (MQ=255)
                           ttGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGCGAAGTCAGCACGTCTGTAATCTCAGTAATAGATTTGTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGa                                                                                                                   <  3:361599/150‑1 (MQ=255)
                                ggTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCt                                                                                                              <  2:871026/150‑1 (MQ=255)
                                  gggCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTc                                                                                                            <  2:1012933/149‑1 (MQ=255)
                                     ccAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                                                                                                         >  3:340830/1‑150 (MQ=255)
                                     ccAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                                                                                                         <  4:409680/150‑1 (MQ=255)
                                     ccAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                                                                                                         <  2:282252/150‑1 (MQ=255)
                                        gTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCAcc                                                                                                                              >  4:113325/1‑126 (MQ=255)
                                         ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGAt                                                                                                                                  <  2:524559/121‑1 (MQ=255)
                                         ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGAt                                                                                                                                  >  1:524559/1‑121 (MQ=255)
                                         ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATc                                                                                                                                 <  2:722007/122‑1 (MQ=255)
                                         ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa                                                                                                     >  2:1015107/1‑150 (MQ=255)
                                          tttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa                                                                                                    <  3:102238/150‑1 (MQ=255)
                                          tttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa                                                                                                    <  1:1118722/150‑1 (MQ=255)
                                              cAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa                                                                                                <  1:1592667/150‑1 (MQ=255)
                                              cAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa                                                                                                <  3:13387/150‑1 (MQ=255)
                                              cAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTAAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCt                                                                                                           <  3:545819/139‑1 (MQ=255)
                                               aGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg                                                                                               >  4:377115/1‑150 (MQ=255)
                                                 gcgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt                                                                                             <  3:377115/150‑1 (MQ=255)
                                                 gcgcACGCCCAGAGACATTCCACGGGAAGCCAGCACGTGTTTAATCACAGAAAGAGCAGTTATACCAAGATAAGTCGTTTTAAGCAACTTAACCTCGATACGCCGAACCAGATCACCGAAAAAGGGCATAGCTTCTGCGTGAAGGCAGtt                                                                                             >  2:1222413/1‑150 (MQ=255)
                                                   gcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAg                                                                                           <  4:303360/150‑1 (MQ=255)
                                                    cATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAg                                                                                                                       <  2:1607382/121‑1 (MQ=255)
                                                    cATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAg                                                                                                                       >  1:1607382/1‑121 (MQ=255)
                                                     aTGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATa                                                                                                                                               >  3:623159/1‑96 (MQ=255)
                                                       gCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCaga                                                                                       >  3:599867/1‑150 (MQ=255)
                                                       gCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCaga                                                                                       >  1:1256358/1‑150 (MQ=255)
                                                         ccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTg                                                                                                          <  4:1015874/129‑1 (MQ=255)
                                                           agagACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCgg                                                                                   >  2:823628/1‑150 (MQ=255)
                                                            gagaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGa                                                                                  <  2:354123/150‑1 (MQ=255)
                                                             agaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGAc                                                                                 >  3:115706/1‑150 (MQ=255)
                                                              gaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGATAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                >  2:335999/1‑150 (MQ=255)
                                                               aCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGTTAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAg                                                                               <  3:506081/150‑1 (MQ=255)
                                                                     cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCa                                                                                                                                >  2:1595812/1‑95 (MQ=255)
                                                                     cACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAAtt                                                                         >  1:1187509/1‑150 (MQ=255)
                                                                      aCGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                >  2:991709/1‑142 (MQ=255)
                                                                      aCGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                <  1:991709/142‑1 (MQ=255)
                                                                      aCGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTc                                                                        <  2:254177/150‑1 (MQ=255)
                                                                      aCGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCAACGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTc                                                                        >  2:1506946/1‑150 (MQ=255)
                                                                       cGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTcc                                                                       <  1:732719/150‑1 (MQ=255)
                                                                        gggAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCa                                                                      <  1:1224214/150‑1 (MQ=255)
                                                                        cactgaccATCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTGAAGCGTTTTGAGGAGCTCAGCCGCGGTACGCTGTAGCAGATGAGCGATGTAGTGGATGGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGGCAATTCCa                                                                      <  4:881905/144‑1 (MQ=255)
                                                                        aggAAGCCAGCATGTCTGTAATCTCAGGAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCa                                                                      <  4:115706/149‑1 (MQ=255)
                                                                          gAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCACTTCCCGa                                                                    <  2:1358779/150‑1 (MQ=255)
                                                                           aaGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGAt                                                                   <  1:1479934/150‑1 (MQ=255)
                                                                            aGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATc                                                                  <  1:1015107/150‑1 (MQ=255)
                                                                              ccAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGt                                                                <  2:1256793/150‑1 (MQ=255)
                                                                                   aCGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACa                                                           <  1:43457/150‑1 (MQ=255)
                                                                                        tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTc                                                               <  4:709708/141‑1 (MQ=255)
                                                                                        tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTc                                                               >  3:709708/1‑141 (MQ=255)
                                                                                        tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgcg                                                      >  2:615275/1‑150 (MQ=255)
                                                                                         ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCAcc                                                                                                                              >  3:32683/1‑77 (MQ=255)
                                                                                         ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTc                                                                                                            >  1:691512/1‑95 (MQ=255)
                                                                                         ttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTc                                                                                                            <  2:691512/95‑1 (MQ=255)
                                                                                           aaTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgcg                                                      <  2:354868/147‑1 (MQ=255)
                                                                                           aaTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgcg                                                      >  1:354868/1‑147 (MQ=255)
                                                                                           aaTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGAGCACCGTTACCGGTGTGGGCAACGGCTGGTATGCATGCAGCAGGTCGTACACGCAATTCCAGATCGTCAACAGGGCgcag                                                   >  1:1426866/1‑150 (MQ=255)
                                                                                             tctcAGTAAGAGATTTTTTACTATGTTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCgcagca                                                 >  4:659576/1‑150 (MQ=255)
                                                                                             tctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGt                                                                                                                                         >  1:1347655/1‑62 (MQ=255)
                                                                                                 aGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggat                                             <  2:607695/150‑1 (MQ=255)
                                                                                                   tAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcg                                           >  1:283330/1‑150 (MQ=255)
                                                                                                     tgagaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg                                                                                                   >  2:646322/2‑92 (MQ=255)
                                                                                                     tgagaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg                                                                                                   <  1:646322/91‑1 (MQ=255)
                                                                                                       agaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                       <  1:490109/150‑1 (MQ=255)
                                                                                                         aTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGgatcggatcgg                                     >  4:633992/1‑149 (MQ=255)
                                                                                                          ttttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGaa                                    <  2:283330/150‑1 (MQ=255)
                                                                                                           tttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAAc                                   >  3:383487/1‑150 (MQ=255)
                                                                                                           tttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAAc                                   >  1:482295/1‑150 (MQ=255)
                                                                                                            ttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                >  3:775449/1‑104 (MQ=255)
                                                                                                            ttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                <  4:775449/104‑1 (MQ=255)
                                                                                                                aCCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTgg                              >  3:668980/1‑150 (MQ=255)
                                                                                                                 ccAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCCGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                >  2:1410880/1‑99 (MQ=255)
                                                                                                                 ccAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCCGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                <  1:1410880/99‑1 (MQ=255)
                                                                                                                  cAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAgg                                                         >  4:164944/1‑121 (MQ=255)
                                                                                                                  cAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGtt                            >  2:1202846/1‑150 (MQ=255)
                                                                                                                       ttAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGt                                                                                                                      >  3:930651/1‑55 (MQ=255)
                                                                                                                         aGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGt                                                                                              >  4:868927/1‑77 (MQ=255)
                                                                                                                          ggCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtctt                      >  1:564022/1‑148 (MQ=255)
                                                                                                                           gCGTTTTAAGGAGCTCAACCTCGATACGCTGTACTAGATCACCGATATAGTGGATAGc                                                                                                               <  2:409691/58‑1 (MQ=255)
                                                                                                                           gCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTcttctt                   >  1:331607/1‑150 (MQ=255)
                                                                                                                               tttAAGGAGCTCAATCTCGATTACCTGTAACAGAGTACCTATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCAc               <  1:907018/150‑1 (MQ=255)
                                                                                                                                   aGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCa                <  1:676011/145‑1 (MQ=255)
                                                                                                                                   aGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCa                >  2:676011/1‑145 (MQ=255)
                                                                                                                                   aGGAGCTCAACCTCGATACGCTGTACCAGATCACCGAAATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCa           >  2:1187383/1‑150 (MQ=255)
                                                                                                                                    ggAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAg          >  4:798610/1‑150 (MQ=255)
                                                                                                                                     gAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGa                                                                                                                   <  4:1033611/44‑1 (MQ=255)
                                                                                                                                     gAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAgg         <  3:1034612/150‑1 (MQ=255)
                                                                                                                                        cTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATc                                                                  >  4:476477/1‑90 (MQ=255)
                                                                                                                                        cTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATc                                                                  <  3:476477/90‑1 (MQ=255)
                                                                                                                                        cTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTg      >  2:1304425/1‑150 (MQ=255)
                                                                                                                                        cTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTg      >  4:599995/1‑150 (MQ=255)
                                                                                                                                           aaCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACg   >  2:1251695/1‑150 (MQ=255)
                                                                                                                                           aaCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACg   >  3:213782/1‑150 (MQ=255)
                                                                                                                                            aCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACgc  >  4:320027/1‑149 (MQ=255)
                                                                                                                                             ccTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCa                                                                            <  3:1027426/75‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGT  >  NC_000913_3_pae_pgi/3440030‑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: