Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,353,932 C→T 7.8% H338H (CAC→CAT yrpB → 2‑nitropropane dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,353,9320CT 7.8% 283.2 / 9.5 99H338H (CAC→CATyrpB2‑nitropropane dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (35/56);  new base T (4/4);  total (39/60)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 102 total aligned reads displayed.

ATGCGCGAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAG  >  NC_002947/4353783‑4354078
                                                                                                                                                     |                                                                                                                                                  
aTGCGCGAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGcac                                                                                                                                                    >  2:547472/1‑150 (MQ=255)
  gcgcGAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACAt                                                                                                                                                  <  1:468066/150‑1 (MQ=255)
    gcgAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGc                                                                                                                                                >  2:104355/1‑150 (MQ=255)
        gCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCCGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCcggc                                                                                                                                            >  2:1561290/1‑150 (MQ=255)
           ggggCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCATGTGCcggccta                                                                                                                                         <  1:1735317/150‑4 (MQ=255)
              gCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGa                                                                                                                                      >  1:849098/1‑150 (MQ=255)
              gCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGa                                                                                                                                      <  2:849098/150‑1 (MQ=255)
                cGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGc                                                                                                                                    <  2:1724487/150‑1 (MQ=255)
                cGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGc                                                                                                                                    <  2:1775303/150‑1 (MQ=255)
                 gATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCt                                                                                                                                   <  2:165623/150‑1 (MQ=255)
                        gagcTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCgg                                                                                                                            <  1:653844/150‑1 (MQ=255)
                         agcTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCggg                                                                                                                           <  2:1532297/150‑1 (MQ=255)
                          gcTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGga                                                                                                                          <  2:1602016/150‑1 (MQ=255)
                          gcTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGga                                                                                                                          <  2:114064/150‑1 (MQ=255)
                              gCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGTCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGAtt                                                                                                                      <  1:1751681/150‑1 (MQ=255)
                              gCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGAtt                                                                                                                      >  1:1139733/1‑150 (MQ=255)
                                aCCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGc                                                                                                                    >  1:1228999/1‑150 (MQ=255)
                                aCCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGc                                                                                                                    >  2:24018/1‑150 (MQ=255)
                                 ccGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGcc                                                                                                                   <  1:587455/150‑1 (MQ=255)
                                 ccGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATNGcc                                                                                                                   <  2:1198118/150‑1 (MQ=255)
                                      cTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTaa                                                                                                              <  2:595329/150‑1 (MQ=255)
                                      cTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTaa                                                                                                              <  2:1783735/150‑1 (MQ=255)
                                      cTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTaa                                                                                                              <  2:1784040/150‑1 (MQ=255)
                                       ttCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGAACGCTCCCTTTGCCGCGGGGTATTTATGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCggagat                                                                                                             >  2:65560/1‑146 (MQ=255)
                                        tCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAgg                                                                                                            <  1:1801697/150‑1 (MQ=255)
                                         cccGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGc                                                                                                           >  1:1920310/1‑150 (MQ=255)
                                           cGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCAt                                                                                                         >  1:27991/1‑150 (MQ=255)
                                            gCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCAtt                                                                                                        <  2:1139733/150‑1 (MQ=255)
                                              tGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTgg                                                                                                      >  2:1663555/1‑150 (MQ=255)
                                              tGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGAACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTgg                                                                                                      <  1:1460837/150‑1 (MQ=255)
                                                 ccggctggGCATTGATGCCTTTACGCGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAACGCATTGGCAg                                                                                                   <  1:1617606/148‑1 (MQ=255)
                                                  gggcgggGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGt                                                                                                  <  1:608515/150‑1 (MQ=255)
                                                  gggcgggGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGt                                                                                                  <  1:22349/150‑1 (MQ=255)
                                                     cgggGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGAt                                                                                               <  2:1503199/150‑1 (MQ=255)
                                                       gggCCTTGATGCCATTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCCTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGACGGTAAGGCATTGGCAGTGATCg                                                                                             <  2:946164/150‑1 (MQ=255)
                                                       gggCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCg                                                                                             <  1:854446/150‑1 (MQ=255)
                                                        ggCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGa                                                                                                                                 <  1:45005/113‑1 (MQ=255)
                                                        ggCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGa                                                                                                                                 >  2:45005/1‑113 (MQ=255)
                                                         gCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGGCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGt                                                                                           >  1:1567738/1‑150 (MQ=255)
                                                         gCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTCCAGCAATCTGTGGTCGGGGCAGGCGTTACGTCTAGGCCGGCACATGCCGGCGCGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGt                                                                                           <  2:978627/150‑1 (MQ=255)
                                                               aTGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCgg                                                                                                                            <  1:378502/111‑1 (MQ=255)
                                                               aTGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCgg                                                                                                                            >  2:378502/1‑111 (MQ=255)
                                                               aTGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAg                                                                                     <  1:319019/150‑1 (MQ=255)
                                                               aTGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGCGATCGGTCACCAg                                                                                     >  1:912028/1‑150 (MQ=255)
                                                                   cTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGc                                                                                                     <  1:671968/130‑1 (MQ=255)
                                                                   cTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGc                                                                                                     >  2:671968/1‑130 (MQ=255)
                                                                    tttACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCATGTGCcggccta                                                                                                                                         <  2:945196/93‑4 (MQ=255)
                                                                    tttACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCATGTGCcggccta                                                                                                                                         >  1:945196/1‑90 (MQ=255)
                                                                         gggCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGAc                                                                           >  2:443082/1‑150 (MQ=255)
                                                                          ggCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGAca                                                                          >  1:925319/1‑150 (MQ=255)
                                                                             gATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGAcaggc                                                                       >  1:1677587/1‑150 (MQ=255)
                                                                              aTCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCa                                                                                      <  1:1003559/134‑1 (MQ=255)
                                                                              aTCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCa                                                                                      >  2:1003559/1‑134 (MQ=255)
                                                                               tCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCATGTGCcggcc                                                                                                                                           >  2:1680534/1‑79 (MQ=255)
                                                                               tCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCATGTGCcggcc                                                                                                                                           >  2:1680205/1‑79 (MQ=255)
                                                                               tCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCATGTGCcggcc                                                                                                                                           <  1:1680534/80‑2 (MQ=255)
                                                                               tCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCATGTGCcggcc                                                                                                                                           <  1:1680205/80‑2 (MQ=255)
                                                                                   cGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCTGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCAt                                                                 >  1:1569535/1‑150 (MQ=255)
                                                                                   cGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGAACATGCGCCAGGCATTCTGACAGGCAGGCAt                                                                 <  1:988674/150‑1 (MQ=255)
                                                                                    gagctgCAGTGTAAGGGTGATTGCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATg                                                                <  2:832950/145‑1 (MQ=255)
                                                                                      tCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTGCGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGc                                                              <  2:1892628/150‑1 (MQ=255)
                                                                                      tCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGc                                                              <  2:246690/150‑1 (MQ=255)
                                                                                      tCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGc                                                              <  2:27991/150‑1 (MQ=255)
                                                                                      tCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGc                                                              <  1:104355/150‑1 (MQ=255)
                                                                                      tCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGc                                                              <  1:547472/150‑1 (MQ=255)
                                                                                              gTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGa                                                                            <  2:596646/128‑1 (MQ=255)
                                                                                              gTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGa                                                                            >  1:596646/1‑128 (MQ=255)
                                                                                                aaGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACgg                                                    <  1:343119/150‑1 (MQ=255)
                                                                                                    gTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGga                                                                                                                          <  2:1452984/76‑1 (MQ=255)
                                                                                                    gTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGga                                                                                                                          >  1:1452984/1‑76 (MQ=255)
                                                                                                         ttCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCa                                                                                         <  2:192724/104‑1 (MQ=255)
                                                                                                         ttCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCa                                                                                         >  1:192724/1‑104 (MQ=255)
                                                                                                           cAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTTGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTa                                         <  1:777527/150‑1 (MQ=255)
                                                                                                            aGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTAt                                        <  1:1399084/150‑1 (MQ=255)
                                                                                                            aGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTAt                                        >  2:1439259/1‑150 (MQ=255)
                                                                                                                 tCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCta                                   >  1:821015/1‑150 (MQ=255)
                                                                                                                 tCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCta                                   >  1:999586/1‑150 (MQ=255)
                                                                                                                      ggTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAg                              <  2:130401/150‑1 (MQ=255)
                                                                                                                       gTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGt                             >  1:1775539/1‑150 (MQ=255)
                                                                                                                       gTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGt                             >  2:255833/1‑150 (MQ=255)
                                                                                                                          ggggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTaa                          <  2:508233/150‑1 (MQ=255)
                                                                                                                          ggggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTaa                          <  1:373281/150‑1 (MQ=255)
                                                                                                                           gggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAg                         <  2:925319/150‑1 (MQ=255)
                                                                                                                            ggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                        <  1:443082/150‑1 (MQ=255)
                                                                                                                             gCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                        >  2:1709652/1‑149 (MQ=255)
                                                                                                                             gCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                        <  1:1709652/149‑1 (MQ=255)
                                                                                                                                 gCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAg                                                                                                   >  1:1300174/1‑70 (MQ=255)
                                                                                                                                 gCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAg                                                                                                   <  2:1300174/70‑1 (MQ=255)
                                                                                                                                    ttACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACg                >  1:1088906/1‑150 (MQ=255)
                                                                                                                                        ggCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATaa            <  2:507497/150‑1 (MQ=255)
                                                                                                                                         gCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAAc           <  1:24018/150‑1 (MQ=255)
                                                                                                                                          cTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAAcc          >  2:175752/1‑150 (MQ=255)
                                                                                                                                           tAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                        >  1:1600200/1‑135 (MQ=255)
                                                                                                                                           tAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                        <  2:1600200/135‑1 (MQ=255)
                                                                                                                                            ccgCCGGCATGTGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTc        >  2:316507/3‑150 (MQ=255)
                                                                                                                                               ccGGCACATGCCGGCGGGTGCGCAGACCAGGGAGATCTTCTCTAAGGCGTTGGCATTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAg     <  2:1647314/150‑1 (MQ=255)
                                                                                                                                                  gCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAg  >  1:403778/1‑150 (MQ=255)
                                                                                                                                                  gCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAg  <  1:1062022/150‑1 (MQ=255)
                                                                                                                                                    acaTGCCGGCGGGTGAGCTGACCCGTGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTaa                          <  1:964020/124‑1 (MQ=255)
                                                                                                                                                    acaTGCCGGCGGGTGAGCTGACCCGTGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGACCTCTTCACGGCCAATGGCTTATCGCTATATAGTTaa                          >  2:964020/1‑124 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATGCGCGAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAG  >  NC_002947/4353783‑4354078

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