Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 2,867,340 A→C L71* (TTA→TGA)  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,3400AC95.4% 674.4 / 15.6 198L71* (TTA→TGA) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base A (6/3);  new base C (98/90);  total (104/94)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

Only 100 of 206 total aligned reads displayed.

TAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAAC  >  NC_000913_3_pae_pgi/2867194‑2867488
                                                                                                                                                  |                                                                                                                                                    
tAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTcaca                                                                                                                                                   <  1:848854/150‑1 (MQ=255)
tAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTcaca                                                                                                                                                   <  1:311336/150‑1 (MQ=255)
tAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGACAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTcaca                                                                                                                                                   <  2:865004/150‑1 (MQ=255)
 aaGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAg                                                                                                                                                  >  2:1554255/1‑150 (MQ=255)
   ggTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtg                                                                                                                                                >  1:1214755/1‑150 (MQ=255)
    gTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCACCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtgg                                                                                                                                               >  2:1611197/1‑150 (MQ=255)
     tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtggt                                                                                                                                              <  2:346330/150‑1 (MQ=255)
     tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtggt                                                                                                                                              <  3:366892/150‑1 (MQ=255)
      ccAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtggtg                                                                                                                                             <  4:208026/150‑1 (MQ=255)
        aGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTCCACAGTGGTGaa                                                                                                                                           >  3:427834/1‑150 (MQ=255)
          cAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATa                                                                                                                                         >  2:150583/1‑150 (MQ=255)
           aaCGCCAGCACACGATCGCCATAACGGCGGGCGATTTATCACACCATGCGAGAGTTTCTCTTGGTCATCCGGCGGCGAGAGGCGACATCTCCGCGCAGTGCGCGACGCGAAAAATAAACTTCTGCTGCGGCCGTTCACAGTGGTGAATaa                                                                                                                                        <  2:1543099/150‑1 (MQ=255)
           aaCGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATaa                                                                                                                                        >  3:458850/1‑150 (MQ=255)
               ccAGACCACGATTGCCATAACGGCGGGCAATGTTTTCCACCAGACGCAAGTTCCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCTCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCTTTCACAGTGGTGAATAACCaa                                                                                                                                    <  1:1510574/150‑1 (MQ=255)
                cAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGAGGCGAGCGACGAAATATACACGACGTGCGCGAAGCGCAAAATAAACCTCTTCTTCGGCCGTTCACAGTGGTGCATAACCaaa                                                                                                                                   >  3:782921/1‑149 (MQ=255)
                  gACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAAtct                                                                                                                                 >  2:979714/1‑150 (MQ=255)
                    ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCa                                                                                                                               >  3:762345/1‑150 (MQ=255)
                    ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCa                                                                                                                               <  1:931570/150‑1 (MQ=255)
                     cACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCAc                                                                                                                              >  3:541386/1‑150 (MQ=255)
                        gATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCaa                                                                                                                           <  3:226380/150‑1 (MQ=255)
                        gATTGCCATAACGGCGGGCAATTTTTACAACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGCTCACAGTGGTGAATAACCAATCTCACCaa                                                                                                                           >  2:523958/1‑150 (MQ=255)
                             ccATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTaa                                                                                                                      >  2:934006/1‑150 (MQ=255)
                              cATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCAcc                                                                                                                             >  4:71244/1‑142 (MQ=255)
                              cATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCAcc                                                                                                                             <  3:71244/142‑1 (MQ=255)
                               aTAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCAACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAg                                                                                                                    <  1:1116615/150‑1 (MQ=255)
                                  aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTg                                                                                                                 <  3:944393/150‑1 (MQ=255)
                                  aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTg                                                                                                                 <  3:887398/150‑1 (MQ=255)
                                   cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGa                                                                                                                <  2:704971/150‑1 (MQ=255)
                                    ggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAg                                                                                                               >  2:225119/1‑150 (MQ=255)
                                      cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTc                                                                                                             >  3:625821/1‑150 (MQ=255)
                                         gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTcgc                                                                                                           >  2:324660/1‑149 (MQ=255)
                                             tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTACGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCa                                                                                                      >  1:587239/1‑150 (MQ=255)
                                              ttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCaa                                                                                                     <  1:475013/150‑1 (MQ=255)
                                                  accaccAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAacaca                                                                                                 <  3:252437/150‑1 (MQ=255)
                                                   ccaccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtggt                                                                                                                                              >  1:1397212/1‑104 (MQ=255)
                                                   ccaccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtggt                                                                                                                                              <  2:1397212/104‑1 (MQ=255)
                                                          aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt                                                                                         >  4:948493/1‑150 (MQ=255)
                                                              aaGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCt                                                                                     >  2:437734/1‑150 (MQ=255)
                                                                   aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTg                                                                                <  4:493353/150‑1 (MQ=255)
                                                                   aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTg                                                                                <  3:734141/150‑1 (MQ=255)
                                                                    ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGc                                                                               <  1:934006/150‑1 (MQ=255)
                                                                     tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtgg                                                                                                                                               <  1:1545036/85‑1 (MQ=255)
                                                                     tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTcgc                                                                                                           <  1:1500959/121‑1 (MQ=255)
                                                                      ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                             <  1:1428900/150‑1 (MQ=255)
                                                                      ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGCTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                             <  3:309111/150‑1 (MQ=255)
                                                                       tcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACg                                                                                               >  1:163291/1‑131 (MQ=255)
                                                                       tcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt                                                                            >  2:1019419/1‑150 (MQ=255)
                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTcgc                                                                                                           >  2:1639218/1‑117 (MQ=255)
                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCg                                                                              <  1:1059159/146‑1 (MQ=255)
                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCg                                                                              >  2:1059159/1‑146 (MQ=255)
                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                                                          >  2:1516858/1‑150 (MQ=255)
                                                                          atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                         <  1:1305512/150‑1 (MQ=255)
                                                                           tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACa                                                                        <  4:1034556/150‑1 (MQ=255)
                                                                             atcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGt                                                                      >  2:559067/1‑150 (MQ=255)
                                                                              tcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGtt                                                                     <  4:541386/150‑1 (MQ=255)
                                                                              tcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGtt                                                                     >  2:1545524/1‑150 (MQ=255)
                                                                                cggcggCGAGAGGCGCCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGACTCGCGTCCAACACACCCTGTGTGGCTCCCTGCGAGAACAGTTcc                                                                   >  4:518025/1‑150 (MQ=255)
                                                                                cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTcaca                                                                                                                                                   >  1:497923/1‑70 (MQ=255)
                                                                                cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTcc                                                                   >  4:839913/1‑150 (MQ=255)
                                                                                cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTcc                                                                   >  2:974636/1‑150 (MQ=255)
                                                                                  gcggcgAGACGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCctc                                                                 >  1:1143313/1‑150 (MQ=255)
                                                                                    ggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCtt                                                               <  4:762345/150‑1 (MQ=255)
                                                                                     gcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTGCGGCCGTTCACAGTGGTGAATAACCACTCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                              <  1:70741/150‑1 (MQ=255)
                                                                                      cgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCg                                                             <  2:1489530/150‑1 (MQ=255)
                                                                                       gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCgg                                                            >  1:1684452/1‑150 (MQ=255)
                                                                                        agagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGc                                                           >  3:703682/1‑150 (MQ=255)
                                                                                         gagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGcc                                                          >  4:685340/1‑150 (MQ=255)
                                                                                           ggCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCaa                                                        >  3:899258/1‑150 (MQ=255)
                                                                                            gCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                             >  2:496441/1‑128 (MQ=255)
                                                                                              gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATc                                                     <  3:954034/150‑1 (MQ=255)
                                                                                               aCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCg                                                    >  4:305756/1‑150 (MQ=255)
                                                                                                  tctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTa                                                 >  4:896543/1‑150 (MQ=255)
                                                                                                  tctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTa                                                 <  4:625821/150‑1 (MQ=255)
                                                                                                   ctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAt                                                >  2:1167248/1‑150 (MQ=255)
                                                                                                   ctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAt                                                >  3:181053/1‑150 (MQ=255)
                                                                                                    tcCACGCAGTGCGCGACGCGCAAAATAAAGTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt                                                                                         <  2:1142340/108‑1 (MQ=255)
                                                                                                    tcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATc                                               >  3:518220/1‑150 (MQ=255)
                                                                                                     ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCa                                              <  1:393386/150‑1 (MQ=255)
                                                                                                      cACGCAGTGCGCGACGCGCAAAATAANCTTCTTCTTCGGCCGTTCANAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCAc                                             >  3:665340/1‑150 (MQ=255)
                                                                                                         gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTgg                                          <  1:325864/150‑1 (MQ=255)
                                                                                                          cAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTggg                                         <  1:559067/150‑1 (MQ=255)
                                                                                                             tGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTc                                      <  4:720840/150‑1 (MQ=255)
                                                                                                                gcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTg                                   >  2:1285162/1‑150 (MQ=255)
                                                                                                                  gtcgcgCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGtt                                 >  1:758227/3‑150 (MQ=255)
                                                                                                                  gACGCGCAAAATAAACTTCTTATTCGGCCGTTGACAGTGGTGAATAAACAATCTCACTAAGGTATAGCTGAGTCGTGTCCAACTCACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGtt                                 <  4:273173/150‑1 (MQ=255)
                                                                                                                      cgcAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt                                                                                         >  3:976069/1‑90 (MQ=255)
                                                                                                                       gcAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGt                                                   >  2:270153/1‑127 (MQ=255)
                                                                                                                       gcAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctact                            >  2:1298278/1‑150 (MQ=255)
                                                                                                                          aaaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTc                                <  4:566350/143‑1 (MQ=255)
                                                                                                                          aaaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTc                                >  3:566350/1‑143 (MQ=255)
                                                                                                                           actAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAgg                        >  1:427911/3‑150 (MQ=255)
                                                                                                                           aaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAgg                        >  1:286588/1‑150 (MQ=255)
                                                                                                                           aaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAgg                        >  2:745618/1‑150 (MQ=255)
                                                                                                                             tAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGcc                      >  2:589268/1‑150 (MQ=255)
                                                                                                                                aCTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCttt                   >  4:389973/1‑150 (MQ=255)
                                                                                                                                 cttcttcttCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCtttt                  <  3:685340/150‑1 (MQ=255)
                                                                                                                                  ttcttcttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTc                 >  1:1186019/1‑150 (MQ=255)
                                                                                                                                     ttcttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTc              >  1:1291207/1‑150 (MQ=255)
                                                                                                                                       cttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaa            <  1:113999/150‑1 (MQ=255)
                                                                                                                                          cGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaaaa         >  1:1580754/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAAC  >  NC_000913_3_pae_pgi/2867194‑2867488

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: