Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,234,667 Δ1 bp intergenic (‑1434/‑499) 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,234,6660A.96.0% 301.2 / 11.6 79intergenic (‑1433/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base . (37/34);  minor base G (0/3);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_pae_pgi/4234523‑4234809
                                                                                                                                               |                                                                                                                                               
atataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatcatc                                                                                                                                           <  1:117515/149‑1 (MQ=255)
  ataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAAATGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgt                                                                                                                                         <  1:167806/149‑1 (MQ=255)
   taCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCCTTGAAGAGGCGCTTGCAGTGATCAtcgtc                                                                                                                                        >  1:468761/1‑149 (MQ=255)
      gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGat                                                                                                                                     <  1:888905/149‑1 (MQ=255)
      gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGag                                                                                                                                     >  1:128030/1‑148 (MQ=255)
         gTAGAGATGCACCCCAAGGTCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGGTTAGAGGCGTAAATCCCCTGGCCACCCCGAAAATAGATGCATTGAAGCGGCGCTTGCGGTGTTCATCGTCGatgaa                                                                                                                                  >  2:1023616/1‑146 (MQ=255)
          tAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGt                                                                                                                                 >  2:847917/1‑149 (MQ=255)
              gaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGc                                                                                                                             <  1:1164232/149‑1 (MQ=255)
                tGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATATGGGCATCGCATCTGTAAAAATGACAAAGATCCAAGGCGTAAATCTCCTGCCCACCCCGAGAATCGACGCACTGCAGCGGCGCATGGAGTGATCTTCGTCGATTTGTaagtat                                                                                                                           >  2:648870/1‑144 (MQ=255)
                   cacccAAACGCGTTTAATCCTTATCACGAGATGTGCGCACCGCCTCTGCAAAACTAGCAGAGATTAGAGGCGTAAACCCCCCGGTCACCCCGAGAATAGATGCATTGAATAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGAt                                                                                                                        <  2:407996/147‑1 (MQ=255)
                     cctaaGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      <  1:1053039/146‑1 (MQ=255)
                     cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAAATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGACGCTTGCAGTGATCACCGTCGCTATGTAGGCCGGATaa                                                                                                                      >  1:229193/1‑149 (MQ=255)
                     accAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTCATCATCGTCGCTATATAGGCCGGATaa                                                                                                                      >  2:64028/2‑149 (MQ=255)
                      ccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     >  1:511878/1‑149 (MQ=255)
                      ccAAGGGCGTTTTATCCTTGTTACGAGATGGGGGACTCGCCCCTGTAAAACTGGAAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGTATTTAAGCGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     <  1:225495/149‑1 (MQ=255)
                           ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                <  1:847917/149‑1 (MQ=255)
                           ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                <  2:1521/149‑1 (MQ=255)
                           ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTCTAAAAATGGCATAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTCGTCATGGTTGATATGTAGGCAGGATATGGCGtt                                                                                                                >  2:707786/1‑149 (MQ=255)
                               ttttGTCCGTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcg                                                                                                            <  2:792432/149‑1 (MQ=255)
                               ttttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcg                                                                                                            >  1:8020/1‑149 (MQ=255)
                                    tCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGTGATTAGAGGCGTAAATCCCCTGGCCACTCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGTCGGTTTAGGCGTTCACGGCGCa                                                                                                       >  1:866680/1‑149 (MQ=255)
                                    tCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                       >  1:315058/1‑149 (MQ=255)
                                      cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                     >  1:837109/1‑149 (MQ=255)
                                             cGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              <  1:1210022/149‑1 (MQ=255)
                                                gaTGTGGCGATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                           <  1:409616/149‑1 (MQ=255)
                                                 aTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTTGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCTTCGATATGTAGACCGGATTAGTCGTTCACGCCGCATCCGGCAACCGAt                                                                                          >  2:409616/1‑149 (MQ=255)
                                                   gtgGGCATCGCCTCTGTAAAACTGGCAAAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGTTGc                                                                                        >  2:1189448/1‑149 (MQ=255)
                                                     gggCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCTGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                      <  1:1189837/149‑1 (MQ=255)
                                                       gagtcctctctGTAAAACTGGCATGGGTTTGAGGCGTAAATCCCCTGGCCACCTCGAAAATAGATGCATTGAAGAGGCGTTTGCAGTCATCATCGTCGATATGTAGGCCGGATAAGACGTTCACGCCGCGTCCGGCAACCGATGCCTGa                                                                                    <  2:81325/142‑1 (MQ=255)
                                                         aTCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                  <  2:314717/149‑1 (MQ=255)
                                                          tCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGc                                                                                 >  2:143920/1‑149 (MQ=255)
                                                             ccTCTGTAAACCTGGCAGAGATTAGAGGCGAAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                              <  1:1014929/149‑1 (MQ=255)
                                                             ccTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTGGCAGTGATCATGGCCGCTATGTTGGCCGGAGCAGGCGTTCACGCCACATACGGCAACCGATGACCGAttctga                                                                              >  2:459198/1‑144 (MQ=255)
                                                              ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgtc                                                                                                                                        <  1:315898/90‑1 (MQ=255)
                                                              ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgtc                                                                                                                                        >  2:315898/1‑90 (MQ=255)
                                                               acagaaaaaCCGGCAAAGATAAGACGCGTAAATCCCCTGACCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGAAGTGATCATCGACCATATGAAGGCCAGCAAATTCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                            <  2:866680/144‑1 (MQ=255)
                                                                   tAAAACTGGCAGGTATTAGAGGCGTAAATCTCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                        <  1:674001/149‑1 (MQ=255)
                                                                        cTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                                   <  2:315058/149‑1 (MQ=255)
                                                                           gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCTTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  2:381201/1‑149 (MQ=255)
                                                                           gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  1:468654/1‑149 (MQ=255)
                                                                             agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgcatc                                                              <  2:814669/149‑6 (MQ=38)
                                                                             agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCAGCAACCGATGCCTGATGCAACGCGGTCGCGCCTTATc                                                              >  2:684735/1‑149 (MQ=38)
                                                                             agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                              >  1:1036142/1‑149 (MQ=255)
                                                                             agagATTAGAAGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAATGGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                              <  2:740595/149‑1 (MQ=255)
                                                                              gagaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             >  1:1001013/1‑149 (MQ=255)
                                                                                gaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                           <  1:143920/149‑1 (MQ=255)
                                                                                 aTTAGAGGCGTTAATCCCTTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                          <  1:918439/149‑1 (MQ=255)
                                                                                 aTTAGAGGCGTAAACCCCCTGGGCACCCCGAGAATCGGTGCATTGAAGAGGCGCTTGCAGTTATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                          <  2:1036142/149‑1 (MQ=37)
                                                                                          gTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                 >  1:788164/1‑149 (MQ=255)
                                                                                            aaaTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               <  2:837109/149‑1 (MQ=255)
                                                                                            aaaGCCCCTGGCCACTCCGATAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGACAACCGAGGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               <  1:474580/149‑1 (MQ=37)
                                                                                                  ccTGGCGACCCCGGGAATAGATGCATTTCAGAGTCTCCTGCAGGGGTCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                         <  1:707786/149‑1 (MQ=21)
                                                                                                    tGGCCACCGCGAGAATAGATGCATTGAAGAGGCGTTTGCAGTGACCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                       <  2:534930/149‑1 (MQ=255)
                                                                                                    tGGCCACCCCGAGAATAGATGCATTGGAGAGGCGCTGGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                       <  2:591351/149‑1 (MQ=255)
                                                                                                       ccACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGAAAAGGAGTTCACGCCGCATCCGGCAACCGATGCCTTATGCGACGCGGCCGCGTCTTACCAGGCCTACAGGTCGATGCCCattttt                                    >  2:428256/1‑145 (MQ=255)
                                                                                                          ccccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                                 >  2:1138773/1‑149 (MQ=255)
                                                                                                           cccGAGAATAGATGCATTGACGAGGCGCTTGCAGGGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                <  2:637248/149‑1 (MQ=255)
                                                                                                           cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                >  2:1124581/1‑149 (MQ=255)
                                                                                                            ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                               <  2:932975/85‑1 (MQ=255)
                                                                                                            ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                               >  1:932975/1‑85 (MQ=255)
                                                                                                             cGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                              >  1:334404/1‑149 (MQ=255)
                                                                                                               agaATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTa                            <  1:1138773/149‑1 (MQ=255)
                                                                                                                    aGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTCTCAGGCCTACAGGTCGATGCCCATATGTCAATCGTATGCgg                       >  1:498256/1‑149 (MQ=25)
                                                                                                                      aTGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                                                               <  2:356470/91‑1 (MQ=255)
                                                                                                                      aTGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                                                               >  1:356470/1‑91 (MQ=255)
                                                                                                                      aTGCATTGAAGAGGCGCTTGCAGTGATCATCGGCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCTTGATGCGACGCGGTCGCGTCTTATCACGCTTAAAGGTTGATGACATTTTGTAAATCGGATTTggga                     >  2:683551/1‑147 (MQ=21)
                                                                                                                       cgCATTGAAGGGGTTCTTGCAGGGATCATCGGCGATATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                    <  1:356309/148‑1 (MQ=25)
                                                                                                                        gCATTGAAGGGGCGCTTGCGGTGATCATAGTCGATGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCGGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAt                   <  1:64821/149‑1 (MQ=25)
                                                                                                                         cATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTAGGTAAtt                  >  1:152131/1‑149 (MQ=255)
                                                                                                                         cATTGAAGAGGCACTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                  <  1:811693/149‑1 (MQ=255)
                                                                                                                          aTTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                 >  2:619311/1‑149 (MQ=255)
                                                                                                                          aTTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                 >  1:575542/1‑149 (MQ=255)
                                                                                                                          aTTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTAACGCCGCCTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                 >  1:480897/1‑149 (MQ=255)
                                                                                                                              aaGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTCCATCGTATTCGGCAATTAATAc             >  1:355758/1‑149 (MQ=255)
                                                                                                                                 gggCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa          <  2:1001013/148‑1 (MQ=255)
                                                                                                                                 agggccTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa          <  2:137213/144‑1 (MQ=255)
                                                                                                                                  ggCGCTTGCAGTGAACATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg         <  1:1065360/149‑1 (MQ=255)
                                                                                                                                   gcgcTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGc        >  1:275090/1‑149 (MQ=255)
                                                                                                                                        tGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGGGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat   <  2:509539/149‑1 (MQ=255)
                                                                                                                                         gCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgatt  >  1:659250/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_pae_pgi/4234523‑4234809

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: