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% 582.8 / ‑4.7 163T285I (ACC→ATC) rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (83/79);  total (83/80)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 163 total aligned reads displayed.

TGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACA  >  NC_000913_3_pae_pgi/3440030‑3440322
                                                                                                                                                  |                                                                                                                                                  
tgtgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATAc                                                                                                                                                 >  2:1700117/1‑150 (MQ=255)
   gATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCt                                                                                                                                              <  2:117925/150‑1 (MQ=255)
     tCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGt                                                                                                                                            <  3:549151/150‑1 (MQ=255)
       cGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTAc                                                                                                                                          <  1:1640103/150‑1 (MQ=255)
        ggTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTAcc                                                                                                                                         >  1:1121754/1‑150 (MQ=255)
          ttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAg                                                                                                                                       <  3:591618/150‑1 (MQ=255)
          ttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAg                                                                                                                                       <  1:598361/150‑1 (MQ=255)
            aCTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGAt                                                                                                                                     <  2:166962/150‑1 (MQ=255)
            aCTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGAt                                                                                                                                     >  2:1493841/1‑150 (MQ=255)
               cGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCAc                                                                                                                                  >  1:266600/1‑150 (MQ=255)
                     cGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGatat                                                                                                                            >  1:836711/1‑150 (MQ=255)
                         gCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTCCCACGTTTGAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTg                                                                                                                        <  1:1746245/150‑1 (MQ=255)
                           ttGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGa                                                                                                                      >  4:1007568/1‑150 (MQ=255)
                             gCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATa                                                                                                                    >  1:202789/1‑150 (MQ=255)
                              ccGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTAc                                                                                                                                          <  4:936840/127‑1 (MQ=255)
                              ccGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTAc                                                                                                                                          >  3:936840/1‑127 (MQ=255)
                               cGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGc                                                                                                                  >  3:916701/1‑150 (MQ=255)
                                    gCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTg                                                                                                             <  4:687178/150‑1 (MQ=255)
                                      cAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                                                                                                           >  1:716211/1‑150 (MQ=255)
                                       aGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGa                                                                                                                      <  4:681734/138‑1 (MQ=255)
                                         ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATa                                                                                                                                                  <  2:1615140/108‑1 (MQ=255)
                                         ttttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa                                                                                                        <  4:267663/150‑1 (MQ=255)
                                           ttCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg                                                                                                      <  3:1007568/150‑1 (MQ=255)
                                              cAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa                                                                                                   >  1:913465/1‑150 (MQ=255)
                                              cAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa                                                                                                   <  1:602675/150‑1 (MQ=255)
                                                 gcgcATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt                                                                                                <  2:757003/150‑1 (MQ=255)
                                                     aTGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa                                                                                            >  1:839950/1‑150 (MQ=255)
                                                       gCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCaga                                                                                          >  1:279526/1‑150 (MQ=255)
                                                        cccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCagag                                                                                         >  2:248902/1‑150 (MQ=255)
                                                         ccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAATAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGa                                                                                                                                      >  1:268096/1‑104 (MQ=255)
                                                         ccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGTTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc                                                                                        >  1:1149687/1‑150 (MQ=255)
                                                         ccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTAc                                                                                                                                          <  2:683568/100‑1 (MQ=255)
                                                         ccAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc                                                                                        >  1:1275455/1‑150 (MQ=255)
                                                          cAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCg                                                                                       <  4:821085/150‑1 (MQ=255)
                                                           agagACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCt                                                                                                                                              >  2:1216948/1‑94 (MQ=255)
                                                             agaCAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGtt                                                                                                <  3:675659/138‑1 (MQ=255)
                                                                cAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTAGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGt                                                                                 <  2:434866/150‑1 (MQ=255)
                                                                   tCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCaa                                                                              >  1:1484861/1‑150 (MQ=255)
                                                                   tCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCaa                                                                              <  2:938347/150‑1 (MQ=255)
                                                                    ccACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAAt                                                                             <  3:889339/150‑1 (MQ=255)
                                                                       cGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTcc                                                                          >  3:712173/1‑150 (MQ=255)
                                                                       cGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTcc                                                                          >  3:712492/1‑150 (MQ=255)
                                                                        gggAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCa                                                                         <  1:719981/150‑1 (MQ=255)
                                                                        gggAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCa                                                                         >  3:159390/1‑150 (MQ=255)
                                                                          gAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGa                                                                       >  1:807671/1‑150 (MQ=255)
                                                                            aGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATc                                                                     <  2:807671/150‑1 (MQ=255)
                                                                            aGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATc                                                                     <  3:403202/150‑1 (MQ=255)
                                                                            aGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATc                                                                     <  1:1612941/150‑1 (MQ=255)
                                                                              ccAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGt                                                                   <  2:839950/150‑1 (MQ=255)
                                                                                  cACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCa                                                                                            >  4:630783/1‑121 (MQ=255)
                                                                                      tCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAggg                                                           <  1:227154/150‑1 (MQ=255)
                                                                                      tCTTTAATCTCAGCAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAggg                                                           <  1:1493841/150‑1 (MQ=255)
                                                                                       cTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa                                                                                                        >  2:252510/1‑104 (MQ=255)
                                                                                        tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTAc                                                                                                                                          <  1:619680/69‑1 (MQ=255)
                                                                                        tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTAc                                                                                                                                          >  2:619680/1‑69 (MQ=255)
                                                                                        tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACa                                                                                   >  3:206049/1‑124 (MQ=255)
                                                                                        tttAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgcg                                                         >  3:574881/1‑150 (MQ=255)
                                                                                           aaTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCgcag                                                      <  2:35766/150‑1 (MQ=255)
                                                                                           aaTCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCgcag                                                      <  1:1014308/150‑1 (MQ=255)
                                                                                             tctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg                                                                                                  >  1:1055355/1‑104 (MQ=255)
                                                                                             tctcAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg                                                                                                  <  2:1055355/104‑1 (MQ=255)
                                                                                                  gTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatc                                               <  1:1700117/150‑1 (MQ=255)
                                                                                                     agagATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCa                                                                                                                                        <  1:536708/58‑1 (MQ=255)
                                                                                                     agagATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGTTCGTCAACAGGGCGCAGCAggatcgga                                            <  3:788773/150‑1 (MQ=255)
                                                                                                      gagaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa                                                                                                       <  4:876425/90‑1 (MQ=255)
                                                                                                      gagaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggat                                           <  3:755006/150‑1 (MQ=255)
                                                                                                       agaTTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCATCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTTTGCTTTAAGGCAGTTAGCAGAGCGGAGAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                          <  4:915112/150‑1 (MQ=255)
                                                                                                            ttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgc                                                          <  3:111134/129‑1 (MQ=255)
                                                                                                            ttttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGgc                                                          >  4:111134/1‑129 (MQ=255)
                                                                                                             tttACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAActc                                    <  2:1275455/150‑1 (MQ=255)
                                                                                                                 ccAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCt                                                                                                                 <  1:1672972/69‑1 (MQ=255)
                                                                                                                  cAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCAc                                                                                                                                  <  2:542286/51‑1 (MQ=255)
                                                                                                                  cAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCAc                                                                                                                                  >  1:542286/1‑51 (MQ=255)
                                                                                                                  cAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGtt                               >  2:1782635/1‑150 (MQ=255)
                                                                                                                  cAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGtt                               >  3:591590/1‑150 (MQ=255)
                                                                                                                   aaGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAggatcggatc                                          <  3:347993/138‑1 (MQ=255)
                                                                                                                    aGGTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtc                             <  4:354345/150‑1 (MQ=255)
                                                                                                                      gTTAGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTtctc                           <  3:853349/150‑1 (MQ=255)
                                                                                                                         aGGCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttc                        >  2:604191/1‑150 (MQ=255)
                                                                                                                           gCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCa                                                                         <  4:725599/99‑1 (MQ=255)
                                                                                                                           gCGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTcttctt                      >  1:362497/1‑150 (MQ=255)
                                                                                                                            cGTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCttt                     <  2:1121754/150‑1 (MQ=255)
                                                                                                                             gTTTTAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTc                    >  2:1723237/1‑150 (MQ=255)
                                                                                                                              ttttAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCa                   <  4:712173/150‑1 (MQ=255)
                                                                                                                              ttttAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCa                   <  3:166294/150‑1 (MQ=255)
                                                                                                                                ttAAGGAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTAACt                 >  1:1359371/1‑150 (MQ=255)
                                                                                                                                    ggAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAg             >  2:1526373/1‑150 (MQ=255)
                                                                                                                                     gAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAgg            >  1:1177223/1‑150 (MQ=255)
                                                                                                                                     gAGCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAgg            <  1:1578084/150‑1 (MQ=255)
                                                                                                                                       gCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc                                                                                        <  2:1274553/72‑1 (MQ=255)
                                                                                                                                       gCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGc                                                                                        >  1:1274553/1‑72 (MQ=255)
                                                                                                                                       gCTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCt          >  2:652177/1‑150 (MQ=255)
                                                                                                                                        cTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCtcttct                       <  2:1417251/136‑1 (MQ=255)
                                                                                                                                        cTCAACCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTg         >  3:278002/1‑150 (MQ=255)
                                                                                                                                           aaCCTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACg      <  4:159390/150‑1 (MQ=255)
                                                                                                                                              cTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAc   >  4:207464/1‑150 (MQ=255)
                                                                                                                                              cTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAc   >  4:729451/1‑150 (MQ=255)
                                                                                                                                              cTCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTAc   >  2:269439/1‑150 (MQ=255)
                                                                                                                                               tCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                                                                                                           >  4:585896/1‑45 (MQ=255)
                                                                                                                                               tCGATACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACa  >  2:1559551/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGCAGAGCGGACAGTCAATTCCAGATCGTCAACAGGGCGCAGCAGGATCGGATCGAACTCTGGTTTCTCTTCTTTCACTTCAGGCTGACGTACA  >  NC_000913_3_pae_pgi/3440030‑3440322

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