Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_pgi 3,815,809 Δ1 bp 19.4% intergenic (‑41/+25) pyrE ← / ← rph orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,815,8090G.19.4% 147.2 / 46.4 67intergenic (‑41/+25)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base . (6/7);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

GAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCC  >  NC_000913_3_pae_pgi/3815664‑3815946
                                                                                                                                                 |                                                                                                                                          
gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCg                                                                                                                                         >  2:21873/1‑149 (MQ=255)
 aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTACTCTTACTTTTCTCCAGCCAAAAAAAAGGCGa                                                                                                                                        >  2:487562/1‑149 (MQ=255)
 aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTAAGGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAACAGGGGa                                                                                                                                        >  2:914135/1‑149 (MQ=255)
   gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTCAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCTTCTTACTTTTCTCCATCACAAAAAAAGTCGACt                                                                                                                                      >  1:167991/1‑149 (MQ=255)
     aaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtca                                                                                                                                    <  1:609190/149‑1 (MQ=255)
        aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatca                                                                                                                                 >  1:836577/1‑149 (MQ=255)
              atttGCGCCCGTGTTTCAGCGAAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTGTAACGTCACAAAAAAGGCGAATCATCAGTcccc                                                                                                                           >  2:373612/2‑149 (MQ=255)
                     cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaaa                                                                                                                    >  2:83698/1‑149 (MQ=255)
                      ccGGATTTCGGCGGATGCTCGCGAGACTTTAAAACCTGCTTGCTAAGCGCAAATGCAATAAACTGGCGCTGATGTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAACAAGGCGCCTAATCCGTCCCACTAAAAAt                                                                                                                   >  2:241607/1‑149 (MQ=255)
                      ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGACTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                                                   >  2:474026/1‑149 (MQ=255)
                       cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCAACTCATCAGTCGCCTTAAAAATc                                                                                                                  >  1:283264/1‑149 (MQ=255)
                            ttCTGCGTAAACTCTCCAAATTTTAACATCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                                                             <  1:980504/149‑1 (MQ=255)
                               aGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                         <  1:21873/149‑1 (MQ=255)
                               aGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                         >  2:788396/1‑149 (MQ=255)
                                            cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCT                                                                                            >  1:600336‑M1/1‑149 (MQ=255)
                                             aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCgc                                                                                           >  2:224867‑M1/1‑147 (MQ=255)
                                             aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                                                           >  2:978352‑M1/1‑149 (MQ=255)
                                             aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                          >  1:971026‑M1/1‑149 (MQ=255)
                                                    aaaacCTGCTTGCTAAGCCCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTa                                                                                   <  2:288565/146‑1 (MQ=255)
                                                    aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCGGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACTAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTa                                                                                   <  1:777245/149‑1 (MQ=255)
                                                       acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCCCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                                 >  1:362882‑M1/1‑140 (MQ=255)
                                                        ccTGTTTGCTAAGCGCAAAGTCTATAGACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTATAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCTAGCGTCGCCTTCTGCGTCGCTACaa                                                                                <  2:820160/149‑1 (MQ=255)
                                                        acTGCTCGCAAAGCGCAAATTCAACAAACTGGCGCTGATACTGTTTCATTCCTTCGCCCCTCATCTTACTTTTCTACATACCAAAAAAAGGCGACTCATCAGTCGCATTAAAAATCAGTTCCCCAGCGCCGCCTTCTGCGTCGCTACaa                                                                                <  2:491646/148‑1 (MQ=255)
                                                         cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt                                                                               >  2:580595/1‑149 (MQ=255)
                                                           gcttgctAAGCGCAAATTCAATACACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATgg                                                                             <  1:644511/149‑1 (MQ=255)
                                                             tgtgtaagcCCAACTTCAGTAAACTAGGGCTTATATTGTGTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATTAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                                                           <  2:1001009‑M1/145‑16 (MQ=255)
                                                              tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                          <  1:580595/149‑1 (MQ=255)
                                                               gctAAGCGCAAATTAAATAAACTGGGGCTGATATGGTTTCATGGCTTCGCGCCTCATCTTACTTTTCTAAAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                         <  2:18125/149‑1 (MQ=255)
                                                                ctAAGCTCAAATTCAATTGACTGGCGCTGATATGTTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                        <  1:621924‑M1/149‑19 (MQ=255)
                                                                 tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGATAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                       <  1:755976/149‑1 (MQ=255)
                                                                 tAAGCGCAAATTCAATAAACTGGCGCTGACATGGTTTCATGCCTTCGCACCTCATCGTACTTTTCTACAGACAAAAAAAAGCGACCCATCAGTCGCCTTAAAAATCAGCTTGCCAGCGCCGCCCTCTGCGTCGCTACCATGgatccgat                                                                      >  1:202154/1‑149 (MQ=255)
                                                                     cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAAt                                                                                                                   >  2:513104/1‑101 (MQ=255)
                                                                     cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAAt                                                                                                                   <  1:513104/101‑1 (MQ=255)
                                                                       cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                        >  1:84205‑M1/1‑124 (MQ=255)
                                                                       cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                        <  2:84205‑M1/142‑19 (MQ=255)
                                                                        aaaTTCAATAAACTGGCGCTGATATGGTTCCAAACCCTGGGTCCTCATCTTACTCTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGCTTCGATTCCCCTc                                                               >  2:468143/1‑149 (MQ=255)
                                                                           ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTCTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg                                                             <  1:984358/149‑1 (MQ=255)
                                                                            tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtc                                                            >  1:148371‑M1/1‑119 (MQ=255)
                                                                             caactaccTGTCGCTTACATGGTTGCATGCCTTGGCTCCTCATCTTTCTTTTCTACAGACAAAAAAAAGGCGACTCATCCGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTTCGTCGCTACATTGGATTCGATTCCCCTCCGGc                                                           <  2:407317/142‑1 (MQ=255)
                                                                                aaacatggcgcTGATCTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                        <  1:226491‑M1/144‑35 (MQ=255)
                                                                                  aaCTGGCGCTGTTATGGTTTCATGCCTTCGCTCCTCATCTTATTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCagag                                                      <  1:347032/149‑1 (MQ=255)
                                                                                      ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaa                                                  <  2:635208/149‑1 (MQ=255)
                                                                                         tcTGATATCGTGTCATGGCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATTCCCCTCGGGCCAGAGCcaacaa                                               <  2:316072/148‑1 (MQ=255)
                                                                                         gcTGATATGGTTACATAACTTAGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctggattctaattttccgcccggcttcggtcatcact                                               >  2:460204‑M1/1‑106 (MQ=255)
                                                                                          gtgtgatGGTTTTACTCCTCCGCTGCTCATCTTACATTTCTATAGACAAAAGAAAGCGACTCATCAGTCGCCTCAAAAATCCGTTTGCCAGCGCCGCCTTCTGCGTGGCTACAATGTATTCGATTCCCCTCGGGCCAGAGCAAACAAGa                                             <  2:952585/144‑1 (MQ=255)
                                                                                           tGATATGGTTTTATGCCTTCGCTCCTCATCTTACTTTTCTACAGAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGa                                             <  2:354708/149‑1 (MQ=255)
                                                                                            gATATTGTGTTATGACTTCGCTCCTCAACTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAt                                            <  2:271547/149‑1 (MQ=255)
                                                                                             ttatTGTTTCCAGCCTTCGCTGCTTATCTTAAGTTTCTACAGACATAAAAAATGCGATTTATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCTATTTCCCTCGGGCCAGAGCCAACAAGATg                                           <  2:630379/148‑1 (MQ=255)
                                                                                                        tGCCTTCGCTCTTCATCTTTCTTTTCTACTGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGACAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                <  1:260017/149‑1 (MQ=255)
                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGt                                        <  2:297227/137‑1 (MQ=255)
                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGt                                        >  1:297227/1‑137 (MQ=255)
                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGCGGCAGAGCCAACAAGATGAGTAGCTCTTCATgg                            <  2:971026/149‑1 (MQ=255)
                                                                                                                ctcctcATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGaa                       <  2:202154/149‑1 (MQ=255)
                                                                                                                  cgtcATCTTACTTTTCTACAGAAAAAAAAAAGGCGACTCATCAGCTGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTTTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAc                      <  1:136763/147‑1 (MQ=255)
                                                                                                                   ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                     <  2:148371‑M1/149‑70 (MQ=255)
                                                                                                                   ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATTGATTCGATTCCCCTCGGGCCAGAGCCACCAAGCTGATTAGCTCTTCATGGGTGCAcc                     >  2:236526/1‑148 (MQ=255)
                                                                                                                   ctcATCTTACTTTTCTACAGAAAAAATAAAGGCGACTCATCAGTCGCCTTAAAATCAGGCCGTGAAGGGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                     <  1:484157‑M1/149‑70 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                        >  2:539692‑M1/1‑70 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                        <  1:539692‑M1/104‑35 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCACCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCAGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc           >  1:282988/1‑149 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCAAGCGCCGCCTTATGAcgaccactaaacttcaaggaatcccccgtcttcggtcatcactacgttcatgtcggtctctgcggcaaagtcttcaacg           >  2:362806‑M1/1‑70 (MQ=255)
                                                                                                                             ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCGCAAAAATCCGTTTGCCAGCGCCACCTACTGCGTAGCTACAATGGATTCGATTACCCTCGGGCCAGAGCCAACAAGATGAGTACCTCTTCACCGGTCACCGCCCCGCCtac           >  2:223097/1‑147 (MQ=255)
                                                                                                                               atCTACAGAGAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCCCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg         <  2:927041/148‑1 (MQ=255)
                                                                                                                                  tACAGACAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAGTGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGGTCTTCATGGGTGAACGGCTCGCCTTCTGCCGt     >  2:518848/1‑149 (MQ=255)
                                                                                                                                   aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGt     <  2:525043/149‑1 (MQ=255)
                                                                                                                                     aGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAACTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTccc  >  1:197271/1‑149 (MQ=255)
                                                                                                                                      gACAAAAAAAAGGCGACTCATCAGCCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc  <  2:1000636/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
GAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCC  >  NC_000913_3_pae_pgi/3815664‑3815946

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: