Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,110,747 G→T intergenic (+62/+784) cdh → / ← yiiQ CDP‑diacylglycerol phosphotidylhydrolase/DUF1454 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,110,7470GT95.7% 205.6 / ‑3.8 70intergenic (+62/+784)cdh/yiiQCDP‑diacylglycerol phosphotidylhydrolase/DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (38/28);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCC  >  NC_000913_3_pae_tpiA/4110601‑4110888
                                                                                                                                                  |                                                                                                                                             
gCTGGCGACACAACGAAATCTACTGACGCTTAACCGTTGTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAggctgc                                                                                                                                             <  2:422233/149‑1 (MQ=255)
  tGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGccc                                                                                                                                           <  2:70944/149‑1 (MQ=255)
          caACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGa                                                                                                                                   <  1:521219/149‑1 (MQ=255)
          caACGAAAGCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGa                                                                                                                                   <  1:1946/149‑1 (MQ=255)
           aaCGAAATCTCCTTACGCTGAACCGGGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAg                                                                                                                                  <  2:264589/149‑1 (MQ=255)
           aaCGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTTTGAGATGTGGCGTTCAGGACACGACATAAAGAAATGAGGGGCCTGATATGAAATCCGGACCCTGTCAAACTTAAGGCTGCCCCGCCGCCAg                                                                                                                                  >  2:62183/1‑149 (MQ=255)
           aaCGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAg                                                                                                                                  >  1:242862/1‑149 (MQ=255)
                  tCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGAGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGACCGCTTTc                                                                                                                           >  2:477734/1‑149 (MQ=255)
                    tACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCTAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTGGc                                                                                                                         >  1:224891/1‑149 (MQ=255)
                    tACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGc                                                                                                                         >  1:224886/1‑149 (MQ=255)
                    tACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGc                                                                                                                         >  1:421034/1‑149 (MQ=255)
                      cTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCa                                                                                                                       >  1:351541/1‑149 (MQ=255)
                       tGACGCTTAACCGTGCTTCAGCCGAGGAAATGCAGGATCATCAGTGTGAGATTTTTCGGTAAGGGGAAGTGTGGAGGAAAGTAAGTGCCGGATATGAAATCAGGCAGCTGTCAGACTTAAGGCTGGCCCGCCGCGAGCGGTTTCGCCaa                                                                                                                      >  2:466822/1‑149 (MQ=255)
                            cTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTa                                                                                                                 <  2:351541/149‑1 (MQ=255)
                             ttAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTaa                                                                                                                >  1:88002/1‑149 (MQ=255)
                              tAACGGTGCTTCAGCCGAGGAAATTCAGGATCCTCAGTGGGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAg                                                                                                               <  2:387839/149‑1 (MQ=255)
                                   gTGCTTCAGCCCAGCAAATCTAGTATCATCAGTGTCAGATTTTGCGTTAGGGATTAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCCCCTGTCAGAACTAAGGCTGCCCCGCCGCGTGCGGTTTCGCCCCTGCTAAGAtagt                                                                                                          >  2:226246/1‑146 (MQ=255)
                                       ttCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGa                                                                                                      >  2:327103/1‑149 (MQ=255)
                                          aGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGatat                                                                                                   >  2:44412/1‑149 (MQ=255)
                                          aGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGatat                                                                                                   <  2:232111/149‑1 (MQ=255)
                                           gCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATg                                                                                                  <  1:359200/149‑1 (MQ=255)
                                           cccgaacaaaTTCAGGAACATCAATGTCATATTTTGCATTAAGGCTAAGAGTTAAGGAAATTAAGTGCCGGATATGAAATCCGTCACCTGTCAGACTTAACACCGCCCAGCGTCAAGCGGTTTCGCCAATTCTAAGATCTTAAGATATg                                                                                                  <  1:278303/142‑1 (MQ=255)
                                            ccGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCATCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATgg                                                                                                 >  1:505405/1‑149 (MQ=255)
                                            ccGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAGGATATgg                                                                                                 >  2:562003/1‑149 (MQ=255)
                                                ggAAATTCAGGATCATCAGTGTGAGATGTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCt                                                                                             <  1:295459/149‑1 (MQ=255)
                                                  aaaTTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCttt                                                                                           <  2:26404/149‑1 (MQ=255)
                                                     ttCAGGATCATCAGTGTGAGCTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTACGATATGGGTCTTTTGc                                                                                        >  2:493000/1‑149 (MQ=255)
                                                               tcaGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCt                                                                              >  2:317629/1‑149 (MQ=255)
                                                                caGTGTGAGATTTTGCGTTAAGGCGAAGAGGTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCCGTCAGCCTTAAGGCTGCCCCACCGCGAGCGGTTTCGCCAATTCTAAGACCTTAAGATATGGGGCTTTTGCCTTCCCGGCTg                                                                             >  2:159947/1‑149 (MQ=255)
                                                                  gtgtGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTgcg                                                                           <  1:285150/149‑1 (MQ=255)
                                                                   tgtgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTgcgc                                                                          <  2:421034/149‑1 (MQ=255)
                                                                     tgAGATTTTGCGTTAAGGCGGAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTgcgcgt                                                                        >  1:151435/1‑147 (MQ=255)
                                                                     tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGc                                                                                        <  1:293489/133‑1 (MQ=255)
                                                                     tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGc                                                                                        >  2:293489/1‑133 (MQ=255)
                                                                     tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCtt                                                                        >  1:465298/1‑149 (MQ=255)
                                                                     tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGAAGAAGTAAGATAGAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCtt                                                                        >  1:250918/1‑149 (MQ=255)
                                                                      gagaTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTg                                                                       >  1:177472/1‑149 (MQ=255)
                                                                            ttGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCa                                                                 >  2:22132/1‑149 (MQ=255)
                                                                               cGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAaca                                                              <  1:288496/149‑1 (MQ=255)
                                                                                 ttAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACt                                                            >  1:402709/1‑149 (MQ=255)
                                                                                 gaagggCGAAGAGTTAAGGAAAGTAAGTGCGTGATATGAAATCAGGTACCTGTCAGACTTAAGGCTGTCCTGCCGCGATCGGTTTCGCCAATTCTAAGATCTTAATATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACt                                                            <  2:339695/145‑1 (MQ=255)
                                                                                  tAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTc                                                           <  1:327103/149‑1 (MQ=255)
                                                                                    gggCGAAGAGTTAAGGAAAGTAAGGGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGGGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTccc                                                         <  1:42960/148‑1 (MQ=255)
                                                                                      gCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCt                                                       >  1:355629/1‑149 (MQ=255)
                                                                                           gagTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGc                                                  <  1:317629/149‑1 (MQ=255)
                                                                                              ttAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTCAGGCTGCCCCGCCGCGCGCGGCTTCGCCAATTCTAAGATCTACAGATATGGGGCTTTTGGCTCTACGTCTTCTCTTTCCCGCAGCCCTCCCCCCGTGCCGa                                               >  2:192436/1‑149 (MQ=255)
                                                                                                 aGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTCAGGCTGCACAGCCGAGAGCGGTTTCGCCAATTCTAGGATCTTCAGATAGGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGtt                                            <  1:22132/149‑1 (MQ=255)
                                                                                                    aaaGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTaa                                                                                                        >  2:441431/1‑86 (MQ=255)
                                                                                                    aaaGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTaa                                                                                                        <  1:441431/86‑1 (MQ=255)
                                                                                                    aaaGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTCAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTaa                                         >  1:22337/1‑149 (MQ=255)
                                                                                                        taagtGCCGGATATGAAATCCGGCACCTGTCAGACTTGAGGGGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAATATCTTAAGATATGGGTCTTCTCCCTACCAGGATGCGATTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAgcgc                                     >  2:428937/1‑149 (MQ=255)
                                                                                                        taagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAgcgc                                     <  2:88002/149‑1 (MQ=255)
                                                                                                            tGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAgcgc                                     >  2:775/1‑145 (MQ=255)
                                                                                                            tGCCGGATTTCATATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAgcgc                                     <  1:775/145‑1 (MQ=255)
                                                                                                               cGGATATGAAATCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTGAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCa                              >  1:345422/1‑149 (MQ=255)
                                                                                                                 gtaatGGAATCCGGCACCTGTCAGACTTAAGGCTGCCCTGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACCCTCTCCTCGTGCCGAGTTTAAGCGCGGCTTCAAc                            <  2:325064/146‑1 (MQ=255)
                                                                                                                          tCCGGCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTcc                   <  2:215423/149‑1 (MQ=255)
                                                                                                                              gCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTc                                                      <  1:25578/110‑1 (MQ=255)
                                                                                                                              gCACCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTc                                                      >  2:25578/1‑110 (MQ=255)
                                                                                                                                aCCTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCGTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCCCCCCCCtctc             >  2:144992/1‑149 (MQ=255)
                                                                                                                                aCCTGTCAGACTTAAGGCTGCCCCACCGCGAGCGGTTTCGCCAATTCCAAGAACTTAAGATACGGGTCTTTTGATTTAACGGCTGCGCTTGCCAGCAACTCTCCCCTCGTGCTGAGTTTAAGAGCGGCTTCAACGTCCCCTcccccaac             >  2:477146/1‑143 (MQ=255)
                                                                                                                                  ccgTCAGAATTAAGTCTGCCCTGCCGCGAGGCGTTTCGCCAAGTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGTGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCaa           <  2:504759/147‑1 (MQ=255)
                                                                                                                                  cTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCa                                                                                   <  1:404381/77‑1 (MQ=255)
                                                                                                                                  cTGTCAGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCa                                                                                   >  2:404381/1‑77 (MQ=255)
                                                                                                                                    gTCAGACTTTAGGCGGCCCCGCCGCGGGCGGTTTCGCCCATTCTAAGTTCTTAAGGTATGGGTCTTTTGCCTTCACGGGTGCGCTTGCCAGCAACACTCCCTTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATa         <  2:152815/149‑1 (MQ=255)
                                                                                                                                     tcgagcTTAAGGCTGCCCCGCCACGACCGGGTTCGCCAATTCTAAGATCTTACGATATGGGTCTTTTGCCTTCACGGCGGCGCTTGAAAGCAACCCTCCCCCCGTGCTGAGTTTAAGCGCGGCTTCATCGTCCTCTCCGCTCTCAATAc        <  2:32990/144‑1 (MQ=255)
                                                                                                                                       aGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCGGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATAccc      >  1:367915/1‑149 (MQ=255)
                                                                                                                                       aGACTTAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATAccc      >  1:379175/1‑149 (MQ=255)
                                                                                                                                         aCTTAAGGCTGCCCCGCAGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGc    <  1:93517/149‑1 (MQ=255)
                                                                                                                                          cTTAACGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACTTCCTCTCCGATCTCACTACCGGcc   >  2:248756/1‑149 (MQ=255)
                                                                                                                                           ttAAGGCTGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGGCTTTTGCCTTCACGGCTGCGCTTGCAAGCAACACTCCCCTAGTGCCGAGATTAAGGACGGCTTCAACGTCCTCTCCGCTCTCAATACCGGccc  >  2:472225/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTCAGCCGAGGAAATTCAGGATCATCAGTGTGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGGCGGCCCCGCCGCGAGCGGTTTCGCCAATTCTAAGATCTTAAGATATGGGTCTTTTGCCTTCACGGCTGCGCTTGCCAGCAACACTCCCCTCGTGCCGAGTTTAAGCGCGGCTTCAACGTCCTCTCCGCTCTCAATACCCGCCC  >  NC_000913_3_pae_tpiA/4110601‑4110888

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