Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,530,252 C→A intergenic (‑9/+3) sgcC ← / ← sgcB putative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,530,2520CA88.6% 169.6 / 6.6 70intergenic (‑9/+3)sgcC/sgcBputative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB
Reads supporting (aligned to +/- strand):  ref base C (6/2);  new base A (41/21);  total (47/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAC  >  NC_000913/4530005‑4530498
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
ccTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAttt                                                                                                                                                                                                                                                      >  1:607342/1‑251 (MQ=255)
  tGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAtttt                                                                                                                                                                                                                                                     >  1:189702/1‑250 (MQ=255)
        aaCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGCTCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGAACGGATGTAATCAAACATAATGCCTCATTTTATTGAGt                                                                                                                                                                                                                                              >  1:403656/1‑251 (MQ=255)
        aaCGTAAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCACAACCGTTACCCATGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGCTAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGt                                                                                                                                                                                                                                              >  1:194302/1‑251 (MQ=255)
         aCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGtt                                                                                                                                                                                                                                             >  1:188481/1‑251 (MQ=255)
                  cctggaagggttgggcagAAACACTCAATCCTCATTTTGATGGGGGGACGGAGGCTGCCGATCCACATAACGAGACCCGTCTCCATACCGACAAAACCAATCCCAACGGTTACCCCGTCATTAAGCGCCGGAAGCCAGGGGAGGCGAAGAATAAGGCCAATCAGGAGTATGTTGATAGGGACGAAAACGGTATCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAATAGAGcc                                                                                                                                                                                                                                   <  2:514754/234‑1 (MQ=255)
                   cATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCt                                                                                                                                                                                                                                  >  1:91302/1‑251 (MQ=255)
                              gggCCAAAACGCTCAATCATCCTTTTGAGAGGGGAACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATCCCGACAAAACCAACCCCACCCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTtgttgt                                                                                                                                                                                                                       <  1:8603/250‑1 (MQ=255)
                                         cTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGCCAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTCTTGAGTTAACAGAGCCTTGGTTTGTTGTTTTAACTcgtc                                                                                                                                                                                                            >  2:131965/1‑251 (MQ=255)
                                          tCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAACGCCAATCACGATCATGATGATAGGGACGAAAACAGTACCGCCGAGAGACAGAATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTCAcgctgtg                                                                                                                                                                                                           >  2:120122/1‑247 (MQ=255)
                                              tcatcaCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGcgtcgtcgt                                                                                                                                                                                                       >  1:529630/1‑251 (MQ=255)
                                              tcatcaCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGcgtcgtcgt                                                                                                                                                                                                       >  1:30721/1‑251 (MQ=255)
                                               caacatttgtagggggggccagAGGCTGTCAACCCCCATAACAATCGCCAGTCCCATACCGCAAAACCAATACCCAACCTTTCCCCCGGCTTTAATCGCCGGAAGCCAGGGGATGCGGAAGATAGGGCCAATCACGATCATGATGAGAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGtt                                                                                                                                                                                                      <  2:191839/230‑1 (MQ=255)
                                                 tcaCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCTTCGTCGTTg                                                                                                                                                                                                     >  1:414105/1‑250 (MQ=255)
                                                     tttGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTATTTGAGTTAACAGAGCCTTGATTTGTGGTTTTAAAGCGTCGTCGTTGATccc                                                                                                                                                                                                >  2:310611/1‑251 (MQ=255)
                                                       tgaaaGGGGGACTGAGGCTGTCGATCGCCATGCCAATCACCAGTCCCATACCGACAAACCCAATCCCAACCGTTACCCCTGCCTTAACCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCAAACCGCTATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCt                                                                                                                                                                                               <  1:598758/246‑1 (MQ=255)
                                                          taagaggACTGAGGCTGGCGTTCGCCATGACGATCACCAGTCACATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCAGGAGGATAGGGACGACAACGGTACCGCCGATGGACAGGATGTAAACGAACAAAATGCATCATTATATTGAGTTAACAGGCACTTGGTTGGTGGTTTTAAAGCGTAGGCGTTGACCCCagtgt                                                                                                                                                                                           >  2:21247/6‑250 (MQ=255)
                                                                  cTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAAAATAATACATCATTTTATTGATTTAACTGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGGGGGCAGggcg                                                                                                                                                                                   >  2:141740/1‑251 (MQ=255)
                                                                    gggacTGTCGATCGACGTGCCGATCACCAGCTCTATCCGGACAAACCAAATCTCACCTGTTCCCACTGCCTTAAACGCTTGCAGCCGTGGGATGCGGAAGATAAGGCAAATCACGATCATGATGATAGGGACGAAAACGTAACCGCCGAGAGACAGGATGTAAGCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATGTGTTGTTTTAACGCTTCGTCGTTGATCCCTGTGAGCAGggcggc                                                                                                                                                                                 <  2:518334/247‑1 (MQ=255)
                                                                    gAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGggcggc                                                                                                                                                                                 >  1:484360/1‑251 (MQ=255)
                                                                          gTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCAt                                                                                                                                                                            <  1:216475/250‑1 (MQ=255)
                                                                            cGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTa                                                                                                                                                                          >  1:157130/1‑250 (MQ=255)
                                                                                   cATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAgg                                                                                                                                                                  >  1:62154/1‑251 (MQ=255)
                                                                                         gATCACCAGTCCCATACCGAAAAAAACAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAACAATAATGCCTCATTTTATTGAGATAACAGGGCATTGATTTGTTGTTTTAACGCGTCGTCGTTGATCACTGTGAACAGGGCGGCGCCATTAAggggtgggaagcc                                                                                                                                                            >  2:380386/1‑239 (MQ=255)
                                                                                           tCACCAGTCCCATCCCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAAGCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTgttgtt                                                                                                                                                                                                                      <  1:429414/191‑1 (MQ=255)
                                                                                           tCACCAGTACAATAACGACAAAACCAATCCCAAACGTTACACCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTgttgtt                                                                                                                                                                                                                      >  2:429414/1‑191 (MQ=255)
                                                                                               catccccATACCGACAAAACCAATCCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAACCCCGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCTCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATc                                                                                                                                                      <  1:179897/247‑1 (MQ=255)
                                                                                                   cccATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTg                                                                                                                                                  >  1:494341/1‑251 (MQ=255)
                                                                                                      cATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGtt                                                                                                                                                >  1:545170/1‑251 (MQ=255)
                                                                                                        tACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGAAGGGATGCCGTAATCACTGTTgg                                                                                                                                              >  2:623286/1‑251 (MQ=255)
                                                                                                              cAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGc                                                                                                                                         >  1:558659/1‑250 (MQ=255)
                                                                                                                     aaTCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACCGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGTTTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTACGCGTAGGGTTGCCGTAATCACTGTTGGTTCgactgacgat                                                                                                                                 >  2:482524/1‑242 (MQ=255)
                                                                                                                      aaccccAACCGTGACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAAAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATgg                                                                                                                                     <  1:482524/244‑1 (MQ=255)
                                                                                                                      aTCCCAACCGTTACCCCTGCCTTAATCGCCTCCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAACCGGTACCGCCGAGAGACAGGATGTAATCAAACATACTGCCTCATTTTATTGTGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTACGCGTAGGGATGCCGTACTCACTGTTGGTTCGCATGGAGGtt                                                                                                                                >  2:348164/1‑251 (MQ=255)
                                                                                                                         cccAACCGTTCCCCCTGCCTTAATCGCCAGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTAc                                                                                                                              <  1:387913/251‑1 (MQ=255)
                                                                                                                          ccAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCACACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTAc                                                                                                                              >  1:381915/1‑250 (MQ=255)
                                                                                                                            aacctttcccccTGCTTAAATCGCCTGCAGCCAGGGGATGCGGAAGAAAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTTCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAAt                                                                                                                           <  1:680142/243‑1 (MQ=255)
                                                                                                                              ccGTTACCCCTGCCATAACCGCCTGCAGCCAGGGGATGCGGAGGATAAGGCCAATCACGAGCAGGAGGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAAGCACTGTTGGTTCGCATGGAGGTTACAATaa                                                                                                                         <  1:302554/251‑1 (MQ=255)
                                                                                                                                           cTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGtt                                                                                                            >  1:568564/1‑251 (MQ=255)
                                                                                                                                             tAATCGCCTGCAGCCAGGGGCTGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGAACCGCCGAGAGAACGGATGTAATCAAACATAATGCCTCATATTATTGAGTGAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCCGGGCGGCGCCATTAAGCGTAGGGATGACGTAATCACTGTTGGGTCGCATGGAGGATACAATAAGGCCCATACCGTTAc                                                                                                          >  2:355785/1‑251 (MQ=255)
                                                                                                                                                   ccTGCAGCAAGGGGATGAGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGCGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACACAGCCTTGATTTGTTGTTTTAACGCGTCGGCGTTGATCCCTGTGCGCAGGTCGGCGCCATTAAGCGTAGGGCTCCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATACGGTCCATACCGTTACAGTTa                                                                                                     >  2:676107/1‑250 (MQ=255)
                                                                                                                                                     ccccccccgacgctcttccgatctAAGGCAAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTATGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAag                                                                                                   <  2:490193/228‑1 (MQ=255)
                                                                                                                                                        aGCCAGGGGATGCGGAAGATAAGGCCAATCACAATCATGATGATAAGGTTGAAAACGGTGCCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAggg                                                                                               <  1:635498/251‑1 (MQ=255)
                                                                                                                                                         gCCAGGGGAGGTAGAAGATAAGGCAAATTCAGTTCAAGATGATAGGGACGAAATCGGTACCGCCGAGAGACAGGATTTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGa                                                                                              <  2:549440/251‑1 (MQ=255)
                                                                                                                                                              gagaggCGGAAGAAAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGGACCGCAGAGAGACAGGATGAAATAAAACATAATGACTCATTTTATTGAGTTAAACGAGCCTTGATTTGTTGATTTAAAGCGTGGTCGTTGAACGATGTGAGCAGGGCGGCGCCATTAAACGAAGGGATGCCGAAATCACTGTTGGTTGGCATGGAGGTGACAATAAGGGCAATACCGtaccagttcaggggggtcgc                                                                                          >  2:479984/6‑231 (MQ=255)
                                                                                                                                                                     ggAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTACTTTTATTGAGGTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATACCTGGGAGCAGGGCGGCGCCATTCAGCGTAGGGATGCCGTCATCACTGTTGGATCGCATGGGGGTTACAATAAGGTCCATACCGTAACAGTTAAGAGGGATCTCATTCAGAc                                                                                  >  2:500899/1‑251 (MQ=255)
                                                                                                                                                                      gAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGAc                                                                                  >  1:225224/1‑250 (MQ=255)
                                                                                                                                                                             aaGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTATGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAag                                                                                                   >  1:490193/1‑227 (MQ=255)
                                                                                                                                                                             aaGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTg                                                                           >  1:60209/1‑251 (MQ=255)
                                                                                                                                                                                  cAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTGACAGTTAAGAGGGCTCTCATTACGACAGACTTGGGCgg                                                                      >  2:600500/1‑251 (MQ=255)
                                                                                                                                                                                   aaTCACGAACATGCTGATAGGATCGAAAACGGTACCGCAGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGCTCCATGTGTGTAGGGCGGCGCCTTTAAGAGTAGGGAGGCCGTAATAACTGTTTGATCGAGTGGAGGTTACAATAAGGTTCATACCGTTACCGTTAAGAGGGGtctgatgacgcaagcattgagcggg                                                                     >  2:637980/1‑229 (MQ=255)
                                                                                                                                                                                    aaccacgataaggagtatggggACGAAACCGGTGCCGCCGAGAGACAGGATGTAGTCAAACATAATGCCTCATTTTATTTAGTTAACAGTGCCTTGATTTGTTTTTTTAACGCGTCGTCGTTGGTCCCTGTTAGCAGGGCGGCTCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCTTACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCgg                                                                      <  1:637980/231‑1 (MQ=255)
                                                                                                                                                                                          gaacgtaatgatAGGGACGAAACCGTTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGc                                                               <  2:62154/244‑1 (MQ=255)
                                                                                                                                                                                          gATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGc                                                               >  1:231245/1‑251 (MQ=255)
                                                                                                                                                                                            tCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGAAATCAAACATAATGCCTCATTTTCTTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTg                                                             >  2:120549/1‑251 (MQ=255)
                                                                                                                                                                                                      gggaaaaaaacatgaactgccgAGAGTCAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGAAAGCATTGGGCGGTAGTGGCTGAAATGccgtg                                                   <  2:115541/233‑4 (MQ=255)
                                                                                                                                                                                                      aGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACCTAATGCCACATTTTGTGGAGTTAAACGAGCCGTGATTTGTTGTTTTAACGCGTAGTCGTTGATCACTGTGAGCAGGGCGGAGACATTAAGAGTAGGGATGCCGTCATAACAGTTGTTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGATAAGCGGGGTCTCATTCAGACAGCATTGGGCgggaggagctgcaataccctg                                                   >  2:154291/1‑232 (MQ=255)
                                                                                                                                                                                                        ggACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCCTTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACCATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGAAAGCATTGGGCGGTAGTGGCTGAAATGCCCTGtt                                                 >  2:528165/1‑251 (MQ=255)
                                                                                                                                                                                                                       cgccgAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCtt                                  >  1:545470/1‑251 (MQ=255)
                                                                                                                                                                                                                        tccgAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGCGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATAAAGCGTAGGGAGGCCATAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACATTTAAGAGGGATCTACTTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTg                                 >  2:313472/2‑251 (MQ=255)
                                                                                                                                                                                                                              agaCAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTAGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACCGTTAAGAGGGATCTCATTCAGACAGCATTGGGAGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGCACTCTTGCAGCtt                           >  1:197555/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATa                                                                                                                          >  2:118405/1‑150 (MQ=255)
                                                                                                                                                                                                                                    gATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATa                                                                                                                          <  1:118405/150‑1 (MQ=255)
                                                                                                                                                                                                                                       gTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCACCTTGTGTGCAAt                  <  1:122259/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGtt           <  1:78475/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aCATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATATCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGACCTGTTCAGTCAGGAACTCTTTAAGCTTGTGTGCAAACAGGGTTg          >  2:164496/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cATAATGCCTCATTTTATTGGGTTAACAGAGACTTGCTTTGTTGTTTTAACGCGTCGTCTTTGATCCCTGTGAGCAGAGCGGCGCCATTAAGCGTCGGGATCCCGTACTCACTGTTGGTTCGCATGGAGGTTACAATATGGTACCTACCGTTCCAGTTAAGGAGGACCTAATTCAGACACCATTGGACGGTAGTGGCTGTAATGCCCTGTTCAGTCAGGATATCTTGCAGATTGTGTGccatactgtatg          >  2:229796/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 aTAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGTATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGTGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcagtc                                            <  2:490944/215‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGTGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcagtc                                            >  1:490944/1‑215 (MQ=255)
                                                                                                                                                                                                                                                   aaTGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGt      <  1:533454/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       cctcatTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAc  <  1:446500/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
CCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAC  >  NC_000913/4530005‑4530498

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