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% 441.8 / NA 125L71* (TTA→TGA) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (57/67);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

TAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTC  >  NC_000913_3_pae_pgi/2867194‑2867485
                                                                                                                                                  |                                                                                                                                                 
tAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTcaca                                                                                                                                                >  2:788882/1‑150 (MQ=255)
   ggTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtg                                                                                                                                             <  2:866701/150‑1 (MQ=255)
    gTCCAGCCACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtgg                                                                                                                                            <  4:561857/150‑1 (MQ=255)
    gTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtgg                                                                                                                                            >  2:501223/1‑150 (MQ=255)
    gTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtgg                                                                                                                                            <  2:1208878/150‑1 (MQ=255)
     tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtggt                                                                                                                                           <  1:668712/150‑1 (MQ=255)
     gccAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtggt                                                                                                                                           <  1:635249/149‑1 (MQ=255)
         gCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAAt                                                                                                                                       >  1:64141/1‑150 (MQ=255)
         gCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAAt                                                                                                                                       >  2:937218/1‑150 (MQ=255)
             cGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAAcc                                                                                                                                   <  1:678430/150‑1 (MQ=255)
             cGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAAcc                                                                                                                                   <  2:872746/150‑1 (MQ=255)
              gCCAGACCACGATTGCCATAACGGCGGGCATTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGCGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCa                                                                                                                                  >  2:632100/1‑150 (MQ=255)
                     cACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCAc                                                                                                                           >  2:275846/1‑150 (MQ=255)
                         aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAg                                                                                                                       <  3:436730/150‑1 (MQ=255)
                         aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAg                                                                                                                       <  1:615429/150‑1 (MQ=255)
                          ttGCCATAACGGCGGGCAATTTTTACCCCCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAgg                                                                                                                      <  3:255120/150‑1 (MQ=255)
                              cATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTaaa                                                                                                                  <  3:166092/150‑1 (MQ=255)
                                tAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGc                                                                                                                >  2:126851/1‑150 (MQ=255)
                                 aaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCt                                                                                                               >  4:14631/1‑150 (MQ=255)
                                   cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGa                                                                                                             <  2:668326/150‑1 (MQ=255)
                                   cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGa                                                                                                             <  2:362216/150‑1 (MQ=255)
                                   cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGa                                                                                                             <  1:670288/150‑1 (MQ=255)
                                     gcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGt                                                                                                           >  4:314242/1‑150 (MQ=255)
                                      cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTc                                                                                                          <  1:431114/150‑1 (MQ=255)
                                        ggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACATTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTcgc                                                                                                        >  4:384381/1‑150 (MQ=255)
                                          cAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGt                                                                                                      >  4:541054/1‑150 (MQ=255)
                                          cAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGGAAAGCTGAGTCGCGt                                                                                                      >  4:541064/1‑150 (MQ=255)
                                             tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtgg                                                                                                                                            >  1:1115392/1‑109 (MQ=255)
                                             tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtgg                                                                                                                                            <  2:1115392/109‑1 (MQ=255)
                                                ttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGGCATCTCCACGCAGTGCGCGACGCTCAAAATAAACTTCTTCTGCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAaca                                                                                                <  2:830704/150‑1 (MQ=255)
                                                ttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAaca                                                                                                <  2:332847/150‑1 (MQ=255)
                                                   ccaccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAacacac                                                                                             <  4:338119/150‑1 (MQ=255)
                                                     accaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCACAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGc                                                                                           >  1:153441/1‑150 (MQ=255)
                                                      ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                                                                                          >  1:1217041/1‑150 (MQ=255)
                                                      ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                                                                                          >  3:307965/1‑150 (MQ=255)
                                                      ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                                                                                          >  3:529249/1‑150 (MQ=255)
                                                       caGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtg                                                                                         >  2:1048304/1‑150 (MQ=255)
                                                        aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgt                                                                                        <  1:160493/150‑1 (MQ=255)
                                                          aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAACCTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt                                                                                      <  1:1222965/150‑1 (MQ=255)
                                                          aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt                                                                                      <  2:1148730/150‑1 (MQ=255)
                                                          aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtgt                                                                                      <  1:501223/150‑1 (MQ=255)
                                                            gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTgg                                                                                    >  4:362681/1‑150 (MQ=255)
                                                            gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTgg                                                                                    >  4:653415/1‑150 (MQ=255)
                                                            gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGGGCGCGACTCTCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTgg                                                                                    <  1:111533/150‑1 (MQ=255)
                                                                gTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTcc                                                                                <  1:1048304/150‑1 (MQ=255)
                                                                    ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGc                                                                            <  2:313492/150‑1 (MQ=255)
                                                                     tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCg                                                                           >  3:424607/1‑150 (MQ=255)
                                                                      ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                          <  4:283671/150‑1 (MQ=255)
                                                                      ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                          <  3:35382/150‑1 (MQ=255)
                                                                      ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCCCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                          <  1:1157857/150‑1 (MQ=255)
                                                                       gcGATCATCCGGCGCCGAGAGGCAACATCTCCACGCAGTGCGCGGCGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCACGGTAAAGCTCAGTCGCGTCCAACACCCGCCGTGTGGCTCCCTGCgct                                                                         >  2:102427/2‑148 (MQ=255)
                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                                                       >  3:636591/1‑150 (MQ=255)
                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                                                       <  3:404044/150‑1 (MQ=255)
                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                                                       >  3:321747/1‑150 (MQ=255)
                                                                           tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACa                                                                     <  1:788882/150‑1 (MQ=255)
                                                                             atcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGt                                                                   >  1:400246/1‑150 (MQ=255)
                                                                                cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTcc                                                                <  3:27443/150‑1 (MQ=255)
                                                                                    ggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCtt                                                            <  3:366070/150‑1 (MQ=255)
                                                                                     gcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                           <  1:936599/150‑1 (MQ=255)
                                                                                       gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCTGCCGTTCACAgtggtg                                                                                                                                          <  4:741692/69‑1 (MQ=255)
                                                                                       gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAgtggtg                                                                                                                                          >  3:741692/1‑69 (MQ=255)
                                                                                       gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAAt                                                                                                                                <  4:748091/79‑1 (MQ=255)
                                                                                       gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAAt                                                                                                                                >  3:748091/1‑79 (MQ=255)
                                                                                          agGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCa                                                      >  3:726877/1‑150 (MQ=255)
                                                                                            gCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTcc                                                                                <  2:4997/122‑1 (MQ=255)
                                                                                            gCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTcc                                                                                >  1:4997/1‑122 (MQ=255)
                                                                                              gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATc                                                  >  4:75236/1‑150 (MQ=255)
                                                                                                cATCTCCATGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATCTCTAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGt                                                <  3:384381/150‑1 (MQ=255)
                                                                                                cATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGt                                                <  1:92044/150‑1 (MQ=255)
                                                                                                 aTCTCCACGCAGTGCGCGACGCGCAAAATAATCTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGtt                                               <  4:494537/150‑1 (MQ=255)
                                                                                                 aTCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCAGTTCACAGAGGTGAATGACAAATCTCACCTAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGtt                                               <  3:541054/150‑1 (MQ=255)
                                                                                                 aTCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCAGCACGCAGTGCAGCATCCCCGCCCTAACGTGGGTAATGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGtt                                               <  3:541064/150‑1 (MQ=255)
                                                                                                  tctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTa                                              <  4:307965/150‑1 (MQ=255)
                                                                                                   ctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAt                                             <  1:632100/150‑1 (MQ=255)
                                                                                                       ttattttGCGCGTCGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACt                                         >  2:800823/7‑150 (MQ=255)
                                                                                                       aCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACt                                         <  2:64141/150‑1 (MQ=255)
                                                                                                       aCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCAAAGTGGGGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACt                                         >  1:32493/1‑150 (MQ=255)
                                                                                                        tattttGCGCGTCGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTg                                        >  2:1197768/6‑150 (MQ=255)
                                                                                                        cGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTg                                        <  1:126851/150‑1 (MQ=255)
                                                                                                         gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTgg                                       >  4:202410/1‑150 (MQ=255)
                                                                                                            gTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGtt                                    <  4:636591/150‑1 (MQ=255)
                                                                                                             tGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTc                                   >  4:640974/1‑150 (MQ=255)
                                                                                                               cgcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCt                                 <  1:39721/150‑1 (MQ=255)
                                                                                                                       gcAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACATACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctact                         <  1:137724/150‑1 (MQ=255)
                                                                                                                       gcAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctact                         >  3:630348/1‑150 (MQ=255)
                                                                                                                         aaaaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTaa                       >  2:448114/1‑150 (MQ=255)
                                                                                                                         aaaaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTaa                       <  2:1188726/150‑1 (MQ=255)
                                                                                                                         aaaaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTaa                       <  4:630348/150‑1 (MQ=255)
                                                                                                                          aaaTAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAg                      <  3:653415/150‑1 (MQ=255)
                                                                                                                              aaaCTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCt                  <  3:640974/150‑1 (MQ=255)
                                                                                                                                 cttcttcttCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCtttt               <  2:620646/150‑1 (MQ=255)
                                                                                                                                       cttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTc                                                                                                     <  3:187118/58‑1 (MQ=255)
                                                                                                                                       cttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTc                                                                                                     >  4:187118/1‑58 (MQ=255)
                                                                                                                                       cttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaa         >  2:940952/1‑150 (MQ=255)
                                                                                                                                        ttcGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaa        >  4:199844/1‑150 (MQ=255)
                                                                                                                                          cGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaaaa      >  1:181823/1‑150 (MQ=255)
                                                                                                                                          cGGCCGTTCACAGTGGTGAAGAACCCATCGCACCAAGGTAAAGCTGAGTCGCGTCCAACTCTCGCTGTGTGGCTCCCTGCGCTAACCGTTCCTCTTCGGGCAAATCGTTATCACTGGGTTCCTGTCCTACTAAGGCCTTTTCGTCaaaga      <  4:437893/150‑3 (MQ=255)
                                                                                                                                           ggCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAAc     >  1:464551/1‑150 (MQ=255)
                                                                                                                                           ggCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAAc     >  1:143202/1‑150 (MQ=255)
                                                                                                                                              cGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTc  <  1:873840/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTC  >  NC_000913_3_pae_pgi/2867194‑2867485

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: