Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,235,485 T→C 6.2% intergenic (+79/‑420) pgi → / → yjbE glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,235,4850TC 6.2% 182.7 / 3.4 65intergenic (+79/‑420)pgi/yjbEglucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (26/35);  new base C (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

TAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAA  >  NC_000913/4235337‑4235632
                                                                                                                                                    |                                                                                                                                                   
tAAAGACATCAGCGGCCACGATAGCTCGACCAATGGTATGATTAACCGCTATAAAGCGTGGCGCGGTTAATCGTCGGCGCTGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                                                                                                     <  2:126767/149‑1 (MQ=37)
  aaGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTGACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTGTTAt                                                                                                                                                   >  2:44767/1‑149 (MQ=255)
         cagcagCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcat                                                                                                                                            <  1:161040/149‑3 (MQ=255)
         cagcagCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcat                                                                                                                                            <  1:211374/149‑3 (MQ=255)
         cagcagCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcat                                                                                                                                            <  1:223488/149‑3 (MQ=38)
            cagcCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAGGCCTGATAAGACGCGACCGCGTCGCATCAGGcatcgg                                                                                                                                         <  2:255755/149‑6 (MQ=37)
               ccACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                                                                                                      <  2:286757/149‑1 (MQ=255)
               ccACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                                                                                                      <  2:318484/149‑1 (MQ=255)
                cACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                                                                                                     >  1:301219/1‑149 (MQ=255)
                     tAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGCCGCGTCTTATCAGGCCTACAGGTCGATGc                                                                                                                                <  1:476420/149‑1 (MQ=255)
                       gCTCGACCAATGGTCGGATTAACCGCGATAAAGCGTGGCGCGGTCCGGGCGCTTGGACGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCg                                                                                                                              <  1:343358/149‑1 (MQ=21)
                           gACCAATTGTCTGATTAACCGCTATAACGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                                                                                                          <  1:211286/149‑1 (MQ=255)
                           gACCAATGGTCTGATTAACCGCTATAAGGCGTGGCGTCGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGTTCGATGCCGattt                                                                                                                          <  2:748156/149‑3 (MQ=255)
                            aCCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGCCGATGCCGatagg                                                                                                                         >  1:244010/1‑147 (MQ=255)
                                  ggTCGGATCAACCGTGAGAAGGCGTGGCGAGGTTAATCAAAGAGGATATATAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAGGCGTGAGGCGACGCGGTCGCGTCTTATCAGGCTTACAGGTCGATGCCGATATGTACATc                                                                                                                   <  2:512064/149‑1 (MQ=21)
                                    tCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGt                                                                                                                 >  2:543359/1‑149 (MQ=255)
                                    tCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGt                                                                                                                 >  2:550557/1‑149 (MQ=255)
                                       gATTAACCGCTATAAAGCCTGGCGCGGTTAACCATCGTCGACACGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                                                                                                              <  1:320646/149‑1 (MQ=37)
                                        aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                                                                                                             >  2:322389/1‑149 (MQ=255)
                                           aaCCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                                                                                                          <  2:261615/149‑1 (MQ=255)
                                            aCCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                                                                                                         <  2:592423/149‑1 (MQ=255)
                                            aCCGCGATAAAGCGGGGCGCGGTTAAGCATCGTCGATATGTGGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                                                                                                         <  2:72956/149‑1 (MQ=255)
                                             ccGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                                                                                                        >  1:49056/1‑149 (MQ=255)
                                                   taAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTCCAGGTCGATGCCGATATGTACATCGTATTCGGCACTTAATa                                                                                                  >  1:144994/1‑149 (MQ=255)
                                                   taAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGCTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGGCGGGTCTTATCAGGCATAAAGGTCGATGCCGATAGGTACATGGTAGTCGGCAATTAATa                                                                                                  >  2:510964/1‑149 (MQ=37)
                                                    aaaGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc                                                                                                 <  1:90303/149‑1 (MQ=255)
                                                    aaaGCGGGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGGGTCTTATCAGGCCTACAGGTCGATGCCGAGATGTACATCGTATTGGGCAATTCAtgc                                                                                                 >  2:586433/1‑147 (MQ=37)
                                                       gcggtgcgcGGTTAACCAGCGGCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCGTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa                                                                                              <  2:127392/144‑1 (MQ=25)
                                                       gCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa                                                                                              <  2:725513/149‑1 (MQ=255)
                                                       gCGGGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACTTCGTATTCGGCTATTAATACATa                                                                                              <  2:656558/149‑1 (MQ=255)
                                                              gcgGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat                                                                                       >  2:52251/1‑149 (MQ=255)
                                                                         tcgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATaa                                                                            >  1:135528/1‑149 (MQ=255)
                                                                         tcgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAACTACCACGATTGATTAAATaa                                                                            >  2:179329/1‑149 (MQ=255)
                                                                         tcgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAAATCGTATTCGGCAATTAATACATAGCCGGATTGATTAAATaa                                                                            >  2:9282/1‑149 (MQ=255)
                                                                          ggtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAc                                                                           <  2:660912/148‑1 (MQ=255)
                                                                            tcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATACCACGATTGATTAAATAACCt                                                                         >  1:636743/1‑149 (MQ=255)
                                                                               atatGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTtaa                                                                      <  1:400243/149‑1 (MQ=255)
                                                                                    tAGGCAGGATCAGGCGGTGGCGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACa                                                                 >  2:289377/1‑149 (MQ=38)
                                                                                        ccGGATAAGGCGTTCACGCGGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGc                                                             <  1:195169/149‑1 (MQ=255)
                                                                                                 gCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTAt                                                    >  2:658057/1‑149 (MQ=255)
                                                                                                 gCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAAGTTAATAACACTGCGGAGGTTAt                                                    >  2:214723/1‑149 (MQ=38)
                                                                                                  cGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAATCTTAATAACAATGCCGACGTTATg                                                   <  1:583598/149‑1 (MQ=255)
                                                                                                  cGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATg                                                   <  1:231556/149‑1 (MQ=255)
                                                                                                  cGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATg                                                   <  1:582756/149‑1 (MQ=255)
                                                                                                   gTTGACGCCGCAGCCGGCAACCCATGCCTGACGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAGCTGATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGt                                                  <  2:497382/149‑1 (MQ=38)
                                                                                                   gTTCACGCCGCAGCCGGCAACCGATGCCTGACGCGACGCGGTCGCGTTTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGt                                                  <  2:236878/149‑1 (MQ=255)
                                                                                                    ttCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATAAGGCATACAAGGGAAGGCCGAAATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGGCGTTCTGTc                                                 >  2:285070/1‑149 (MQ=255)
                                                                                                      cACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCgg                                               >  1:550589/1‑149 (MQ=255)
                                                                                                         gccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATCACAATGCCGACGTTATGTCGGcct                                            >  1:279124/1‑147 (MQ=255)
                                                                                                              aTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAtttttt                                       >  2:734540/1‑149 (MQ=255)
                                                                                                              aTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATCCCGACGTTATATCGGCAtttttt                                       >  2:486829/1‑149 (MQ=255)
                                                                                                              aTCCAGCAACCGATGCCTGATGCGACGCAGTCGCGGCTTATCAGGCCTACAGGTCGATGCCGGTATATACATCGTATTCGGCAATTAATACATAGCACGATTGTTTAAACAACCTTAATAACCATGCCGCCGTTATGTCGGCAtgtttt                                       >  2:604740/1‑149 (MQ=255)
                                                                                                                ccGGCAAACGATGCCTGATGCGACGCGGTCGCCTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTa                                     <  1:44767/149‑1 (MQ=255)
                                                                                                                     gaCCGATGCCTGATGCGACGCGGGCGCGGCGTATCAGGCCTACAGGTCGATGGCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAAGGCCGACGTTATGTCGGCATTTTTTTATCAGa                                <  2:99060/148‑1 (MQ=255)
                                                                                                                        cGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGAGATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATACCGACATTATGTCGGCATTTTTTTATCAGATaa                             >  1:445526/1‑149 (MQ=255)
                                                                                                                         taggCCTGATAAGACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATaaa                            >  1:605918/4‑149 (MQ=38)
                                                                                                                         taggCCTGATAAGACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATaaa                            >  1:564412/4‑149 (MQ=38)
                                                                                                                              cTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATcccc                       <  1:543359/149‑1 (MQ=255)
                                                                                                                              cTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTACGTCGGCATTTTTTTATCAGATAAATcccc                       <  1:550557/149‑1 (MQ=255)
                                                                                                                              cTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGCACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATcccc                       <  1:227361/149‑1 (MQ=255)
                                                                                                                                        gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACCTCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACCATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTaa             <  2:50713/149‑1 (MQ=255)
                                                                                                                                         cgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACCTCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGAAGATATGGCGGCATTTTTTTATCAGATAAATCTCCTTGTTTGTAAt            >  2:377899/1‑149 (MQ=255)
                                                                                                                                          ggTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAAtt           <  2:660009/149‑1 (MQ=255)
                                                                                                                                             cgcgTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTaa        <  2:66941/149‑1 (MQ=255)
                                                                                                                                               tgTCTTATCAGGCCTACCGGTCGAGGCCGACATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATCACCTTAATAACAATGCCGACGTTATGCCGGCATTTTTTTACCAGATAAATCCCCTTGTCTGTAATTTAACg      <  1:156984/148‑1 (MQ=255)
                                                                                                                                                  cTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGaa   <  2:301219/149‑1 (MQ=255)
                                                                                                                                                   ttATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTATCTGTAATTTAAAGGaaa  >  2:254773/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAA  >  NC_000913/4235337‑4235632

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: