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.

CGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTT  >  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)
                                                 cGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttc‑t                                                                                                    >  2:718472/1‑150 (MQ=255)
                                                  ggATAAGGCACTTGTGCCGCATCCGTCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttc‑tt                                                                                                   <  2:1087078/150‑1 (MQ=255)
                                                   gATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttc‑ttt                                                                                                  <  1:884230/150‑1 (MQ=255)
                                                   gATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCACATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttc‑ttt                                                                                                  >  2:548808/1‑150 (MQ=255)
                                                    aTACGGCACTGGTGCCGCATCCGGCAATCAATGCCTGCTGCGACGCTGTCGCGTCTTAGCAGGCCTACCACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAAtttc‑tttc                                                                                                 <  1:1103586/150‑1 (MQ=18)
                                                     tAAGGCCCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATttc‑tttct                                                                                                <  2:226092/150‑1 (MQ=255)
                                                     tAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATttc‑tttct                                                                                                <  2:1413145/150‑1 (MQ=255)
                                                      aaGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTtctc                                                                                               <  1:1310879/150‑1 (MQ=255)
                                                       aGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCa                                                                                              >  2:627190/1‑150 (MQ=255)
                                                       aGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCa                                                                                              <  1:370341/150‑1 (MQ=255)
                                                        ggCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCAt                                                                                             >  1:1235031/1‑150 (MQ=255)
                                                        ggCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCAt                                                                                             >  2:566601/1‑150 (MQ=255)
                                                        ggCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGGTTACCCGCGGCGACTCAAAATTTC‑GTTCTCaa                                                                                             >  1:1088113/1‑149 (MQ=255)
                                                         gCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATa                                                                                            >  1:238376/1‑150 (MQ=255)
                                                         gCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATa                                                                                            >  2:352629/1‑150 (MQ=255)
                                                          cACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGCCTCAAAATTTC‑TTTCTCATaa                                                                                           >  1:629644/1‑150 (MQ=255)
                                                          cACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATaa                                                                                           >  2:1057559/1‑150 (MQ=255)
                                                           aCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAg                                                                                          <  1:1167936/150‑1 (MQ=255)
                                                                cgtcGCATCAGGCATTGATTGCCGGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGc                                                                                     >  2:186725/4‑150 (MQ=21)
                                                                 gCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCaaa                                                                                                           <  1:1240232/128‑1 (MQ=255)
                                                                   cGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACg                                                                                  <  2:915884/150‑1 (MQ=255)
                                                                   cGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACg                                                                                  <  1:318997/150‑1 (MQ=255)
                                                                   cGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGGCTTTTCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACg                                                                                  <  1:719958/150‑1 (MQ=255)
                                                                    gCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGc                                                                                 >  2:1387029/1‑150 (MQ=255)
                                                                        ccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCtctc                                                                             >  1:740575/1‑150 (MQ=255)
                                                                         cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTcc                                                                            >  2:1004645/1‑150 (MQ=255)
                                                                         cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCACAAGCCCGCACGCTCTcc                                                                            <  2:1543628/150‑1 (MQ=255)
                                                                              atcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCCCGCGCTCCAGccc                                                                       >  1:739021/1‑149 (MQ=255)
                                                                                caatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCagt                                                                     >  1:1232498/1‑148 (MQ=255)
                                                                                caatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCAACCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCCtt                                                                     >  2:727738/1‑150 (MQ=255)
                                                                                 aatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTc                                                                    >  2:1184582/1‑150 (MQ=255)
                                                                                  atGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCg                                                                   <  1:1183028/150‑1 (MQ=255)
                                                                                     ccTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCAACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTa                                                                >  2:1053258/1‑150 (MQ=255)
                                                                                        gATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACct                                                             >  2:1703152/1‑150 (MQ=255)
                                                                                         aTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACctg                                                            <  2:223747/150‑1 (MQ=255)
                                                                                              aCGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTCTTTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTgg                                                        <  2:1689232/150‑1 (MQ=255)
                                                                                                 cTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGCCCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTGTACCCGCGGCGACGCAAAATTTC‑GTTCTCATAAGCCCGCACGCTCTCCAGCCCTTCGCGACCTGCTGGCGTa                                                    >  2:1720737/1‑150 (MQ=37)
                                                                                                      gcgTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTg                                               <  2:1653767/150‑1 (MQ=255)
                                                                                                       cgTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGa                                              >  2:1404884/1‑150 (MQ=255)
                                                                                                             aTCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTgg                                        <  2:1026358/150‑1 (MQ=255)
                                                                                                              tCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGc                                       <  1:1057559/150‑1 (MQ=255)
                                                                                                               cAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCa                                      >  2:1036100/1‑150 (MQ=255)
                                                                                                                aGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCaa                                     >  2:477146/1‑150 (MQ=255)
                                                                                                                 ggCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAAt                                    <  1:929212/150‑1 (MQ=255)
                                                                                                                  gCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATa                                   <  1:662749/150‑1 (MQ=255)
                                                                                                                  gCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCAGCACGCTATACAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATa                                   >  2:1592547/1‑150 (MQ=255)
                                                                                                                  gCCTACAACTGTTGCCGAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATa                                   <  1:1662512/150‑1 (MQ=255)
                                                                                                                   ccTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATAc                                  <  1:50362/150‑1 (MQ=255)
                                                                                                                   ccTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCACATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATAc                                  <  2:669000/150‑1 (MQ=255)
                                                                                                                    cTACAACTGTTGCCGACCCGTAGGCCGGATTAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGCCGTTGGCAATACa                                 >  2:1620478/1‑150 (MQ=255)
                                                                                                                      aCAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACAtt                               <  1:191026/150‑1 (MQ=255)
                                                                                                                                ccGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCg                     >  1:868689/1‑150 (MQ=255)
                                                                                                                                 cGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGc                    >  2:930467/1‑150 (MQ=255)
                                                                                                                                   aCCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGccg                  >  1:1429115/1‑149 (MQ=255)
                                                                                                                                    cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTg                 <  2:740575/150‑1 (MQ=255)
                                                                                                                                    cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTg                 <  2:749382/150‑1 (MQ=255)
                                                                                                                                    cccGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTg                 <  1:718472/150‑1 (MQ=255)
                                                                                                                                        actgCCGGATGCGGCGCGAGCGCCTTATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGcccc             >  2:1584280/4‑148 (MQ=255)
                                                                                                                                        actgCCGGATGCGGCGCGAGCGCCTTATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCAc             >  1:1213730/4‑150 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACgg           <  1:778054/150‑1 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGCTCGCGCAGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCGGGCGTATCGTGCCGTTGGCAAGACATTTCCCAGCCCGCCTGCCACgg           >  1:334212/1‑150 (MQ=255)
                                                                                                                                           gCCGGATGCGGCGCGAGCGCCTTATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGc          >  1:535870/1‑150 (MQ=255)
                                                                                                                                           gCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGc          <  2:621459/150‑1 (MQ=255)
                                                                                                                                           gCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGc          <  1:684744/150‑1 (MQ=255)
                                                                                                                                            ccGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGct         <  1:352629/150‑2 (MQ=255)
                                                                                                                                              ggATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAAc       <  1:627190/150‑1 (MQ=255)
                                                                                                                                               gATAAGGCGCTCGCGCAGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACg      <  2:683097/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGTTTACCCGCGGCGACTCAAAATTTC‑TTTCTCATAAGCCCGCACGCTCTCCAGCCATTCGCTACCTGCTGGCGTATCGTGACGTTGGCAATACATTTCCCAGACCGCCTGCCACGGCAACGATTT  >  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: ‑