Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,233,690 A→G intergenic (‑457/‑1476) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,233,6900AG88.4% 185.5 / 14.8 69intergenic (‑457/‑1476)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (7/1);  new base G (28/33);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 5.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

TATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAGCCAATTATTACAAGATTACCTCAATTGGGAAAACTACATCTTACGCAGAGTGGATTTTCCGACTAGC  >  NC_000913_3_pae_pgi/4233543‑4233829
                                                                                                                                                   |                                                                                                                                           
tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATcgca                                                                                                                                           <  2:577499/150‑1 (MQ=255)
  taaataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCAtt                                                                                                                                         <  4:151951/150‑1 (MQ=255)
     ataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTAGCGGTCAGTAGTTCAGCTAATAAATGAATCACTGCGCTAAGGGATTACACTCAACATTACGCTAACGGCACTAAAACCATCGCAttttt                                                                                                                                      >  1:503704/1‑150 (MQ=255)
       aaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCt                                                                                                                                    <  4:768323/150‑1 (MQ=255)
       aaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCt                                                                                                                                    <  1:1405611/150‑1 (MQ=255)
       aaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCt                                                                                                                                    >  1:340117/1‑150 (MQ=255)
       aaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCt                                                                                                                                    <  3:572782/150‑1 (MQ=255)
       aaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCt                                                                                                                                    <  3:349830/150‑1 (MQ=255)
         aGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCtgt                                                                                                                                  <  1:24595/150‑1 (MQ=255)
            gTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGAc                                                                                                                               <  4:667274/150‑1 (MQ=255)
               aGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTgg                                                                                                                            <  2:548981/150‑1 (MQ=255)
                 cGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGgcg                                                                                                                          >  3:647172/1‑150 (MQ=255)
                  ggggCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGgcgc                                                                                                                         <  3:649451/150‑1 (MQ=255)
                   gggCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCt                                                                                                                        <  2:883029/150‑1 (MQ=255)
                      cGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACa                                                                                                                     <  3:639174/150‑1 (MQ=255)
                      cGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACa                                                                                                                     >  3:243479/1‑150 (MQ=255)
                        gTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAAt                                                                                                                   >  1:1200390/1‑150 (MQ=255)
                           gTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCtt                                                                                                                >  2:292820/1‑150 (MQ=255)
                            tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTc                                                                                                               <  1:632370/150‑1 (MQ=255)
                            tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTc                                                                                                               <  2:680746/150‑1 (MQ=255)
                              aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCa                                                                                                             >  1:264992/1‑150 (MQ=255)
                              aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCa                                                                                                             >  4:253034/1‑150 (MQ=255)
                              aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCa                                                                                                             <  2:162616/150‑1 (MQ=255)
                              aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCa                                                                                                             >  2:251901/1‑150 (MQ=255)
                                 gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGt                                                                                                         >  1:733475/1‑150 (MQ=255)
                                    ggggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTca                                                                                                       <  3:673682/150‑1 (MQ=255)
                                    ggggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTca                                                                                                       >  1:1492175/1‑150 (MQ=255)
                                     gggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTcac                                                                                                      >  4:915208/1‑150 (MQ=255)
                                        tCATGCGGATTTTTTTTCCTTTTTTGGCAGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAAt                                                                                                   >  4:971248/1‑150 (MQ=255)
                                          aTGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCCGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtc                                                                                                 >  1:1322074/1‑150 (MQ=255)
                                            gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                                                               <  3:915208/150‑1 (MQ=255)
                                            gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                                                               >  4:706097/1‑150 (MQ=255)
                                             cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCa                                                                                              <  1:895807/150‑1 (MQ=255)
                                              ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                             >  1:1606720/1‑150 (MQ=255)
                                                 tttttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAAt                                                                                          >  3:359697/1‑150 (MQ=255)
                                                 tttttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTAACAAt                                                                                          <  4:735784/150‑1 (MQ=255)
                                                   tttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                        <  3:706097/150‑1 (MQ=255)
                                                   tttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                        <  2:1328491/150‑1 (MQ=255)
                                                    ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                                                                       >  2:1356728/1‑150 (MQ=255)
                                                     tCATCCACTGCTGGCGGTCAGTAGGTCAGCTGATAAGTGCTTCACGGCGCTAAGGGTTTACACTCAACATTTCGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCACAGTCACAATTCTCAAGATCaga                                                                                      <  2:1258137/150‑1 (MQ=255)
                                                     tCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCaga                                                                                      <  4:359697/150‑1 (MQ=255)
                                                        tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                   <  2:1642029/150‑1 (MQ=255)
                                                           aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                                                                >  1:1012340/1‑150 (MQ=255)
                                                           aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                                                                <  2:1321773/150‑1 (MQ=255)
                                                           aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                                                                >  3:1026023/1‑150 (MQ=255)
                                                           aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                                                                >  1:476006/1‑150 (MQ=255)
                                                             tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTGACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                                                              >  2:1130502/1‑150 (MQ=255)
                                                             tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCtgt                                                                                                                                  >  2:1308262/1‑98 (MQ=255)
                                                             tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCtgt                                                                                                                                  <  1:1308262/98‑1 (MQ=255)
                                                              ccTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAttt                                                                             <  1:74623/150‑2 (MQ=255)
                                                               cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtttc                                                                            >  2:1057553/1‑148 (MQ=255)
                                                                tGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTTCt                                                                           <  2:1248888/150‑1 (MQ=255)
                                                                tGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCt                                                                           >  3:969396/1‑150 (MQ=255)
                                                                 ggCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACGCTCAACATTACGCTAAGGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTACCAAAGTCACAATTCTGAGAAGCAGAAGAGTATTGCTa                                                                          >  2:1433887/1‑150 (MQ=255)
                                                                  gCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                           <  1:677539/116‑1 (MQ=255)
                                                                  gCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                           >  2:677539/1‑116 (MQ=255)
                                                                  gCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTaa                                                                         >  1:801747/1‑150 (MQ=255)
                                                                   cGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTTCTAAt                                                                        <  1:340834/150‑1 (MQ=255)
                                                                    ggTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATg                                                                       >  2:6929/1‑150 (MQ=255)
                                                                           agtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTTCTAATggccgctc                                                                <  2:25294/150‑8 (MQ=255)
                                                                                     ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                           <  1:8000/97‑1 (MQ=255)
                                                                                     ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                           >  2:8140/1‑97 (MQ=255)
                                                                                     ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                           >  2:8000/1‑97 (MQ=255)
                                                                                     ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                           <  1:8140/97‑1 (MQ=255)
                                                                                          tGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGgc                                                                                                                           <  3:333774/76‑1 (MQ=255)
                                                                                          tGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCGCATTTTTCTGTGACTGgc                                                                                                                           >  4:333774/1‑76 (MQ=255)
                                                                                                       tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAGCCAATTATTACAAGATTACCTCAATTGGGaaaa                                    <  2:801747/150‑1 (MQ=255)
                                                                                                                      aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAGCCAATTATTACAAGATTACCTCAATTGGGAAAACTACATCTTACGCag                     >  4:589278/1‑150 (MQ=255)
                                                                                                                                       aCTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAGCCAATTATTACAAGATTACCTCAATTGGGAAAACTACATCTTACGCAGAGTGGATTTTCCGACTa    >  3:595236/1‑150 (MQ=255)
                                                                                                                                         tAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAGCCAATTATTACAAGATTACCTCAATTGGGAAAACTACATCTTACGCAGAGTGGATTTTCCGACTAGc  >  1:218461/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAGCCAATTATTACAAGATTACCTCAATTGGGAAAACTACATCTTACGCAGAGTGGATTTTCCGACTAGC  >  NC_000913_3_pae_pgi/4233543‑4233829

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: