Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 298,093 A→G 5.5% S212P (TCA→CCA)  paoD ← moco insertion factor for PaoABC aldehyde oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,0930AG 5.5% 221.4 / 3.3 72S212P (TCA→CCA) paoDmoco insertion factor for PaoABC aldehyde oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (28/40);  new base G (2/2);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

GACGTAACGTGTGGGTTCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTGAGATTCCAGCCCGT  >  NC_000913/297946‑298240
                                                                                                                                                   |                                                                                                                                                   
gACGTAACGTGTGGGTTCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGa                                                                                                                                                     <  8:96738/148‑1 (MQ=255)
 aCGTAACGTGTGGGTTCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAg                                                                                                                                                    <  3:184052/148‑1 (MQ=255)
       cGTGTGGGTTCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCgg                                                                                                                                               >  2:66221/1‑147 (MQ=255)
           tgGGTTCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGaaaa                                                                                                                                           <  2:102245/147‑1 (MQ=255)
            gggTTCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAg                                                                                                                                         <  4:328702/148‑1 (MQ=255)
             ggTTCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGa                                                                                                                                        <  5:193863/148‑1 (MQ=255)
             ggTTCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGa                                                                                                                                        <  4:171688/148‑1 (MQ=255)
              gTTCGATAGCTGCCCAATGCGCCGAGCTAACAGGGTTTTGCGTCTCGCGCGGCCTGCAACACTGGCAGCGCCCTGGTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGGATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGAt                                                                                                                                       <  3:78225/148‑1 (MQ=255)
               ttCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCCACACTGGCAGCTCCCGGTTGAGCTCATGGCACAGCAAAATGACCGCCGTATCGGTATCGAGCTGCGCGCTGGCTGAGGCCGGaaa                                                                                                                                            >  3:44866/1‑142 (MQ=255)
                  gATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGa                                                                                                                                    <  5:309554/147‑1 (MQ=255)
                  gATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCACGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGa                                                                                                                                    <  5:132694/147‑1 (MQ=255)
                     aGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGat                                                                                                                                <  5:29541/148‑1 (MQ=255)
                      gCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGata                                                                                                                               >  5:210675/1‑148 (MQ=255)
                      gCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGat                                                                                                                                <  3:127090/147‑1 (MQ=255)
                      gCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAAGGACCGCCGTAGCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAACAGATCGAAGat                                                                                                                                <  6:326792/147‑1 (MQ=255)
                        tGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATg                                                                                                                             >  6:63618/1‑148 (MQ=255)
                          cccAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATg                                                                                                                             <  1:293059/146‑1 (MQ=255)
                               tGCGCCGAGATTATAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCa                                                                                                                      <  4:61706/148‑1 (MQ=255)
                                        aTAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCtgtg                                                                                                             <  2:235821/148‑1 (MQ=255)
                                         tAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTgg                                                                                                            >  3:15424/1‑148 (MQ=255)
                                                 ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCTGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCa                                                                                                    <  4:44866/148‑1 (MQ=255)
                                                 ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCtgctgc                                                                                                     >  5:332891/1‑147 (MQ=255)
                                                 ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCtgctgc                                                                                                     >  7:257034/1‑147 (MQ=255)
                                                            gcgcgGCCTGCAACACTGGCATCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTcgctg                                                                                          <  3:190658/147‑3 (MQ=255)
                                                            gcgcgGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCTCAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTcgcg                                                                                           <  2:97293/146‑1 (MQ=255)
                                                                         cacTGGCAGCTCCCGGTTGAGATCATTGCACATCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCa                                                                            <  7:10951/148‑1 (MQ=255)
                                                                         cacTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCa                                                                            <  5:298287/148‑1 (MQ=255)
                                                                         cacTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCa                                                                            <  5:63618/148‑1 (MQ=255)
                                                                           cTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGCCCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAag                                                                          <  4:15424/148‑1 (MQ=255)
                                                                               cAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGaa                                                                      >  5:376637/1‑148 (MQ=255)
                                                                               cAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGaa                                                                      >  4:5395/1‑148 (MQ=255)
                                                                               cAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGaa                                                                      <  1:23171/148‑1 (MQ=255)
                                                                                 gCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACg                                                                    >  2:335703/1‑148 (MQ=255)
                                                                                    cccGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAGGAGCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTcc                                                                 >  1:258016/1‑148 (MQ=255)
                                                                                    cccGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTcc                                                                 >  8:184612/1‑148 (MQ=255)
                                                                                     ccGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTc                                                                  <  1:66221/146‑1 (MQ=255)
                                                                                      cGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGGGGCCGGAAAAAGa                                                                                                                                        >  7:280478/1‑75 (MQ=255)
                                                                                      cGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGGGGCCGGAAAAAGa                                                                                                                                        >  7:285693/1‑75 (MQ=255)
                                                                                      cGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGGGGCCGGAAAAAGa                                                                                                                                        <  8:280478/75‑1 (MQ=255)
                                                                                      cGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGGGGCCGGAAAAAGa                                                                                                                                        <  8:285693/75‑1 (MQ=255)
                                                                                             gaTCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAatcat                                                        <  7:1993/148‑1 (MQ=255)
                                                                                             gaTCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAatcat                                                        <  7:203996/148‑1 (MQ=255)
                                                                                              aTCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAatcatc                                                       <  2:199871/148‑1 (MQ=255)
                                                                                               tCATGGCACAGCAAAATGACCGCCTTATCGGTATCGATCTGAGCGCTGGCTTAGGCCTGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGTCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAatcat                                                        >  5:222185/1‑146 (MQ=255)
                                                                                               tCATGGCACAGCAAAATGACCGCCGTATCGGTATCGAGGTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATAGGGCTGTCATCGGGGGTGGCTGCTGCAAGGCTCGCGGTTGCGTGGGGGTGAAGAGGACGTGCGTAGatcgtgc                                                      >  2:156515/1‑144 (MQ=255)
                                                                                               tCATGGCACAGCAAAATGACCGCCGTATCGGTATCGAGCTGAGCGCTGGATGAGGACGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCAGGAAatcatc                                                       <  6:222185/147‑1 (MQ=255)
                                                                                                 aTGGCACACCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAAtcatcac                                                     <  2:260499/147‑2 (MQ=255)
                                                                                                   ggCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCt                                                  >  1:140865/1‑148 (MQ=255)
                                                                                                        aGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACg                                              >  1:163644/1‑147 (MQ=255)
                                                                                                         gCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGTCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCa                                            <  7:253443/148‑1 (MQ=255)
                                                                                                          cAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCAt                                           <  2:234245/148‑1 (MQ=255)
                                                                                                          cAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCAt                                           <  6:203373/148‑1 (MQ=255)
                                                                                                                    gccgTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAAcc                                 >  7:37462/1‑148 (MQ=255)
                                                                                                                         atcggtatcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTgcg                                                                                 >  4:180204/1‑95 (MQ=255)
                                                                                                                         atcggtatcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTgcg                                                                                 <  3:180204/95‑1 (MQ=255)
                                                                                                                            ggtatcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTc                          >  7:62801/1‑147 (MQ=255)
                                                                                                                             gtatcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCaa                        <  5:3262/148‑1 (MQ=255)
                                                                                                                             gtatcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCa                         >  2:278621/1‑147 (MQ=255)
                                                                                                                              tatcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCaaa                       >  8:111488/1‑148 (MQ=255)
                                                                                                                              tatcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCaaa                       <  6:376637/148‑1 (MQ=255)
                                                                                                                              tatcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCa                         >  7:241529/1‑146 (MQ=255)
                                                                                                                              tatcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCa                         >  3:164960/1‑146 (MQ=255)
                                                                                                                                tcgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGc                     <  2:140865/148‑1 (MQ=255)
                                                                                                                                 cgATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGc                     >  7:355191/1‑147 (MQ=255)
                                                                                                                                   aTCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCAt                  >  6:234736/1‑148 (MQ=255)
                                                                                                                                        aGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTgag              <  6:332891/147‑1 (MQ=255)
                                                                                                                                        aGCGCTGGATGAGGACGGAATAAGATCGAAGATATGGCTGTCATAGCCTTTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTgag              <  2:163644/147‑1 (MQ=255)
                                                                                                                                         gcgcTGGCTGAGGCCGGAAAAAGATCGATGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCCGCCTGCCGCATGGCCTGCCCCCCACCTCAAAGCCATTGAGAt            >  1:315635/1‑148 (MQ=255)
                                                                                                                                          cgcTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTGAGAtt           >  7:135132/1‑148 (MQ=255)
                                                                                                                                               gctgAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTGAGATTCCAg       <  3:5395/147‑1 (MQ=255)
                                                                                                                                                ctgAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTGAGATTCCAGcc     <  1:278621/148‑1 (MQ=255)
                                                                                                                                                ctgAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTGAGATTCCAGcc     >  4:339744/1‑148 (MQ=255)
                                                                                                                                                  gAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTGAGATTCCAGCCCg   <  8:355191/148‑1 (MQ=255)
                                                                                                                                                   aGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTGAGATTCCAGCCCGt  <  7:111488/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GACGTAACGTGTGGGTTCGATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCTGCTGCAAGACTCGCGGTTGCCTGCGCCTCAAGAGAACGTCCGTAAATCATCAGCCTGACGCATGGCCTGAACCCCACCTCAAAGCCATTGAGATTCCAGCCCGT  >  NC_000913/297946‑298240

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