Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,863,778 C→G 11.0% S62S (TCC→TCG) ‡ ygbL → putative class II aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,863,7780CG11.0% 289.5 / 22.2 102S62S (TCC→TCG) ‡ygbLputative class II aldolase
Reads supporting (aligned to +/- strand):  ref base C (37/53);  new base G (5/6);  total (42/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

Only 100 of 103 total aligned reads displayed.

GGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAG  >  NC_000913/2863631‑2863926
                                                                                                                                                   |                                                                                                                                                    
ggagATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGTCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCaa                                                                                                                                                    >  1:1100437/1‑150 (MQ=255)
ggagATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCaa                                                                                                                                                    >  2:1387046/1‑150 (MQ=255)
 gagaTGACGCGGATTGCCTGTTCATTCATACAGCGCGGCTATGCAACCGGTTGGGGTGGCAATCTGTCGGTGCTTTTACCTGACGGGTATTTACAGGCGACAGCGACAGGTACATGGCTGGGCTATCTCGATCCGCCGCCGCTGTCCtaa                                                                                                                                                   >  1:1060424/1‑150 (MQ=255)
 gagaTGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCaaa                                                                                                                                                   <  2:961425/150‑1 (MQ=255)
    aTGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTc                                                                                                                                                <  2:586362/150‑1 (MQ=255)
     tGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcg                                                                                                                                               <  2:46122/150‑1 (MQ=255)
      gACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                              >  1:878836/1‑150 (MQ=255)
      gACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                              >  1:1731933/1‑150 (MQ=255)
         gcgGATTGCCAGTTCATTCTTTCAGCTCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGAAACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgccgc                                                                                                                                           <  2:821073/150‑1 (MQ=255)
         gcgGATTCCCAGTTCCTTCTTGCAGCGCGGCTATGCACCCGGCGCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGCGAATTTACTGGCGACACCGACAGGGTCATGCCTGGGCAATCTCGATCCGCGGCGACTTTGGAAAGCcgctgc                                                                                                                                           <  1:272301/150‑1 (MQ=255)
          cgGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCcgcg                                                                                                                                          <  1:878874/150‑1 (MQ=255)
          cgGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGACTTTGGAAAGCCGCtgcg                                                                                                                                          <  2:46433/150‑1 (MQ=255)
             aTTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCGGCGACTTTGGAAAGCCGCTGCGGAt                                                                                                                                       <  1:967134/150‑1 (MQ=255)
             aTTGCCAGTTCATGCTTTCAGCGCGGCTATGCAACCGGCTCGGCTGGCAATCTGTCGCTGCTTTTCCCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCGGCGACTTTGGAAAGCCGCTGCGGAt                                                                                                                                       <  1:87761/150‑1 (MQ=255)
              ttGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATg                                                                                                                                      <  1:468866/150‑1 (MQ=255)
              ttGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATg                                                                                                                                      <  1:1127592/150‑1 (MQ=255)
               tGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCGGCGACTTTGGAAAGCCGCTGCGGAtcg                                                                                                                                     <  2:629288/150‑3 (MQ=255)
                gCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGc                                                                                                                                    >  1:1250820/1‑150 (MQ=255)
                gCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGc                                                                                                                                    >  2:1033196/1‑150 (MQ=255)
                 ccAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCg                                                                                                                                   >  2:175582/1‑150 (MQ=255)
                      gcatccattCAGCACGGCTATCCAACGGGATCGGCTGGCAATCTGTCGCTTCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGATCAAGCCTGGGCAATCTAGATCCGGAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATgg                                                                                                                              <  2:1705535/143‑1 (MQ=255)
                       cATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGt                                                                                                                             >  1:1567068/1‑150 (MQ=255)
                         ttCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTa                                                                                                                           <  2:438070/150‑1 (MQ=255)
                           cTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAg                                                                                                                         >  1:1374005/1‑150 (MQ=255)
                           cTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAg                                                                                                                         >  1:1024518/1‑150 (MQ=255)
                           cTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAg                                                                                                                         >  2:67300/1‑150 (MQ=255)
                           cTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAg                                                                                                                         >  2:627748/1‑150 (MQ=255)
                             ttCAGCGCGGCTATGCAACCGGTTCTGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtg                                                                                                                       <  2:11218/150‑1 (MQ=255)
                             ttCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGGCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtg                                                                                                                       <  1:1519285/150‑1 (MQ=255)
                                 gccgggCTATGCAACAGGTTCGGCTGGCAATCTGTCGCTGCTTTAACCTGACGGGAATATACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGa                                                                                                                   <  2:169156/146‑1 (MQ=255)
                                  cgcgGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGAc                                                                                                                  <  2:410465/150‑1 (MQ=255)
                                   gcgGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACa                                                                                                                 <  1:844514/150‑1 (MQ=255)
                                      gCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAAc                                                                                                              >  2:1508626/1‑150 (MQ=255)
                                        tATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGCCACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAAccc                                                                                                            >  2:1440480/1‑150 (MQ=255)
                                          tGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTc                                                                                                          <  2:70898/150‑1 (MQ=255)
                                           gCAACCGGTTCGTCTGGCAATCGGGCGCTGCGTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCg                                                                                                         <  2:1079995/150‑1 (MQ=255)
                                           gCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCg                                                                                                         >  2:1678159/1‑150 (MQ=255)
                                           gCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCg                                                                                                         >  1:1306380/1‑150 (MQ=255)
                                           gCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCg                                                                                                         <  2:1468289/150‑1 (MQ=255)
                                            cAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGa                                                                                                        >  2:366742/1‑150 (MQ=255)
                                             aaCCGTTTCTGCTGGCAATCTGTCGCTGCTTTTACCTGACGTGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGACGCGGATGGCGAATGGTTAAGTGGTGACAATCACTCGaa                                                                                                       <  2:386144/150‑1 (MQ=255)
                                              aCCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGaaa                                                                                                      <  1:1102247/150‑1 (MQ=255)
                                                cGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAaga                                                                                                    <  1:175582/150‑1 (MQ=255)
                                                 ggTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAagag                                                                                                   >  2:805384/1‑150 (MQ=255)
                                                  gTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAgg                                                                                                  <  2:606382/150‑1 (MQ=255)
                                                         gggCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATTCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCtt                                                                                           <  2:1082186/149‑1 (MQ=255)
                                                          ggCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGCCGCCGCGGATGGCGAATGGTTAAGTGGCGACAAACCCTCGAAAGAGGTGCTCttt                                                                                          <  2:159744/150‑1 (MQ=255)
                                                          ggCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCAGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGAAGGCGAATGGTTAAGTGGTGACACACCCTCGAAAGAGGtccacatt                                                                                          >  2:251518/1‑143 (MQ=255)
                                                            cAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCa                                                                                        >  1:644525/1‑150 (MQ=255)
                                                            cAATCTGGCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCa                                                                                        <  1:1633785/150‑1 (MQ=255)
                                                              aTCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCAAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATc                                                                                      <  2:577903/150‑1 (MQ=255)
                                                               tCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGCACCTCTTTcgagg                                                                                     <  2:631820/150‑5 (MQ=255)
                                                                 tGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTgg                                                                                   <  2:1154491/150‑1 (MQ=255)
                                                                  gTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGgc                                                                                  <  2:352795/150‑1 (MQ=255)
                                                                      ctgctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTg                                                                              <  1:170314/150‑1 (MQ=255)
                                                                      ctgctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTg                                                                              <  2:1291577/150‑1 (MQ=255)
                                                                           tttACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCg                                                                         <  2:1374005/150‑1 (MQ=255)
                                                                           tttACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCg                                                                         >  1:776525/1‑150 (MQ=255)
                                                                            ttACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGc                                                                        >  1:308730/1‑150 (MQ=255)
                                                                               ccTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaac                                                                     <  2:1351764/150‑1 (MQ=255)
                                                                                  gACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAAt                                                                  <  1:1703954/150‑1 (MQ=255)
                                                                                    cGGGAATTTACTGGCGACACCGTCAGGTTCATGCCTGGGCAATCTCGCTCCTCAGCGGCTATCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGGGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATcc                                                                <  2:1629180/150‑1 (MQ=255)
                                                                                      ggAATTTACTGGCTACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCcgc                                                              <  2:90352/150‑1 (MQ=255)
                                                                                       gAATTTACTCGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCGTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCcgcg                                                             <  2:1067605/150‑1 (MQ=255)
                                                                                           ttACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGCTCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGAGAATGGTTACGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGt                                                         <  1:67300/150‑1 (MQ=255)
                                                                                             aCTGGCGCCACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGtta                                                       >  2:1090388/1‑148 (MQ=255)
                                                                                              cTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTaaa                                                      >  1:724893/1‑150 (MQ=255)
                                                                                                ggCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGc                                                    <  1:753961/150‑1 (MQ=255)
                                                                                                     cacCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtgg                                               <  2:1250820/150‑1 (MQ=255)
                                                                                                     cacCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtgg                                               <  1:17383/150‑1 (MQ=255)
                                                                                                       ccGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGgtggtg                                             <  1:347915/150‑1 (MQ=255)
                                                                                                         gACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCa                                           <  2:521116/150‑1 (MQ=255)
                                                                                                           cAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCAtt                                         <  2:1100437/150‑1 (MQ=255)
                                                                                                           cAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCAtt                                         <  2:480733/150‑1 (MQ=255)
                                                                                                           cAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCAtt                                         >  2:1209573/1‑150 (MQ=255)
                                                                                                             ggTTCATGCCTGTGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGATAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTg                                       <  1:979336/150‑1 (MQ=255)
                                                                                                                 cATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcaca                                   <  1:1387046/150‑1 (MQ=255)
                                                                                                                 cATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGAGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcaca                                   >  2:933746/1‑150 (MQ=255)
                                                                                                                   tggcTGGGCAATCTCGATCCGCAGCGGCTTTCCATAGTCGCCGGGGATGGCGTTTGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTGATATGGCGCTTTTTAGAAAGAATCCGCGCTGTAATGGTGTGTTGCATTTgcaaagt                                 >  1:110764/4‑149 (MQ=255)
                                                                                                                   tGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcacagc                                 <  2:324052/150‑1 (MQ=255)
                                                                                                                     ccTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcacagcac                               <  1:1131966/150‑1 (MQ=255)
                                                                                                                      cTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcacagcaca                              >  2:272406/1‑150 (MQ=255)
                                                                                                                         ggCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATgg                           <  1:1033196/150‑1 (MQ=255)
                                                                                                                         ggCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATgg                           <  1:1440480/150‑1 (MQ=255)
                                                                                                                            aaTCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCg                        >  1:1319760/1‑150 (MQ=255)
                                                                                                                             ttcgccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCACCAATCCGCGCTGTAAAGCGGTGGTGCTTTTGACCAGGACATGGTCGa                       >  2:1137960/7‑150 (MQ=37)
                                                                                                                              tctcGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGAc                      <  1:1508626/150‑1 (MQ=255)
                                                                                                                               cgccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACg                     >  1:935518/5‑150 (MQ=39)
                                                                                                                               cgccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATCGCGAATGGTTAAGTGGTGACAAACCCACGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACg                     <  2:928367/146‑1 (MQ=38)
                                                                                                                                gccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGCCgg                    >  1:160556/4‑150 (MQ=255)
                                                                                                                                 ctaTCCTCAGCGGCTTTCCATAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCTCTGTATCGCAACTATCCGAGCTGTAAAGCGGTGGTGCATTTGCACAGCACATTGTCGACGgc                   <  2:392885/148‑1 (MQ=255)
                                                                                                                                 cGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgc                   >  2:48864/1‑150 (MQ=255)
                                                                                                                                 cGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgc                   >  2:1529150/1‑150 (MQ=255)
                                                                                                                                 cGAACCGCAGCGGCGTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgc                   <  1:1153779/150‑1 (MQ=255)
                                                                                                                                  caTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgcg                  >  1:812336/2‑150 (MQ=255)
                                                                                                                                   aTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCAACTGGAGCTGTATCGCAACAATCAGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgcgc                 >  2:200848/1‑150 (MQ=255)
                                                                                                                                     ccGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCtt               <  1:968720/150‑1 (MQ=255)
                                                                                                                                       gCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGATTGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATGCGCGCTGTATAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTc             >  1:198938/1‑150 (MQ=255)
                                                                                                                                       gCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTc             >  2:1636209/1‑150 (MQ=255)
                                                                                                                                                  ccAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAg  >  1:1610495/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAG  >  NC_000913/2863631‑2863926

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: