Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 2,867,269 T→A I95F (ATC→TTC)  rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi2,867,2690TA98.8% 582.3 / ‑3.9 162I95F (ATC→TTC) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (84/75);  total (86/76)
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 163 total aligned reads displayed.

GAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAT  >  NC_000913_3_pae_pgi/2867124‑2867414
                                                                                                                                                 |                                                                                                                                                 
 agaAGCGGAATCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATcc                                                                                                                                              <  3:828052/150‑1 (MQ=255)
   aaGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCcgg                                                                                                                                            <  2:93602/150‑1 (MQ=255)
   aaGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCcgg                                                                                                                                            <  4:141286/150‑1 (MQ=255)
         aaaCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCga                                                                                                                                      >  4:488831/1‑150 (MQ=255)
             cACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGGGGGAAATTTTTAACAACAGACGCCCGCTACCCTCGAACATCCGGCGGCGAGAgg                                                                                                                                  >  1:808029/1‑150 (MQ=255)
             cACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAgg                                                                                                                                  >  1:932044/1‑150 (MQ=255)
             cACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAgg                                                                                                                                  <  2:1248344/150‑1 (MQ=255)
               cGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCg                                                                                                                                <  4:950050/150‑1 (MQ=255)
                     gggTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACAtct                                                                                                                          >  2:514061/1‑150 (MQ=255)
                      ggTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACAtctc                                                                                                                         >  1:1078835/1‑150 (MQ=255)
                      ggTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACAtctc                                                                                                                         <  4:755305/150‑1 (MQ=255)
                       gTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTcc                                                                                                                        >  4:590110/1‑150 (MQ=255)
                          aaaCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGATATCTCCACg                                                                                                                     <  3:847740/150‑1 (MQ=255)
                          aaaCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACg                                                                                                                     >  1:1340622/1‑150 (MQ=255)
                          aaaCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACg                                                                                                                     <  2:269003/150‑1 (MQ=255)
                           aaCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGc                                                                                                                    >  3:599075/1‑150 (MQ=255)
                            aCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCa                                                                                                                   >  4:799072/1‑150 (MQ=255)
                             cTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAg                                                                                                                  >  3:784451/1‑150 (MQ=255)
                              tactcaACCGCGCGGATCAGCCCCAGGTTGCCCACTACGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCACGGTACCCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGt                                                                                                                 >  2:977995/3‑150 (MQ=255)
                                   tACCGCGCTGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGcgcg                                                                                                            <  4:888721/150‑1 (MQ=255)
                                     ccGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGcgcgtc                                                                                                          <  3:71173/150‑3 (MQ=255)
                                       gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCa                                                                                                                   >  2:348687/1‑139 (MQ=255)
                                       gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGcgccgc                                                                                                        >  2:511810/1‑150 (MQ=255)
                                        cgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgcg                                                                                                       >  4:12556/1‑150 (MQ=255)
                                                     cAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtc                                                                                          <  2:765211/150‑1 (MQ=255)
                                                         ttGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttct                                                                                      <  2:1495463/150‑1 (MQ=255)
                                                          tGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttctt                                                                                     <  4:1073420/150‑1 (MQ=255)
                                                             ccTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCgg                                                                                  <  3:945360/150‑1 (MQ=255)
                                                             ccTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCgg                                                                                  >  3:396541/1‑150 (MQ=255)
                                                               tctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCg                                                                                   >  2:1716439/1‑147 (MQ=255)
                                                               tctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCg                                                                                   <  1:1716439/147‑1 (MQ=255)
                                                               tctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGcc                                                                                >  3:353474/1‑150 (MQ=255)
                                                               tctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGcc                                                                                >  1:1138930/1‑150 (MQ=255)
                                                                ctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCg                                                                               >  2:1444665/1‑150 (MQ=255)
                                                                 ttCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGt                                                                              <  1:1552069/150‑1 (MQ=255)
                                                                  tCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGtt                                                                             <  1:1469035/150‑1 (MQ=255)
                                                                    gATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTaa                                                                           >  3:730117/1‑150 (MQ=255)
                                                                      tAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGc                                                                                 <  1:297859/142‑1 (MQ=255)
                                                                      tAAGGTCCAGAAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCGCCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCGCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACa                                                                         <  3:871582/150‑1 (MQ=255)
                                                                          gTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtc                                                                                          <  2:684343/129‑1 (MQ=255)
                                                                           tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATaa                                                                                               <  2:548655/123‑1 (MQ=255)
                                                                           tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgtggt                                                                    <  3:646544/150‑1 (MQ=255)
                                                                           tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgtggt                                                                    >  2:195425/1‑150 (MQ=255)
                                                                              aGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGaa                                                                 <  3:758224/150‑1 (MQ=255)
                                                                              aGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGaa                                                                 >  2:108382/1‑150 (MQ=255)
                                                                              aGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGaa                                                                 >  1:1592664/1‑150 (MQ=255)
                                                                              aGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGaa                                                                 >  1:1513145/1‑150 (MQ=255)
                                                                               gCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCgg                                                                                  <  1:1181772/132‑1 (MQ=255)
                                                                               gCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAt                                                                <  2:1305122/150‑1 (MQ=255)
                                                                                cAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATa                                                               <  3:488831/150‑1 (MQ=255)
                                                                                   cGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAAcc                                                            <  2:1340622/150‑1 (MQ=255)
                                                                                      cAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAt                                                         >  4:1015366/1‑150 (MQ=255)
                                                                                       aGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAtc                                                        <  2:441851/150‑1 (MQ=255)
                                                                                       aGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAtc                                                        >  1:513075/1‑150 (MQ=255)
                                                                                         aCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAtctc                                                      >  1:124248/1‑150 (MQ=255)
                                                                                            aCGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCAcc                                                   >  1:701241/1‑150 (MQ=255)
                                                                                            aCGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCAcc                                                   >  3:976330/1‑150 (MQ=255)
                                                                                             cGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCa                                                  >  4:519600/1‑150 (MQ=255)
                                                                                              gATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCaa                                                 >  1:685686/1‑150 (MQ=255)
                                                                                               aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAg                                                <  1:1444665/150‑1 (MQ=255)
                                                                                                ttGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAgg                                               <  1:589871/150‑1 (MQ=255)
                                                                                                     aTAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAg                                          <  4:784451/150‑1 (MQ=255)
                                                                                                      tAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgc                                                                                                        >  3:1092071/1‑87 (MQ=255)
                                                                                                      tAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgc                                                                                                        <  4:1092071/87‑1 (MQ=255)
                                                                                                      tAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGc                                         >  2:798651/1‑150 (MQ=255)
                                                                                                       aaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCt                                        <  3:519600/150‑1 (MQ=255)
                                                                                                       aaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCt                                        >  3:21748/1‑150 (MQ=255)
                                                                                                        aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGc                                                                                                                                 >  2:973155/1‑60 (MQ=255)
                                                                                                        aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGc                                                                                                                                 <  1:973155/60‑1 (MQ=255)
                                                                                                          ggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAg                                     >  1:1006505/1‑150 (MQ=255)
                                                                                                            cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTc                                   <  1:904242/150‑1 (MQ=255)
                                                                                                             gggCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTcg                                  <  4:829203/150‑1 (MQ=255)
                                                                                                                 aaTTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTc                              <  3:494039/150‑1 (MQ=255)
                                                                                                                  aTTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTcc                             <  1:108382/150‑1 (MQ=255)
                                                                                                                   tttttACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCa                            >  2:1494412/1‑150 (MQ=255)
                                                                                                                    ttttACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCaa                           <  4:166339/150‑1 (MQ=255)
                                                                                                                         ccaccaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACa                                                                         <  1:1008858/99‑1 (MQ=255)
                                                                                                                         ccaccaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACa                                                                         >  2:1008858/1‑99 (MQ=255)
                                                                                                                          caccaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACg                     >  1:1524567/1‑150 (MQ=255)
                                                                                                                          caccaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACg                     >  4:674069/1‑150 (MQ=255)
                                                                                                                           accaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGa                                                                                                                               <  1:23270/43‑1 (MQ=255)
                                                                                                                           accaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGa                                                                                                                               >  2:23270/1‑43 (MQ=255)
                                                                                                                           accaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGc                    >  3:715061/1‑150 (MQ=255)
                                                                                                                           accaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGc                    <  4:353474/150‑1 (MQ=255)
                                                                                                                            ccaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                   <  1:305948/150‑1 (MQ=255)
                                                                                                                            ccaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                   >  1:691282/1‑150 (MQ=255)
                                                                                                                              aGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgtgg                                                                     <  4:908846/98‑1 (MQ=255)
                                                                                                                               gACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtg                <  2:1138930/150‑1 (MQ=255)
                                                                                                                                acccAAGTTACTCTCGAACATCCGCCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGGGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt               <  2:808029/147‑1 (MQ=255)
                                                                                                                                  gCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTgg             >  2:1644543/1‑150 (MQ=255)
                                                                                                                                    aaGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCt           >  2:1317279/1‑150 (MQ=255)
                                                                                                                                    aaGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTATTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCt           >  2:1144080/1‑150 (MQ=255)
                                                                                                                                      gTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCa                                                  >  2:665995/1‑109 (MQ=255)
                                                                                                                                      gTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCa                                                  <  1:665995/109‑1 (MQ=255)
                                                                                                                                         aCTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTg      <  3:225658/150‑1 (MQ=255)
                                                                                                                                          ctctCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGc     <  1:679163/150‑1 (MQ=255)
                                                                                                                                           tctcGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAg                                                <  1:465180/106‑1 (MQ=255)
                                                                                                                                             tcGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt  >  2:1596679/1‑150 (MQ=255)
                                                                                                                                             tcGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt  >  3:566601/1‑150 (MQ=255)
                                                                                                                                             tcGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCCCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt  >  2:1190195/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAT  >  NC_000913_3_pae_pgi/2867124‑2867414

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: