Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,449 T→C 16.5% intergenic (+14/‑78) gpt → / → frsA xanthine phosphoribosyltransferase; xanthine‑guanine phosphoribosyltransferase/fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913256,4490TC16.5% 383.9 / 64.9 169intergenic (+14/‑78)gpt/frsAxanthine phosphoribosyltransferase; xanthine‑guanine phosphoribosyltransferase/fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity
Reads supporting (aligned to +/- strand):  ref base T (72/68);  new base C (16/12);  total (88/81)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 172 total aligned reads displayed.

TGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCG  >  NC_000913/256351‑256548
                                                                                                  |                                                                                                   
tgttgtTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTg                                                                                                    >  1:1474505/1‑100 (MQ=255)
  ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGc                                                                                                  <  2:1996468/100‑1 (MQ=255)
  ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGc                                                                                                  <  2:3421279/100‑1 (MQ=255)
    gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAc                                                                                                >  1:2870855/1‑100 (MQ=255)
    gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAc                                                                                                >  2:1008308/1‑100 (MQ=255)
    gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAc                                                                                                >  2:2046647/1‑100 (MQ=255)
             ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGcn                                                                                       <  1:3219623/100‑2 (MQ=255)
             ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCg                                                                                       >  1:1966652/1‑100 (MQ=255)
               gCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                     >  1:2255061/1‑100 (MQ=255)
               gCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                     >  2:2719150/1‑100 (MQ=255)
                cAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttg                                                                                    <  1:311224/100‑1 (MQ=255)
                 aaGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttga                                                                                   <  2:874217/100‑2 (MQ=39)
                  aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTACTCCCGCCAATCTCCGGTCGCTAATATTTTCAACGCCTGGCACTGCCGGGCgttg                                                                                    >  2:1090503/1‑98 (MQ=255)
                    aTACCTCGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCt                                                                                <  2:1641804/100‑1 (MQ=255)
                     tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                               >  1:1836384/1‑100 (MQ=255)
                     tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                               >  2:3519241/1‑100 (MQ=255)
                      aCCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttt                                                                              <  1:3551197/100‑1 (MQ=255)
                      aCCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttt                                                                              <  1:605986/100‑1 (MQ=255)
                      aCCTGGATTGAACAGCCGTCGGATATGGGGGTCGTATTCGTCCCGCCGAGCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                               <  2:3537169/99‑1 (MQ=255)
                       ccTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtttt                                                                             >  2:3008975/1‑100 (MQ=255)
                       ccTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgaaaagat                                                                             <  1:2927677/100‑8 (MQ=38)
                         tGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTa                                                                           <  2:234674/100‑1 (MQ=255)
                           gATTGACCAACCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATATTTTCAACGCCCGGCAGTGCCAGGCgttgaaaagatt                                                                            <  2:245719/97‑9 (MQ=25)
                           gATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAAc                                                                         <  2:112393/100‑1 (MQ=255)
                            aTTTAACAGCCGTCGGATATGGGTGTCGTATTCGTCCCGCCAATCTCCTGTCGCTAATCTTTTCAACCTCTTGCTCTGCCTGCCGTTGTTCTTTTTaa                                                                          <  2:1722868/98‑1 (MQ=255)
                            aTTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACt                                                                        <  1:1936860/100‑1 (MQ=255)
                            aTTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACt                                                                        >  1:562251/1‑100 (MQ=255)
                             ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACtt                                                                       >  1:2374604/1‑100 (MQ=255)
                               gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                     >  1:1539785/1‑100 (MQ=255)
                               gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                     >  1:3030119/1‑100 (MQ=255)
                                aaCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAg                                                                    >  2:1599875/1‑100 (MQ=255)
                                 aCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAgg                                                                   >  1:1486913/1‑100 (MQ=255)
                                   aGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAgg                                                                   >  2:3090269/1‑98 (MQ=255)
                                    gCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTaagatcggaaga                                                                >  2:905083/1‑90 (MQ=255)
                                    gCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                >  2:1199913/1‑100 (MQ=255)
                                      cGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                              <  1:2869114/100‑1 (MQ=255)
                                        tGGGATATGGGCGTCGTATTCGTCCCGCCCATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTa                                                            <  2:755676/100‑1 (MQ=255)
                                         gggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTAc                                                           <  1:2978359/100‑1 (MQ=255)
                                         gggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTAc                                                           <  1:3334508/100‑1 (MQ=255)
                                          ggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACa                                                          >  1:2991377/1‑100 (MQ=255)
                                          ggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACa                                                          <  1:1733424/100‑1 (MQ=255)
                                           gATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACaa                                                         <  1:1345778/100‑1 (MQ=255)
                                            atatGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAAt                                                        <  2:2199497/100‑1 (MQ=255)
                                                 ggCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGttt                                                   <  1:1008315/100‑1 (MQ=255)
                                                  gCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTc                                                  <  1:119708/100‑1 (MQ=255)
                                                    gtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTcc                                                 >  2:1789609/1‑99 (MQ=255)
                                                    gtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                >  2:1547805/1‑100 (MQ=255)
                                                    gtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                >  2:3325667/1‑100 (MQ=255)
                                                    gtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGACGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTcc                                                 <  2:268211/99‑1 (MQ=255)
                                                      cgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagt                                              >  2:3272727/1‑100 (MQ=255)
                                                        tATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaa                                            >  1:1384601/1‑100 (MQ=255)
                                                        tATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagta                                             <  2:3311451/99‑1 (MQ=255)
                                                         aTTCGTCCCGCCTATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                           <  1:1723713/100‑1 (MQ=255)
                                                         aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                           >  1:772361/1‑100 (MQ=255)
                                                         aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                           >  2:2814108/1‑100 (MQ=255)
                                                           tCGTCCCGACAATCTCTGGTCGCTAATCTTTTCAACGTCTGGCACTGCCGGGCGTTGTTCTTTTTAACGTCAGGCGGGTTACAATAGTTTCCagtaagta                                         <  1:404213/100‑1 (MQ=255)
                                                            cGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAt                                        >  2:81049/1‑100 (MQ=255)
                                                            cGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAt                                        <  1:596867/100‑1 (MQ=255)
                                                            cGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAt                                        >  1:731996/1‑100 (MQ=255)
                                                             gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTGAGGCGGGTTACAATAGTTTCCAGTAAGTAtt                                       >  2:2381349/1‑100 (MQ=255)
                                                             gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAtt                                       >  1:2505690/1‑100 (MQ=255)
                                                              tCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTc                                      <  1:3063133/100‑1 (MQ=255)
                                                                ccGCCTATCTCCGGTCGCTAATCTTTTCACCTTCTGGCACTGCCGGTCGTTGTTCTTTTTAATTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCt                                     <  1:347700/99‑1 (MQ=255)
                                                                ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                    <  1:81049/100‑1 (MQ=255)
                                                                ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                    >  2:3266704/1‑100 (MQ=255)
                                                                 cGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTgg                                   >  2:3214265/1‑100 (MQ=255)
                                                                    cAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAgg                                >  2:2625192/1‑100 (MQ=255)
                                                                     cttccgatctCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGc                               <  1:454909/91‑1 (MQ=255)
                                                                       tctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTg                             <  2:1836375/100‑1 (MQ=255)
                                                                       tctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTg                             >  2:1326954/1‑100 (MQ=255)
                                                                       tctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTg                             <  1:2698590/100‑1 (MQ=255)
                                                                        ctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGc                            <  1:1258692/100‑1 (MQ=255)
                                                                               cGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGcagatcggaag                     >  2:454907/1‑90 (MQ=255)
                                                                               cGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATg                     <  2:1379586/100‑1 (MQ=255)
                                                                                 cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                   >  1:1120246/1‑100 (MQ=255)
                                                                                 cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                   >  2:811494/1‑100 (MQ=255)
                                                                                 cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                   >  1:1611606/1‑100 (MQ=255)
                                                                                   aaaaagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGacac                 >  1:2729291/9‑100 (MQ=37)
                                                                                    aaaagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAGCTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGAcaca                >  2:1985571/8‑100 (MQ=37)
                                                                                    aTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGAcaca                >  1:2977371/1‑100 (MQ=255)
                                                                                      cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg              <  2:2743921/100‑1 (MQ=255)
                                                                                       ttttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg              >  2:1473032/1‑99 (MQ=255)
                                                                                        tttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCa            <  1:20845/100‑1 (MQ=255)
                                                                                        tttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCa            <  1:1031085/100‑1 (MQ=255)
                                                                                        gaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCa            >  1:316392/4‑100 (MQ=38)
                                                                                          tCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa          <  1:3640385/100‑1 (MQ=255)
                                                                                          acAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa          >  1:1690048/2‑100 (MQ=39)
                                                                                          acAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa          >  1:990647/2‑100 (MQ=39)
                                                                                          acAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa          >  2:2344007/2‑100 (MQ=39)
                                                                                          acAACACCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGTCTACAATAGTTTCCCGTAAGTATTCTGGAGTCTGCATCCATGACACAGGc             >  2:638881/2‑97 (MQ=255)
                                                                                           cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCCGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCCAAc         >  1:3382404/1‑100 (MQ=255)
                                                                                           cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAc         >  1:3516932/1‑100 (MQ=255)
                                                                                           cAACGCCTGGCACTGCCGGGAGTTGTTCTTTTTAACTTCAGGCGGATTCAAAAATTTTCCAATAAGTAGTATGCAGGATGCCTCCATGaccacgacagac         >  2:3620235/1‑90 (MQ=255)
                                                                                             aCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCt       >  2:1845059/1‑100 (MQ=255)
                                                                                                ccTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAg    <  2:479154/100‑1 (MQ=255)
                                                                                                 cTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGc   <  1:1849908/100‑1 (MQ=255)
                                                                                                 cTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGc   <  2:25554/100‑1 (MQ=255)
                                                                                                  tGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCg  <  1:2814128/100‑1 (MQ=255)
                                                                                                  tGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCg  <  1:3405785/100‑1 (MQ=255)
                                                                                                  tGGCACTGCCGGGCGTTGTTCGTTTTAACCTCAGGCAGGGTACAGCAGGTGCGAGTGAGTATTCTGGCGGCTGCCTCCCGCGCACAGGCAGAACAGAGCg  >  2:1298176/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCG  >  NC_000913/256351‑256548

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