Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,563,547 C→A 5.1% intergenic (+29/+240) 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,5470CA 5.1% 420.6 / 9.2 138intergenic (+29/+240)glpD/yzgLsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology
Reads supporting (aligned to +/- strand):  ref base C (59/72);  new base A (5/2);  total (64/74)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01

Only 100 of 139 total aligned reads displayed.

GCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTT  >  NC_000913/3563399‑3563696
                                                                                                                                                    |                                                                                                                                                     
gccgccGACGCCCTGTGGCGTCGCACAAACCAAGGCATGTGGCTAAATGCTGAGCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCa                                                                                                                                                      <  1:292194/150‑1 (MQ=255)
 ccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCaa                                                                                                                                                     >  1:676206/1‑150 (MQ=255)
 ccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGGCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCaa                                                                                                                                                     <  2:219633/150‑1 (MQ=255)
   gacgacGCCCTGTGGCGTCGCACAAAACCAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATc                                                                                                                                                   <  1:1105278/150‑1 (MQ=255)
   gacgacGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATc                                                                                                                                                   <  1:1165870/150‑1 (MQ=255)
      gacgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTgg                                                                                                                                                <  1:1732648/150‑1 (MQ=255)
       acgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGc                                                                                                                                               <  1:826791/150‑1 (MQ=255)
       acgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGc                                                                                                                                               >  2:181808/1‑150 (MQ=255)
         gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCt                                                                                                                                             <  1:576306/150‑1 (MQ=255)
          cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGGCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTg                                                                                                                                            <  1:372075/150‑1 (MQ=255)
            cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTgag                                                                                                                                          <  2:707037/150‑1 (MQ=255)
             ggtgGCGTCGCACAAACCAAGCCACGTGCCTAAATGCCGATCAACAATCTCTTCTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGCCAGATTGAAATCTGACCTgatc                                                                                                                                         <  2:253240/149‑3 (MQ=255)
              gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGAc                                                                                                                                        >  2:58129/1‑150 (MQ=255)
                 gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtga                                                                                                                                     >  1:619441/1‑150 (MQ=255)
                   gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                   <  1:19212/150‑1 (MQ=255)
                   gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCGGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatg                                                                                                                                   <  1:1682071/150‑1 (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)
                         caAACCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAac                                                                                                                             <  1:1734037/150‑1 (MQ=255)
                           aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAacac                                                                                                                           >  1:1589725/1‑150 (MQ=255)
                           aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTCGTGGAGTATACGCAGCAGAGGTTATCGCAGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGAATGATGACAAacac                                                                                                                           >  1:1191459/1‑150 (MQ=255)
                             aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACaa                                                                                                                         >  1:383096/1‑150 (MQ=255)
                                  gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTg                                                                                                                    >  2:681847/1‑150 (MQ=255)
                                    aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGcc                                                                                                                  >  2:1107333/1‑150 (MQ=255)
                                    aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGcc                                                                                                                  <  1:484946/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)
                                          cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCATCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgc                                                                                                            <  1:1349729/150‑1 (MQ=255)
                                          cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgc                                                                                                            <  2:618693/150‑1 (MQ=255)
                                          cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgc                                                                                                            <  1:493548/150‑1 (MQ=255)
                                            aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATgcgc                                                                                                          <  2:1476658/150‑1 (MQ=255)
                                            aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCGGAGACTGATGACAAACACAAAACTGCCTGATgcgc                                                                                                          >  2:284317/1‑150 (MQ=255)
                                            aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGATGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACACACACAAAACTGCCTGATgcgc                                                                                                          >  2:288040/1‑150 (MQ=255)
                                               tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTAc                                                                                                       <  2:206419/150‑1 (MQ=255)
                                               tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTAc                                                                                                       >  2:79355/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)
                                                   gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCGATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCtt                                                                                                   <  2:740976/150‑1 (MQ=255)
                                                   gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCtt                                                                                                   >  2:682370/1‑150 (MQ=255)
                                                      caacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATc                                                                                                >  2:937942/1‑150 (MQ=255)
                                                       atcaaGCTCGTGTGAGGCAGAGGCTGGTGGAGGATACCCCGCAGAGGTTTTCGCAGGCGTCGTAAATAACCGTTAGGTGATCAGGTCAGATTTCAATCTGGGCTGAGACTGCTGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                               <  2:1214619/148‑1 (MQ=25)
                                                       aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCa                                                                                               <  1:1787194/150‑1 (MQ=255)
                                                        acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAg                                                                                              <  1:701831/150‑1 (MQ=255)
                                                          aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGc                                                                                            <  2:1618820/150‑1 (MQ=255)
                                                               cGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                                       >  2:1383606/1‑150 (MQ=255)
                                                                   ggAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGt                                                                                   <  1:436350/149‑1 (MQ=255)
                                                                    gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGtt                                                                                  >  2:1137805/1‑150 (MQ=255)
                                                                    gAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGtt                                                                                  <  2:1200938/150‑1 (MQ=255)
                                                                      gTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTa                                                                                <  2:1656984/150‑1 (MQ=255)
                                                                         aGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGc                                                                             >  1:108560/1‑150 (MQ=255)
                                                                          gTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGTCGTCGTAAATTAACGTAATGTGATCAGGTCATATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCa                                                                            <  1:1504640/150‑1 (MQ=255)
                                                                              cTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGCCTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAatat                                                                        >  1:663388/1‑150 (MQ=255)
                                                                                 gtggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATAtt                                                                     <  2:1705366/150‑1 (MQ=255)
                                                                                  tggAGTATACGCAGCAGAGGTAATCCCTGGCGTCGTAAATTAACTTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGCCAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATAATg                                                                    >  2:967663/1‑150 (MQ=38)
                                                                                  gggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTg                                                                    <  2:1648519/149‑1 (MQ=255)
                                                                                      gTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCGGGCCTGCGACTGATGACATACCCAAAACTGCCTGATGCTCTTCGCTTATCAGGCCTACGTGGTTTATGCAAAATATTGAAtt                                                                >  1:276866/1‑150 (MQ=37)
                                                                                        ataCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTg                                                              >  2:1697564/1‑150 (MQ=255)
                                                                                        ataCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGGGGTTTATGCAATATATTGAATTTg                                                              >  2:663991/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)
                                                                                         taCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGc                                                             >  2:417212/1‑150 (MQ=255)
                                                                                          aCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                            >  1:936358/1‑150 (MQ=255)
                                                                                          aCGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCa                                                            <  2:619441/150‑1 (MQ=255)
                                                                                           cGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCAt                                                           >  2:1533077/1‑150 (MQ=255)
                                                                                           cGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCAt                                                           >  1:1843251/1‑150 (MQ=255)
                                                                                           cGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCAt                                                           <  1:1445215/150‑1 (MQ=255)
                                                                                            gcagcaGAGGTTATCGCTGGCGTCGTAAATTAACTTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATTTTGAATTTGCATg                                                          <  1:1349819/150‑1 (MQ=255)
                                                                                                cagAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCt                                                      >  2:434664/1‑150 (MQ=255)
                                                                                                    ggTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTa                                                  >  2:1011064/1‑150 (MQ=255)
                                                                                                     gTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAg                                                 <  1:405511/150‑1 (MQ=255)
                                                                                                        aTCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGcc                                              <  1:181808/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)
                                                                                                                    gtAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGtt                                  <  2:896733/150‑1 (MQ=255)
                                                                                                                     tAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTc                                 <  2:1799361/150‑1 (MQ=255)
                                                                                                                      aaaTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCa                                >  1:1077706/1‑150 (MQ=255)
                                                                                                                        aTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACg                              >  2:1278228/1‑150 (MQ=255)
                                                                                                                        aTTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACg                              >  1:52175/1‑150 (MQ=255)
                                                                                                                           aaCGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCg                           <  1:1391045/150‑1 (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:1289391/4‑150 (MQ=21)
                                                                                                                                    gtctCAGGCCAGATTGAAATCTGACCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCa                  >  2:1021455/4‑150 (MQ=21)
                                                                                                                                    gtctCAGGCCAGATTGAAATCTGACCTGAGACCGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCa                  >  1:1170/4‑150 (MQ=18)
                                                                                                                                      aTCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATg                <  2:1579305/150‑1 (MQ=255)
                                                                                                                                          ggTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACa            >  1:1685418/1‑150 (MQ=255)
                                                                                                                                           gTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACaa           <  1:937942/150‑1 (MQ=255)
                                                                                                                                               gATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAgcg       >  2:513751/1‑150 (MQ=255)
                                                                                                                                                aTTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAgcgc      <  2:1600539/150‑1 (MQ=255)
                                                                                                                                                 tttCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGACTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGATCACGTCGCATCCGGCATGAACAAAGCGCa     >  1:737774/1‑150 (MQ=255)
                                                                                                                                                  ttCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCAc    >  2:202845/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACAAACACAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTT  >  NC_000913/3563399‑3563696

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