Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,563,554 G→A 6.4% intergenic (+36/+170) glpD → / ← yzgL aerobic glycerol 3‑phosphate dehydrogenase/putative uncharacterized protein YzgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,563,5540GA 6.4% 243.3 / 4.8 78intergenic (+36/+170)glpD/yzgLaerobic glycerol 3‑phosphate dehydrogenase/putative uncharacterized protein YzgL
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base A (2/3);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAA  >  NC_000913/3563409‑3563702
                                                                                                                                                 |                                                                                                                                                    
cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTg                                                                                                                                                  >  1:1030430/1‑150 (MQ=255)
   tgtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTGAAATCTGACCTgatc                                                                                                                                               <  1:865081/150‑3 (MQ=255)
   tgtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTGAAATCTGACCTgatc                                                                                                                                               <  1:870315/150‑3 (MQ=255)
         gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                         <  2:755149/150‑1 (MQ=255)
           cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGAc                                                                                                                                       <  2:87135/150‑1 (MQ=255)
              acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACaaa                                                                                                                                    >  1:1010535/1‑150 (MQ=255)
                    cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaaa                                                                                                                              >  1:263229/1‑150 (MQ=255)
                      aGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAAc                                                                                                                            <  1:1166465/150‑1 (MQ=255)
                        gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGGTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTg                                                                                                                          <  2:997963/150‑1 (MQ=255)
                        gCATGTGGCTAAATGCGGATCAACAATCACGTGTCAGTCAGTGGCTGGTGGAGTATACGCACAAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACtt                                                                                                                          >  2:498265/1‑149 (MQ=255)
                            gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTg                                                                                                                      >  2:955609/1‑150 (MQ=255)
                              ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGAt                                                                                                                    >  1:449719/1‑150 (MQ=255)
                                cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgc                                                                                                                  >  1:734545/1‑150 (MQ=255)
                                    aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTa                                                                                                              <  1:1074781/150‑1 (MQ=255)
                                    aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTa                                                                                                              <  1:723245/150‑1 (MQ=255)
                                        ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCt                                                                                                          >  1:1260043/1‑150 (MQ=255)
                                        ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCATGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCt                                                                                                          <  1:1002278/150‑1 (MQ=255)
                                             aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                     >  1:1003797/1‑150 (MQ=255)
                                             aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                     <  1:174087/150‑1 (MQ=255)
                                             aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                     >  2:1059733/1‑150 (MQ=255)
                                              acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAg                                                                                                    >  1:673547/1‑150 (MQ=255)
                                               caaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                   <  2:734915/150‑1 (MQ=255)
                                                aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGc                                                                                                  <  1:1167019/150‑1 (MQ=255)
                                                 aTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGcc                                                                                                 <  2:59522/150‑1 (MQ=255)
                                                 aTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGcc                                                                                                 >  1:286813/1‑150 (MQ=255)
                                                     ggtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTTATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                                             <  2:1098391/149‑1 (MQ=255)
                                                     cGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                                             >  1:1086468/1‑150 (MQ=255)
                                                      gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGTCTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                                            <  1:1059733/150‑1 (MQ=255)
                                                      gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                                            >  2:1051011/1‑150 (MQ=255)
                                                          gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGtt                                                                                        >  1:319907/1‑150 (MQ=255)
                                                           aGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGttt                                                                                       <  2:721771/150‑1 (MQ=255)
                                                           aGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGttt                                                                                       >  1:182537/1‑150 (MQ=255)
                                                            gTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                      >  1:12193/1‑150 (MQ=255)
                                                              cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATg                                                                                    <  2:286813/150‑1 (MQ=255)
                                                                gTGGCTGGTGGAGTATACGCAGCAGAGGTTAACGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                                  >  2:1369034/1‑150 (MQ=255)
                                                                  ggCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAat                                                                                >  1:499840/1‑150 (MQ=255)
                                                                   gCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAata                                                                               <  1:1283336/150‑1 (MQ=255)
                                                                       gtggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATAtt                                                                           >  1:310734/1‑150 (MQ=255)
                                                                           aGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAAt                                                                       <  2:321618/150‑1 (MQ=255)
                                                                                aCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                                  <  2:285350/150‑1 (MQ=255)
                                                                                 cGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCAt                                                                 >  2:458464/1‑150 (MQ=255)
                                                                                  gcagcaGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATg                                                                <  1:1369151/150‑1 (MQ=255)
                                                                                  gcagcaGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATg                                                                <  1:1352993/150‑1 (MQ=255)
                                                                                          ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                        <  1:1321078/150‑1 (MQ=255)
                                                                                          ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                        >  1:824107/1‑150 (MQ=255)
                                                                                          ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                        <  2:1165588/150‑1 (MQ=255)
                                                                                          ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                        >  1:199414/1‑150 (MQ=255)
                                                                                           gTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAg                                                       <  1:619782/150‑1 (MQ=255)
                                                                                               tCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                                            >  2:80492/1‑109 (MQ=255)
                                                                                               tCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                                            <  1:80492/109‑1 (MQ=255)
                                                                                               tCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCa                                                   >  2:1310722/1‑150 (MQ=255)
                                                                                                       gtcgtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTAtg                                                                                                      <  1:1103879/91‑2 (MQ=255)
                                                                                                       gtcgtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGAc                                           <  2:449719/150‑1 (MQ=255)
                                                                                                       gtcgtAAATTAACGTAAGGTCATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATc                                                                                                      >  2:1103879/1‑91 (MQ=255)
                                                                                                          gtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACAAAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGtt                                        >  1:1079355/1‑150 (MQ=255)
                                                                                                           tAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTc                                       >  1:1288541/1‑150 (MQ=255)
                                                                                                            aaaTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCa                                      >  2:562253/1‑150 (MQ=255)
                                                                                                             aaTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCATGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCAc                                     >  2:1214707/1‑150 (MQ=255)
                                                                                                              aTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACg                                    <  2:12193/150‑1 (MQ=255)
                                                                                                                 aaCGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCg                                 <  2:990524/150‑1 (MQ=255)
                                                                                                                    gTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCAt                              >  2:1380334/1‑150 (MQ=255)
                                                                                                                        ggTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCgg                          <  1:816646/150‑1 (MQ=255)
                                                                                                                         gTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGc                         >  2:1058107/1‑150 (MQ=255)
                                                                                                                         gTGATCAGGTCAGATTTCAATATGGACTCAGACAGATGACAAACACAAAACTGCCTGATGCGCTAAGCTTATCAGGAATACCTGGTTTATGAAAAATATTGAATTTGCATGGTCTTGTAGGACAGATAAGACGTTCACGTCGCAtaccgc                         >  2:1041542/1‑145 (MQ=255)
                                                                                                                          tGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCa                        >  1:1199603/1‑150 (MQ=255)
                                                                                                                             tCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGa                     <  1:955609/150‑1 (MQ=255)
                                                                                                                              cAGGCCAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTc                                       <  2:379728/131‑1 (MQ=255)
                                                                                                                              cAGGCCAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTc                                       >  1:379728/1‑131 (MQ=255)
                                                                                                                               aGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAAc                   <  2:1086468/150‑1 (MQ=255)
                                                                                                                               aGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAAc                   <  1:458464/150‑1 (MQ=255)
                                                                                                                                ggccAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACa                  >  2:132028/4‑150 (MQ=255)
                                                                                                                                ggTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACa                  >  1:220690/1‑150 (MQ=255)
                                                                                                                                   cAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAg               >  1:1384085/1‑150 (MQ=255)
                                                                                                                                      aTTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAgcgc            <  2:1199603/150‑1 (MQ=255)
                                                                                                                                            aTCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTg      >  1:1010557/1‑150 (MQ=255)
                                                                                                                                             tCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGt     >  2:339034/1‑150 (MQ=255)
                                                                                                                                             tCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGAATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGt     >  2:344225/1‑150 (MQ=255)
                                                                                                                                              cTGGCCTGAGACTGATGACAAAGACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTc    <  1:1041542/150‑1 (MQ=255)
                                                                                                                                                ggCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCaa  <  2:1010535/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAA  >  NC_000913/3563409‑3563702

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: