Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,491,274:1 +C 6.5% coding (677/1695 nt) oxc ← oxalyl‑CoA decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,491,2741.C 6.5% 378.2 / 21.2 93coding (677/1695 nt)oxcoxalyl‑CoA decarboxylase
Reads supporting (aligned to +/- strand):  ref base . (51/36);  new base C (4/2);  total (55/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTCTTTT  >  NC_000913/2491128‑2491423
                                                                                                                                                   |                                                                                                                                                     
tGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTt‑gtg                                                                                                                                                    <  1:55634/150‑1 (MQ=255)
tGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTt‑gtg                                                                                                                                                    <  2:409581/150‑1 (MQ=255)
   cAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAAt                                                                                                                                                 >  2:978990/1‑150 (MQ=255)
     aCGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAc                                                                                                                                               <  1:1270753/150‑1 (MQ=255)
        cgttcgtttGCCCTGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgcc                                                                                                                                            >  1:1001951/7‑150 (MQ=255)
        cgttcgtttGCCCTGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgcc                                                                                                                                            >  1:31026/7‑150 (MQ=255)
         gttcgtttGCCCTGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgccg                                                                                                                                           >  2:1532077/6‑150 (MQ=255)
         cAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATAcgccg                                                                                                                                           >  2:1314087/1‑150 (MQ=255)
         cAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttattg                                                                                                                                           <  2:353682/150‑7 (MQ=38)
         cAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCAC‑AAGCTGATGAACAGCTTCGTGAAtttattg                                                                                                                                           <  2:419229/150‑7 (MQ=38)
           gCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCcgcg                                                                                                                                         >  2:1237606/1‑150 (MQ=255)
           gCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCcgcg                                                                                                                                         <  1:82055/150‑1 (MQ=255)
                tGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCttt                                                                                                                                    >  1:593809/1‑150 (MQ=255)
                 gCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTg                                                                                                                                   <  2:1263709/150‑1 (MQ=255)
                 gCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTg                                                                                                                                   <  2:373270/150‑1 (MQ=255)
                  ccAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGc                                                                                                                                  >  1:1411995/1‑150 (MQ=255)
                   cAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGcc                                                                                                                                 <  2:300601/150‑1 (MQ=255)
                     gggCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCTATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCaa                                                                                                                               >  2:181861/1‑150 (MQ=255)
                       gCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAgg                                                                                                                             <  2:1179256/150‑1 (MQ=255)
                         aaaCGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAg                                                                                                                              >  2:634561/1‑147 (MQ=255)
                         aaaCGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAg                                                                                                                              <  1:634561/147‑1 (MQ=255)
                          cacgaacgCTCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGCCATTGGAAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGAGCCTTTGCCAAGGAta                                                                                                                          >  2:649690/2‑150 (MQ=255)
                          aacgaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAta                                                                                                                          >  2:1393420/1‑150 (MQ=255)
                            cgaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                                        >  2:1128609/1‑150 (MQ=255)
                            cgaacgCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGAtaat                                                                                                                        <  2:303958/150‑1 (MQ=255)
                                 gcgcAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATg                                                                                                                   <  2:1311147/150‑1 (MQ=255)
                                  cgcAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATgg                                                                                                                  >  2:1450900/1‑150 (MQ=255)
                                    cAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGcc                                                                                                                >  2:200338/1‑150 (MQ=255)
                                       cTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGtt                                                                                                             >  1:1608699/1‑150 (MQ=255)
                                          ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAg                                                                                                          >  2:512213/1‑150 (MQ=255)
                                          ccGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAg                                                                                                          <  2:566492/150‑1 (MQ=255)
                                             cAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCtt                                                                                                       >  1:131197/1‑150 (MQ=255)
                                                  aGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCt                                                                                                  <  2:1056065/150‑1 (MQ=255)
                                                   gTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGc                                                                                                         <  2:692941/142‑1 (MQ=255)
                                                   gTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGc                                                                                                         >  1:692941/1‑142 (MQ=255)
                                                   gTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTa                                                                                                 <  1:489381/150‑1 (MQ=255)
                                                    tGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACATTCAATAAATTCACGAAGCTGTTCATCAGCGT‑GGGAATACGCCGCGCCTTTGCCAAGGCTAATTAAGGGCCGTGCAGCTTTTGCTaa                                                                                                >  2:329366/1‑150 (MQ=255)
                                                      gATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTaaaa                                                                                              >  1:421541/1‑150 (MQ=255)
                                                       aTGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAg                                                                                             <  2:1280218/150‑1 (MQ=255)
                                                         gCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCg                                                                                           >  2:52193/1‑150 (MQ=255)
                                                          cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                                                                          <  1:1237606/150‑1 (MQ=255)
                                                          cGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGa                                                                                          >  1:332261/1‑150 (MQ=255)
                                                              tCTTCAAGGATCCCTTTCGCCATAGACATTGGAAGGAATGGAATCTGAGCACTTTCAATAAATTCAAGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAAtt                                                                                      >  1:1467513/1‑150 (MQ=255)
                                                               cTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCATTTGCCAAGGATAATTAATGGCCGTTAAGCTTTTGCTAAAAGCGAAATTg                                                                                     >  1:24251/1‑150 (MQ=255)
                                                                ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgc                                                                                    >  2:1116276/1‑150 (MQ=255)
                                                                ttCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgc                                                                                    >  1:1331417/1‑150 (MQ=255)
                                                                 tCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTgcg                                                                                   <  1:1532077/150‑1 (MQ=255)
                                                                        tCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtga                                                                            >  2:1379708/1‑150 (MQ=255)
                                                                        tCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtga                                                                            >  2:530236/1‑150 (MQ=255)
                                                                           cTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactg                                                                         >  2:1532073/1‑150 (MQ=255)
                                                                             ttCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGtgactgac                                                                       >  1:107395/1‑150 (MQ=255)
                                                                                   aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCggg                                                                 >  2:1609113/1‑150 (MQ=255)
                                                                                   aTAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCggg                                                                 >  2:1115915/1‑150 (MQ=255)
                                                                                    tAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGc                                                                <  2:1428582/150‑1 (MQ=255)
                                                                                    tAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGc                                                                >  1:498072/1‑150 (MQ=255)
                                                                                     aGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCa                                                               <  2:107395/150‑1 (MQ=255)
                                                                                      gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                            <  2:853441/120‑1 (MQ=255)
                                                                                      gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGc                                                                                            >  1:853441/1‑120 (MQ=255)
                                                                                      gACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCAt                                                              >  1:1415332/1‑150 (MQ=255)
                                                                                       aCATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATg                                                             <  2:972853/150‑1 (MQ=255)
                                                                                         aTTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGc                                                           <  2:360259/150‑1 (MQ=255)
                                                                                         aTTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGc                                                           <  1:1128609/150‑1 (MQ=255)
                                                                                          ttGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCa                                                          <  2:421541/150‑1 (MQ=255)
                                                                                          ttGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCa                                                          <  1:542718/150‑1 (MQ=255)
                                                                                              cAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCaataa                                                      >  2:1052747/1‑150 (MQ=255)
                                                                                                       aatCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGa                                             >  2:1412230/1‑150 (MQ=255)
                                                                                                            gAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGAt                                        >  2:996130/1‑150 (MQ=255)
                                                                                                              gCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGAttt                                      <  2:827602/150‑1 (MQ=255)
                                                                                                                aCTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTc                                    <  1:1052747/150‑1 (MQ=255)
                                                                                                                   ttCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAAc                                 <  1:1532073/150‑1 (MQ=255)
                                                                                                                    tCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACt                                >  1:377456/1‑150 (MQ=255)
                                                                                                                    tCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACt                                <  1:1393420/150‑1 (MQ=255)
                                                                                                                      ggcgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACttt                              >  1:749553/6‑150 (MQ=255)
                                                                                                                      ggcgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACttt                              >  2:1172389/6‑150 (MQ=255)
                                                                                                                      aaTAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACttt                              >  2:342037/1‑150 (MQ=255)
                                                                                                                       gcgtaTTCAC‑AAGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTa                             >  1:317047/5‑150 (MQ=255)
                                                                                                                        tAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTaa                            >  1:988604/1‑150 (MQ=255)
                                                                                                                         aaaTTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAAc                           <  1:181861/150‑1 (MQ=255)
                                                                                                                          aaTTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACa                          <  2:1608699/150‑1 (MQ=255)
                                                                                                                              cACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCg                      >  2:959058/1‑150 (MQ=255)
                                                                                                                               aCGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGt                     <  2:1497213/150‑1 (MQ=255)
                                                                                                                                cGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg                    >  1:554755/1‑150 (MQ=255)
                                                                                                                                cGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg                    >  2:1448306/1‑150 (MQ=255)
                                                                                                                                cGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg                    >  2:278372/1‑150 (MQ=255)
                                                                                                                                acaaGCTGATGAACAGCTTCGTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGt                     >  1:986489/3‑150 (MQ=255)
                                                                                                                                acaaGCTGATGAACAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTg                    >  2:453704/3‑150 (MQ=255)
                                                                                                                                 gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTgg                   <  1:959058/150‑1 (MQ=255)
                                                                                                                                 gAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTgg                   <  2:1432602/150‑1 (MQ=255)
                                                                                                                                   aGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGtt                 <  2:332261/150‑1 (MQ=255)
                                                                                                                                      tGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAAc              <  1:329366/150‑1 (MQ=255)
                                                                                                                                          catcaGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCtt          >  1:204836/1‑150 (MQ=255)
                                                                                                                                            tcaGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCg        >  2:137402/1‑150 (MQ=255)
                                                                                                                                             caGCTT‑GTGAATACGCCGCGCCTTTGCCTAGGATAAGTAATGGCCGTTCAGCTTTTGCTAAACGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGt       >  2:1509658/1‑150 (MQ=255)
                                                                                                                                              aGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTc      >  2:730667/1‑150 (MQ=255)
                                                                                                                                               gCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTCt     <  1:1609113/150‑1 (MQ=255)
                                                                                                                                                  t‑gtgAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTgg                   >  2:1053920/1‑133 (MQ=255)
                                                                                                                                                  t‑gtgAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTgg                   <  1:1053920/133‑1 (MQ=255)
                                                                                                                                                  t‑gtgAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTCtttt  >  1:519283/1‑150 (MQ=255)
                                                                                                                                                  t‑gtgAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTCtttt  <  1:1379708/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGACAACGTCAGCATTTGCCAGGGCAAACGAACGCGCAGCTGCCGCAGAAAGTGGATGCGTATCTTCAAGGATCCCTTTCGCCATAGACATTGGCAGGAATGGAATCTGAGCACTTTCAATAAATTCACGAAGCTGTTCATCAGCTT‑GTGAATACGCCGCGCCTTTGCCAAGGATAATTAATGGCCGTTCAGCTTTTGCTAAAAGCGAAATTGCGCTAGTGACTGACTTCGGGCATGGCAATAATGCTGGCGACGGATTTTCAACTTTAACAATCGTGGTTAACGCTTCGTCTTTT  >  NC_000913/2491128‑2491423

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