Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,696 A→C 5.5% intergenic (‑424/+27) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,6960AC 5.5% 358.1 / 10.4 128intergenic (‑424/+27)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base A (64/57);  new base C (4/3);  total (68/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 = 1.14e-01

Only 100 of 129 total aligned reads displayed.

CGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTT  >  NC_000913/4094549‑4094845
                                                                                                                                                   |                                                                                                                                                      
cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAgg                                                                                                                                                      >  1:669000/1‑150 (MQ=34)
cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAgg                                                                                                                                                      <  2:1504971/150‑1 (MQ=34)
cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGAACCGTAGGCCGGATAAgg                                                                                                                                                      >  1:1501172/1‑150 (MQ=21)
      tggcgcgACTAAGCAGGCTTACTACTTTTGACGCCCTGTAGGCCGGTAATGGCACTTCCGCCGCATCCGTCATTCATTGCCTGATGCGACGCTGTCGCGTCTTATTAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTcg                                                                                                                                                <  2:106947/147‑1 (MQ=11)
         cgcgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCgcgc                                                                                                                                             <  1:1403614/150‑1 (MQ=40)
          gcgTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                            <  1:651599/150‑1 (MQ=255)
          gcgTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCACATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                            <  1:782416/150‑1 (MQ=255)
          acgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGACGTAACAGGCCTACCACGGATGACGACCCGTAGGCCGGATAAGGCGCTCgcgca                                                                                                                                            >  1:683097/2‑149 (MQ=14)
           ggTCTTATCAGGCCTACAACTATTTACGAACTGTAGGCCTGATAAGGCACTTGTGCCGCAGCCGGCATTCAATGCTTGATTCGACGGTGTCGCGTCTTATCAGGCCTACAACTTTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgccg                                                                                                                                           <  2:1648210/149‑1 (MQ=11)
           cttCTTATCAGCCCTACAACTATTGCCGACCTGTAGTCCGGATTAGGCACTTGAGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgccg                                                                                                                                           <  1:465142/148‑1 (MQ=255)
           cgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgccg                                                                                                                                           <  1:148988/150‑1 (MQ=255)
              cTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCAt                                                                                                                                        >  1:1026358/1‑150 (MQ=255)
                tATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCAGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATcc                                                                                                                                      >  1:439246/1‑150 (MQ=255)
                tATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTAGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCGTACGACTGTTGCCGACCCGTAGGCCGCATAAGGCGCTCGCGCCGCAGcc                                                                                                                                      >  2:592899/1‑150 (MQ=14)
                 aTCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCg                                                                                                                                     >  2:884230/1‑150 (MQ=255)
                   cAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                   <  1:757003/150‑1 (MQ=11)
                    aGGCCTTCAACTTTTGCCGACCTGTAGGCCGGATAAGTCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                  <  2:172678/150‑2 (MQ=11)
                     ggCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGcct                                                                                                                                 <  2:693642/150‑3 (MQ=11)
                     ggCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAg                                                                                                                                 <  1:1452647/150‑1 (MQ=255)
                     ggCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCAACAACTGTTGCCGACCCGTAGGCCGGATAATGCGCTCGCGCCGCATCCGGCAg                                                                                                                                 >  2:1556534/1‑150 (MQ=255)
                       ccTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAgtg                                                                                                                               >  1:1689232/1‑150 (MQ=255)
                        cTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAgtgt                                                                                                                              >  2:1155560/1‑150 (MQ=255)
                         tACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGtt                                                                                                                             <  2:173675/150‑1 (MQ=255)
                                aTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACTCTGTAGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCGCGCGCCGCATCCGGCAGTGTTTACCcgc                                                                                                                      >  2:274025/1‑150 (MQ=255)
                                   gCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCgcggc                                                                                                                   >  2:1156266/1‑150 (MQ=255)
                                    ccGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCgcggcg                                                                                                                  >  1:1467098/1‑150 (MQ=255)
                                        ccTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTc                                                                                                              >  2:531520/1‑150 (MQ=255)
                                          tGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCaa                                                                                                            >  1:1373396/1‑150 (MQ=255)
                                           gattgCCGGATGCGGCACAAGTGCCTTATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCaaa                                                                                                           >  2:1596749/5‑150 (MQ=18)
                                              ggCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtt                                                                                                        <  1:686867/150‑1 (MQ=255)
                                               gCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGCCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAttt                                                                                                       >  2:1398465/1‑150 (MQ=255)
                                                ccGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttc                                                                                                      <  1:82018/150‑1 (MQ=255)
                                                ccGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttc                                                                                                      >  2:422793/1‑150 (MQ=255)
                                                 cGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttct                                                                                                    >  2:718472/1‑150 (MQ=255)
                                                  ggATAAGGCACTTGTGCCGCATCCGTCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttctt                                                                                                   <  2:1087078/150‑1 (MQ=255)
                                                   gATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttcttt                                                                                                  <  1:884230/150‑1 (MQ=255)
                                                   gATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCACATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttcttt                                                                                                  >  2:548808/1‑150 (MQ=255)
                                                    aTACGGCACTGGTGCCGCATCCGGCAATCAATGCCTGCTGCGACGCTGTCGCGTCTTAGCAGGCCTACCACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttctttc                                                                                                 <  1:1103586/150‑1 (MQ=18)
                                                     tAAGGCCCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATttctttct                                                                                                <  2:226092/150‑1 (MQ=255)
                                                     tAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATttctttct                                                                                                <  2:1413145/150‑1 (MQ=255)
                                                      aaGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTtctc                                                                                               <  1:1310879/150‑1 (MQ=255)
                                                       aGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCa                                                                                              >  2:627190/1‑150 (MQ=255)
                                                       aGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCa                                                                                              <  1:370341/150‑1 (MQ=255)
                                                        ggCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCAt                                                                                             >  1:1235031/1‑150 (MQ=255)
                                                        ggCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCAt                                                                                             >  2:566601/1‑150 (MQ=255)
                                                        ggCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGGTTACCCGCGGCGACTCAAAATTTCGTTCTCaa                                                                                             >  1:1088113/1‑149 (MQ=255)
                                                         gCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATa                                                                                            >  1:238376/1‑150 (MQ=255)
                                                         gCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATa                                                                                            >  2:352629/1‑150 (MQ=255)
                                                          cACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGCCTCAAAATTTCTTTCTCATaa                                                                                           >  1:629644/1‑150 (MQ=255)
                                                          cACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATaa                                                                                           >  2:1057559/1‑150 (MQ=255)
                                                           aCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAg                                                                                          <  1:1167936/150‑1 (MQ=255)
                                                                cgtcGCATCAGGCATTGATTGCCGGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGc                                                                                     >  2:186725/4‑150 (MQ=21)
                                                                 gCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCaaa                                                                                                           <  1:1240232/128‑1 (MQ=255)
                                                                   cGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACg                                                                                  <  2:915884/150‑1 (MQ=255)
                                                                   cGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACg                                                                                  <  1:318997/150‑1 (MQ=255)
                                                                   cGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGGCTTTTCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACg                                                                                  <  1:719958/150‑1 (MQ=255)
                                                                    gCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGc                                                                                 >  2:1387029/1‑150 (MQ=255)
                                                                        ccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCtctc                                                                             >  1:740575/1‑150 (MQ=255)
                                                                         cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTcc                                                                            >  2:1004645/1‑150 (MQ=255)
                                                                         cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCACAAGCCCGCACGCTCTcc                                                                            <  2:1543628/150‑1 (MQ=255)
                                                                              atcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCCCGCGCTCCAGccc                                                                       >  1:739021/1‑149 (MQ=255)
                                                                                caatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCagt                                                                     >  1:1232498/1‑148 (MQ=255)
                                                                                caatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCAACCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCCtt                                                                     >  2:727738/1‑150 (MQ=255)
                                                                                 aatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTc                                                                    >  2:1184582/1‑150 (MQ=255)
                                                                                  atGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCg                                                                   <  1:1183028/150‑1 (MQ=255)
                                                                                     ccTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCAACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTa                                                                >  2:1053258/1‑150 (MQ=255)
                                                                                        gATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACct                                                             >  2:1703152/1‑150 (MQ=255)
                                                                                         aTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACctg                                                            <  2:223747/150‑1 (MQ=255)
                                                                                              aCGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTgg                                                        <  2:1689232/150‑1 (MQ=255)
                                                                                                 cTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGCCCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTGTACCCGCGGCGACGCAAAATTTCGTTCTCATAAGCCCGCACGCTCTCCAGCCCTTCGCGACCTGCTGGCGTa                                                    >  2:1720737/1‑150 (MQ=37)
                                                                                                      gcgTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTg                                               <  2:1653767/150‑1 (MQ=255)
                                                                                                       cgTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGa                                              >  2:1404884/1‑150 (MQ=255)
                                                                                                             aTCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTgg                                        <  2:1026358/150‑1 (MQ=255)
                                                                                                              tCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGc                                       <  1:1057559/150‑1 (MQ=255)
                                                                                                               cAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCa                                      >  2:1036100/1‑150 (MQ=255)
                                                                                                                aGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCaa                                     >  2:477146/1‑150 (MQ=255)
                                                                                                                 ggCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAAt                                    <  1:929212/150‑1 (MQ=255)
                                                                                                                  gCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATa                                   <  1:662749/150‑1 (MQ=255)
                                                                                                                  gCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCAGCACGCTATACAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATa                                   >  2:1592547/1‑150 (MQ=255)
                                                                                                                  gCCTACAACTGTTGCCGAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATa                                   <  1:1662512/150‑1 (MQ=255)
                                                                                                                   ccTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATAc                                  <  1:50362/150‑1 (MQ=255)
                                                                                                                   ccTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCACATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATAc                                  <  2:669000/150‑1 (MQ=255)
                                                                                                                    cTACAACTGTTGCCGACCCGTAGGCCGGATTAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGCCGTTGGCAATACa                                 >  2:1620478/1‑150 (MQ=255)
                                                                                                                      aCAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACAtt                               <  1:191026/150‑1 (MQ=255)
                                                                                                                                ccGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCg                     >  1:868689/1‑150 (MQ=255)
                                                                                                                                 cGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGc                    >  2:930467/1‑150 (MQ=255)
                                                                                                                                   aCCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGccg                  >  1:1429115/1‑149 (MQ=255)
                                                                                                                                    cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTg                 <  2:740575/150‑1 (MQ=255)
                                                                                                                                    cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTg                 <  2:749382/150‑1 (MQ=255)
                                                                                                                                    cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTg                 <  1:718472/150‑1 (MQ=255)
                                                                                                                                        actgCCGGATGCGGCGCGAGCGCCTTATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGcccc             >  2:1584280/4‑148 (MQ=255)
                                                                                                                                        actgCCGGATGCGGCGCGAGCGCCTTATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCAc             >  1:1213730/4‑150 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACgg           <  1:778054/150‑1 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGCTCGCGCAGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCGGGCGTATCGTGCCGTTGGCAAGACATTTCCCAGCCCGCCTGCCACgg           >  1:334212/1‑150 (MQ=255)
                                                                                                                                           gCCGGATGCGGCGCGAGCGCCTTATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGc          >  1:535870/1‑150 (MQ=255)
                                                                                                                                           gCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGc          <  2:621459/150‑1 (MQ=255)
                                                                                                                                           gCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGc          <  1:684744/150‑1 (MQ=255)
                                                                                                                                            ccGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGct         <  1:352629/150‑2 (MQ=255)
                                                                                                                                              ggATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAAc       <  1:627190/150‑1 (MQ=255)
                                                                                                                                               gATAAGGCGCTCGCGCAGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACg      <  2:683097/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTT  >  NC_000913/4094549‑4094845

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: