Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,353,933 A→G 8.0% M339V (ATG→GTG)  yrpB → 2‑nitropropane dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,353,9330AG 8.0% 284.9 / 9.8 99M339V (ATG→GTG) yrpB2‑nitropropane dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (35/56);  new base G (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 = 8.48e-01

Only 100 of 102 total aligned reads displayed.

GCGCGAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAGCTGT  >  NC_002947/4353785‑4354082
                                                                                                                                                    |                                                                                                                                                     
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:1784040/150‑1 (MQ=255)
                                    cTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTaa                                                                                                                  <  2:1783735/150‑1 (MQ=255)
                                    cTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTaa                                                                                                                  <  2:595329/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                                                                                                         >  2:671968/1‑130 (MQ=255)
                                                                 cTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGc                                                                                                         <  1:671968/130‑1 (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                                                                                             >  1:192724/1‑104 (MQ=255)
                                                                                                       ttCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCa                                                                                             <  2:192724/104‑1 (MQ=255)
                                                                                                         cAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTTGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTa                                             <  1:777527/150‑1 (MQ=255)
                                                                                                          aGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTAt                                            >  2:1439259/1‑150 (MQ=255)
                                                                                                          aGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTAt                                            <  1:1399084/150‑1 (MQ=255)
                                                                                                               tCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCta                                       >  1:999586/1‑150 (MQ=255)
                                                                                                               tCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCta                                       >  1:821015/1‑150 (MQ=255)
                                                                                                                    ggTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAg                                  <  2:130401/150‑1 (MQ=255)
                                                                                                                     gTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGt                                 >  2:255833/1‑150 (MQ=255)
                                                                                                                     gTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGt                                 >  1:1775539/1‑150 (MQ=255)
                                                                                                                        ggggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTaa                              <  1:373281/150‑1 (MQ=255)
                                                                                                                        ggggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTaa                              <  2:508233/150‑1 (MQ=255)
                                                                                                                         gggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAg                             <  2:925319/150‑1 (MQ=255)
                                                                                                                          ggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                            <  1:443082/150‑1 (MQ=255)
                                                                                                                           gCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                            <  1:1709652/149‑1 (MQ=255)
                                                                                                                           gCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                            >  2:1709652/1‑149 (MQ=255)
                                                                                                                               gCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAg                                                                                                       <  2:1300174/70‑1 (MQ=255)
                                                                                                                               gCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAg                                                                                                       >  1:1300174/1‑70 (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                            <  2:1600200/135‑1 (MQ=255)
                                                                                                                                         tAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                            >  1:1600200/1‑135 (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)
                                                                                                                                                    aTGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAGCTGt  <  2:1228999/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCGAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAGCTGT  >  NC_002947/4353785‑4354082

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: