Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,563,539 T→C 5.6% intergenic (+21/+248) glpD → / ← yzgL sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,563,5390TC 5.6% 415.3 / 11.1 142intergenic (+21/+248)glpD/yzgLsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology
Reads supporting (aligned to +/- strand):  ref base T (58/76);  new base C (5/3);  total (63/79)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

Only 100 of 142 total aligned reads displayed.

CCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAA  >  NC_000913/3563392‑3563687
                                                                                                                                                   |                                                                                                                                                    
ccgccgCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACTATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCa                                                                                                                                                    <  1:38182/150‑1 (MQ=255)
 cgccgcgccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAg                                                                                                                                                   >  1:740976/1‑150 (MQ=255)
 cgccgcgccgACGACGCCCTGGGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAg                                                                                                                                                   <  2:298857/150‑1 (MQ=255)
  gccgcgccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGAGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTCGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGa                                                                                                                                                  >  2:289947/1‑150 (MQ=255)
   ccgcgccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGGCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAt                                                                                                                                                 <  2:1305831/150‑1 (MQ=255)
    cgcgccgACGACGCCCTGTGGCGTATCACAAAACAAGTCATGTGGCTAAATGCGGATAAACAATCTCGTGTGAGTCAGTGTCTTGTTGAGTATACGCATCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAtt                                                                                                                                                <  2:1270914/150‑1 (MQ=255)
     gcgccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAttt                                                                                                                                               <  2:1219381/150‑1 (MQ=255)
     gcgccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGAttt                                                                                                                                               >  1:896733/1‑150 (MQ=255)
        ccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCaa                                                                                                                                            >  1:676206/1‑150 (MQ=255)
        ccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGGCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCaa                                                                                                                                            <  2:219633/150‑1 (MQ=255)
         cgacgaCGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAAt                                                                                                                                           >  2:1305571/1‑150 (MQ=255)
         cgacgaCGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAAt                                                                                                                                           >  1:1705366/1‑150 (MQ=255)
          gacgacGCCCTGTGGCGTCGCACAAAACCAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATc                                                                                                                                          <  1:1105278/150‑1 (MQ=255)
          gacgacGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATc                                                                                                                                          <  1:1165870/150‑1 (MQ=255)
             gacgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTgg                                                                                                                                       <  1:1732648/150‑1 (MQ=255)
              acgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGc                                                                                                                                      >  2:181808/1‑150 (MQ=255)
               cgCCCTGTGTCGTCGCACAAAACAAGGCAGGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGcc                                                                                                                                     <  1:83354/150‑1 (MQ=255)
                gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCt                                                                                                                                    <  2:621657/150‑1 (MQ=255)
                 cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTg                                                                                                                                   >  1:1271809/1‑150 (MQ=255)
                 cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGGCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTg                                                                                                                                   <  1:372075/150‑1 (MQ=255)
                   cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTgag                                                                                                                                 <  1:1772009/150‑1 (MQ=255)
                   cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTgag                                                                                                                                 <  2:707037/150‑1 (MQ=255)
                        gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtga                                                                                                                            >  1:619441/1‑150 (MQ=255)
                          gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                          <  1:19212/150‑1 (MQ=255)
                           tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatga                                                                                                                         <  1:673082/150‑1 (MQ=255)
                            cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGAc                                                                                                                        >  2:1491997/1‑150 (MQ=255)
                             gCACAAACCACGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACa                                                                                                                       <  1:91451/150‑1 (MQ=255)
                             gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACa                                                                                                                       >  1:1579305/1‑150 (MQ=255)
                             gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTCTGAGGCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAATTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGGCTGAGACtgaggtca                                                                                                                       >  1:632221/1‑145 (MQ=255)
                               acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACaaa                                                                                                                     <  1:1278017/150‑1 (MQ=255)
                                caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAac                                                                                                                    <  1:51365/150‑1 (MQ=255)
                                    aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaa                                                                                                                >  1:383096/1‑150 (MQ=255)
                                         gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTg                                                                                                           >  2:681847/1‑150 (MQ=255)
                                           aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGcc                                                                                                         <  1:484946/150‑1 (MQ=255)
                                           aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGcc                                                                                                         <  1:580267/150‑1 (MQ=255)
                                           aTGTGGCTAAATGCGGAGCAACAATCTCGTGTGAGTCAGTGGCGGTTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTACCGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGcc                                                                                                         <  2:617927/150‑1 (MQ=255)
                                               ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTTGCCTGAGACTGATGACAAACACAAAACTGCCTGAt                                                                                                     >  1:1602598/1‑150 (MQ=255)
                                               ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGAt                                                                                                     >  2:260379/1‑150 (MQ=255)
                                                gCTAAATGCGGATCAACAATCTCGTGTGAGTCCGTGGCTCGGGGAGTATACGCAGCAGACGGTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATg                                                                                                    <  2:273062/150‑1 (MQ=255)
                                                 cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCATCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgc                                                                                                   <  1:1349729/150‑1 (MQ=255)
                                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcgc                                                                                                 >  2:21040/1‑150 (MQ=255)
                                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcgc                                                                                                 <  2:1476658/150‑1 (MQ=255)
                                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCGGAGACTGATGACAAACACAAAACTGCCTGATgcgc                                                                                                 >  2:284317/1‑150 (MQ=255)
                                                      tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTAc                                                                                              >  2:79355/1‑150 (MQ=255)
                                                        cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGc                                                                                            >  2:1611710/1‑150 (MQ=255)
                                                        cGGAGCCACAATCTCGAGTGAGGCAGGGGCTGGCGGAGTATACCCAGCAGAGGCTCTCGCGGGCGTCGTAAATTAACGTAAGGTTATCAGGTCAGATTTCAATCTGGACTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGc                                                                                            <  2:1151761/150‑1 (MQ=37)
                                                         ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCt                                                                                           <  1:319529/150‑1 (MQ=255)
                                                         ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCt                                                                                           >  2:59764/1‑150 (MQ=255)
                                                          gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCtt                                                                                          >  2:682370/1‑150 (MQ=255)
                                                             caacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATc                                                                                       >  2:937942/1‑150 (MQ=255)
                                                              atcaaGCTCGTGTGAGGCAGAGGCTGGTGGAGGATACCCCGCAGAGGTTTTCGCAGGCGTCGTAAATAACCGTTAGGTGATCAGGTCAGATTTCAATCTGGGCTGAGACTGCTGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                      <  2:1214619/148‑1 (MQ=25)
                                                                caaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                    <  2:573048/150‑1 (MQ=255)
                                                                 aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGc                                                                                   <  2:1618820/150‑1 (MQ=255)
                                                                    ctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTa                                                                                >  2:729074/1‑150 (MQ=255)
                                                                          ggAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGt                                                                          <  1:436350/149‑1 (MQ=255)
                                                                           gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGtt                                                                         >  2:1137805/1‑150 (MQ=255)
                                                                             gTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                       <  2:1656984/150‑1 (MQ=255)
                                                                               cAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATg                                                                     <  2:1691779/150‑1 (MQ=255)
                                                                                aGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGc                                                                    >  1:108560/1‑150 (MQ=255)
                                                                                 gTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGTCGTCGTAAATTAACGTAATGTGATCAGGTCATATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                   <  1:1504640/150‑1 (MQ=255)
                                                                                       ggtggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAatatat                                                             <  2:618014/150‑1 (MQ=255)
                                                                                        gtggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATAtt                                                            <  2:1705366/150‑1 (MQ=255)
                                                                                         tggAGTATACGCAGCAGAGGTAATCCCTGGCGTCGTAAATTAACTTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGCCAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATAATg                                                           >  2:967663/1‑150 (MQ=38)
                                                                                             gTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATTCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAAtt                                                       >  1:1391620/1‑150 (MQ=255)
                                                                                             gTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCGGGCCTGCGACTGATGACATACCCAAAACTGCCTGATGCTCTTCGCTTATCAGGCCTACGTGGTTTATGCAAAATATTGAAtt                                                       >  1:276866/1‑150 (MQ=37)
                                                                                               ataCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTg                                                     >  1:415999/1‑150 (MQ=255)
                                                                                                taCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCGGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGc                                                    <  2:498223/150‑1 (MQ=255)
                                                                                                taCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGc                                                    >  1:1302573/1‑150 (MQ=255)
                                                                                                taCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGc                                                    <  2:1744509/150‑1 (MQ=255)
                                                                                                 aCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                   <  2:619441/150‑1 (MQ=255)
                                                                                                 aCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                   >  2:1041177/1‑150 (MQ=255)
                                                                                                  cGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCAt                                                  >  1:1843251/1‑150 (MQ=255)
                                                                                                   gcagcaGAGGTTATCGCTGGCGTCGTAAATTAACTTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATTTTGAATTTGCATg                                                 <  1:1349819/150‑1 (MQ=255)
                                                                                                       cagAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACTCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCt                                             <  1:289947/150‑1 (MQ=255)
                                                                                                       cagAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCt                                             >  2:434664/1‑150 (MQ=255)
                                                                                                               aTCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGcc                                     <  1:181808/150‑1 (MQ=255)
                                                                                                                 cgcgGGCGACGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAg                                   <  2:190323/150‑1 (MQ=255)
                                                                                                                     ggCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTTGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATaa                               <  1:1611710/150‑1 (MQ=255)
                                                                                                                     ggCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATaa                               >  2:422645/1‑150 (MQ=255)
                                                                                                                     ggCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATaa                               >  2:455444/1‑150 (MQ=255)
                                                                                                                         tcgtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACg                           <  1:1749703/150‑1 (MQ=255)
                                                                                                                         tcgtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACg                           <  2:116053/150‑1 (MQ=255)
                                                                                                                           gtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGtt                         <  2:896733/150‑1 (MQ=255)
                                                                                                                             aaaTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCa                       >  1:1077706/1‑150 (MQ=255)
                                                                                                                               aTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACg                     >  1:52175/1‑150 (MQ=255)
                                                                                                                               aTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACg                     >  2:1278228/1‑150 (MQ=255)
                                                                                                                               aTTAACGTAAGGTGATCAGGTAAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCGTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACg                     <  1:1571000/150‑1 (MQ=255)
                                                                                                                                  aaCGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCg                  <  1:1391045/150‑1 (MQ=255)
                                                                                                                                  aaCGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCg                  >  2:1723761/1‑150 (MQ=255)
                                                                                                                                    cGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCa                <  2:1271809/150‑1 (MQ=255)
                                                                                                                                        tcagtctCAGGCCAGATTGAAATCTGACCTGAGTCTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAg                              >  1:725198/7‑132 (MQ=18)
                                                                                                                                        tcagtctCAGGCCAGATTGAAATCTGACCTGAGTCTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAg                              <  2:725198/126‑1 (MQ=18)
                                                                                                                                          gTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGc          <  2:378297/150‑1 (MQ=255)
                                                                                                                                          agtctCAGGCCAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGc          >  2:415148/5‑150 (MQ=21)
                                                                                                                                           gtctCAGGCCAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCa         >  2:1021455/4‑150 (MQ=21)
                                                                                                                                           gtctCAGGCCAGATTGAAATCTGACCTGAGACCGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCa         >  1:1170/4‑150 (MQ=18)
                                                                                                                                             aTCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATg       <  2:1007617/150‑1 (MQ=255)
                                                                                                                                             aTCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATg       <  2:1579305/150‑1 (MQ=255)
                                                                                                                                                 ggTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACa   >  1:1685418/1‑150 (MQ=255)
                                                                                                                                                  gTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACaa  <  1:937942/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAA  >  NC_000913/3563392‑3563687

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