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,3400AC99.2% 431.7 / ‑4.8 120L71* (TTA→TGA) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (49/70);  total (50/70)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 120 total aligned reads displayed.

ATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCT  >  NC_000913_3_pae_pgi/2867193‑2867484
                                                                                                                                                   |                                                                                                                                                
aTAAGGGCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCACGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCAc                                                                                                                                                >  2:1060579/1‑150 (MQ=255)
    ggTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtg                                                                                                                                            <  1:652254/150‑1 (MQ=255)
      tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTTTTCTTCGGCCGTTCACAgtggt                                                                                                                                          <  2:346343/150‑1 (MQ=255)
      tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtggt                                                                                                                                          <  1:230917/150‑1 (MQ=255)
        cAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGa                                                                                                                                        >  2:1141667/1‑150 (MQ=255)
        cAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGa                                                                                                                                        <  2:1096057/150‑1 (MQ=255)
        cAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGa                                                                                                                                        <  3:628635/150‑1 (MQ=255)
         aGCGACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGaa                                                                                                                                       <  4:53435/150‑1 (MQ=255)
         aGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGaa                                                                                                                                       <  1:687126/150‑1 (MQ=255)
          gCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAAt                                                                                                                                      <  1:958721/150‑1 (MQ=255)
            aaCGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATaa                                                                                                                                    <  1:794908/150‑1 (MQ=255)
            aaCGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATaa                                                                                                                                    <  3:489/150‑1 (MQ=255)
            aaCGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCATTCACAGTGGTGAATaa                                                                                                                                    <  2:548379/150‑1 (MQ=255)
                 cAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGCGAGGCGACATCTCCACGCATGGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAAt                                                                                                                               <  3:535400/150‑1 (MQ=255)
                 cAGACCACGATAGCAATAACGGCGGGAAATGTTTACCACCAAACGCAAGTTACTCTCGATCATCAGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAAt                                                                                                                               <  2:502941/150‑1 (MQ=255)
                   gACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAAtct                                                                                                                             <  4:400565/150‑1 (MQ=255)
                     ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCa                                                                                                                           <  2:817161/150‑1 (MQ=255)
                      cACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCAc                                                                                                                          <  2:1030555/150‑1 (MQ=255)
                      cACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCAc                                                                                                                          <  3:699574/150‑1 (MQ=255)
                      cACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCAc                                                                                                                          <  1:37297/150‑1 (MQ=255)
                          aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACGAATCTCACCAAg                                                                                                                      <  1:874523/150‑1 (MQ=255)
                          aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAg                                                                                                                      <  1:102125/150‑1 (MQ=255)
                               cATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTaaa                                                                                                                 >  4:383653/1‑150 (MQ=255)
                                aTAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAg                                                                                                                <  3:665254/150‑1 (MQ=255)
                                aTAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAg                                                                                                                >  1:12108/1‑150 (MQ=255)
                                 tAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGc                                                                                                               <  4:254945/150‑1 (MQ=255)
                                    cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGa                                                                                                            <  2:220234/150‑1 (MQ=255)
                                     ggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAAcc                                                                                                                                  <  2:159199/127‑1 (MQ=255)
                                     ggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAAcc                                                                                                                                  >  1:159199/1‑127 (MQ=255)
                                       cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTc                                                                                                         >  2:1073474/1‑150 (MQ=255)
                                         ggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCaa                                                                                                                                <  1:491287/125‑1 (MQ=255)
                                         ggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCaa                                                                                                                                >  2:491287/1‑125 (MQ=255)
                                         ggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTcgc                                                                                                       <  3:20760/150‑1 (MQ=255)
                                          gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTcgcg                                                                                                      >  2:387608/1‑150 (MQ=255)
                                            aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGa                                                                                                            <  3:432638/142‑1 (MQ=255)
                                            aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGa                                                                                                            >  4:432638/1‑142 (MQ=255)
                                              tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCa                                                                                                  >  1:1144656/1‑150 (MQ=255)
                                              tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCa                                                                                                  >  3:648971/1‑150 (MQ=255)
                                               ttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCaa                                                                                                 <  2:363327/150‑1 (MQ=255)
                                               ttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCaa                                                                                                 >  3:628656/1‑150 (MQ=255)
                                                     caccaGACGCAGGTGACTCGGGATGAGCCGGGGGCGAGAGGCGACAGCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTCCTTCGGCCGTTCACAGTGGTGAATAGCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACg                                                                                           <  2:1034665/150‑1 (MQ=255)
                                                       ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGc                                                                                                               <  2:993600/128‑1 (MQ=255)
                                                       ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGc                                                                                                               >  1:993600/1‑128 (MQ=255)
                                                       ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                                                                                         <  4:397732/150‑1 (MQ=255)
                                                       ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                                                                                         >  1:1070915/1‑150 (MQ=255)
                                                         aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgt                                                                                       <  2:1070837/150‑1 (MQ=255)
                                                           aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt                                                                                     <  4:16645/150‑1 (MQ=255)
                                                           aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt                                                                                     <  4:628656/150‑1 (MQ=255)
                                                           aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt                                                                                     >  1:674405/1‑150 (MQ=255)
                                                            cGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgtg                                                                                    >  4:156561/1‑150 (MQ=255)
                                                             gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTgg                                                                                   <  1:638356/150‑1 (MQ=255)
                                                             gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTgg                                                                                   <  4:60590/150‑1 (MQ=255)
                                                              cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGc                                                                                  <  2:949672/150‑1 (MQ=255)
                                                                aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTc                                                                                >  2:969076/1‑150 (MQ=255)
                                                                    aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTg                                                                            <  2:12108/150‑1 (MQ=255)
                                                                     ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGc                                                                           <  2:304392/150‑1 (MQ=255)
                                                                      tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTaaa                                                                                                                 >  3:474037/1‑111 (MQ=255)
                                                                      tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTaaa                                                                                                                 <  4:474037/111‑1 (MQ=255)
                                                                        tcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt                                                                        >  1:446892/1‑150 (MQ=255)
                                                                          gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                                                      <  2:1107073/150‑1 (MQ=255)
                                                                          gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                                                      >  1:1122377/1‑150 (MQ=255)
                                                                           atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGc                                                                           <  1:10779/144‑1 (MQ=255)
                                                                           atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                     <  1:1073474/150‑1 (MQ=255)
                                                                           atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                     >  3:651812/1‑150 (MQ=255)
                                                                               tcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGtt                                                                 >  2:571244/1‑150 (MQ=255)
                                                                                ccGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTc                                                                <  1:1130255/150‑1 (MQ=255)
                                                                                      gcgcGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                          <  2:798254/150‑1 (MQ=255)
                                                                                       cgAGAGGCGACATCTCCGCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCg                                                         >  1:136002/1‑150 (MQ=255)
                                                                                        gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCgg                                                        <  3:261674/150‑1 (MQ=255)
                                                                                        gagaGGCGACATCTCCACACAGTGCGCGACGTGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACTCTCTGTGTTGCTCCCTGCGATAACAGTTCCTCTTCgg                                                        >  2:174543/1‑150 (MQ=255)
                                                                                           agGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCa                                                     >  2:234833/1‑150 (MQ=255)
                                                                                                  aTCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGtt                                              <  1:587521/150‑1 (MQ=255)
                                                                                                   tctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTa                                             <  2:1144656/150‑1 (MQ=255)
                                                                                                    ctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCa                                                                                                                           >  1:87127/1‑71 (MQ=255)
                                                                                                    ctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCa                                                                                                                           <  2:87127/71‑1 (MQ=255)
                                                                                                    ctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAt                                            <  1:387608/150‑1 (MQ=255)
                                                                                                    ctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAt                                            >  2:915116/1‑150 (MQ=255)
                                                                                                     tcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                         >  1:166515/1‑120 (MQ=255)
                                                                                                     tcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                         <  2:166515/120‑1 (MQ=255)
                                                                                                     tcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTc                                                                <  4:491626/129‑1 (MQ=255)
                                                                                                     tcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTc                                                                >  3:491626/1‑129 (MQ=255)
                                                                                                       cACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCAc                                         >  1:545359/1‑150 (MQ=255)
                                                                                                        aCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACg                                                                                           <  1:283759/99‑1 (MQ=255)
                                                                                                        aCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACt                                        <  1:1141667/150‑1 (MQ=255)
                                                                                                           ttttGCGCGTCGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTtgg                                     >  2:709149/4‑150 (MQ=255)
                                                                                                                  cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGGAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGt                              >  2:388397/1‑150 (MQ=255)
                                                                                                                    aCGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTAGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTc                            <  3:160823/150‑1 (MQ=255)
                                                                                                                     cgcgCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTct                           <  1:101967/150‑1 (MQ=255)
                                                                                                                           aaaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGa                                                                                                            >  4:620528/1‑63 (MQ=255)
                                                                                                                            aaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAacaca                                                                                             >  1:749275/1‑77 (MQ=255)
                                                                                                                            aaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAacaca                                                                                             <  2:749275/77‑1 (MQ=255)
                                                                                                                                aaCTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCtt                >  3:493509/1‑150 (MQ=255)
                                                                                                                                   ttcttcttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCTAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTc                            <  1:773191/135‑1 (MQ=255)
                                                                                                                                   ttcttcttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCTAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTc                            >  2:773191/1‑135 (MQ=255)
                                                                                                                                   ttcttcttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTc             <  1:809151/150‑1 (MQ=255)
                                                                                                                                         ttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaa       >  3:90629/1‑150 (MQ=255)
                                                                                                                                          tcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaaa      <  3:212308/150‑1 (MQ=255)
                                                                                                                                          tcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaaa      <  2:545359/150‑1 (MQ=255)
                                                                                                                                              ccGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCt  <  1:969076/150‑1 (MQ=255)
                                                                                                                                                gTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCCTCGGCCAAATCGTTATCACTGGGTTCCTg                               >  1:67693/1‑119 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCT  >  NC_000913_3_pae_pgi/2867193‑2867484

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: ‑