Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,202,863 G→T 16.9% R229L (CGG→CTG)  xcpY → type II secretion pathway protein XcpY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,202,8630GT16.9% 159.7 / 21.4 65R229L (CGG→CTG) xcpYtype II secretion pathway protein XcpY
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base T (5/6);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

CGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCTCAAG  >  NC_002947/1202719‑1203011
                                                                                                                                                |                                                                                                                                                    
cGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctgc                                                                                                                                                 >  2:1083147/1‑150 (MQ=255)
        ccGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACCCCCCCCTCCCGCTGTTCGTCCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCTGCCAGCAGCCGCTGCTGCTGGcggga                                                                                                                                         <  1:146370/150‑5 (MQ=255)
             ggCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCt                                                                                                                                    <  1:658397/150‑1 (MQ=255)
                                   cTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCCGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGt                                                                                                              <  1:580609/150‑1 (MQ=255)
                                          gCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGc                                                                                                                                             >  1:521940/1‑112 (MQ=255)
                                          gCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGc                                                                                                                                             <  2:521940/112‑1 (MQ=255)
                                          gCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtgt                                                                                                       >  2:303513/1‑150 (MQ=255)
                                             aTTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGt                                                                                                    <  1:1122845/150‑1 (MQ=255)
                                              ttGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                                   <  2:462418/150‑1 (MQ=255)
                                              ttGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                                   <  2:69839/150‑1 (MQ=255)
                                                ggCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTgg                                                                                                                            <  1:1177681/123‑1 (MQ=255)
                                                ggCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTgg                                                                                                                            >  2:1177681/1‑123 (MQ=255)
                                                 gCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTgg                                                                                                     >  1:679351/1‑145 (MQ=255)
                                                 gCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTgg                                                                                                     <  2:679351/145‑1 (MQ=255)
                                                 gCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAg                                                                                                >  1:1099475/1‑150 (MQ=255)
                                                  ccAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGt                                                                                               <  1:464827/150‑1 (MQ=255)
                                                       cGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGa                                                                                                 <  2:18826/143‑1 (MQ=255)
                                                       cGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGa                                                                                                 >  1:18826/1‑143 (MQ=255)
                                                            ttGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGc                                                                                     >  2:1055753/1‑150 (MQ=255)
                                                                gccgcTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTgg                                                                                                                            <  2:711519/107‑1 (MQ=255)
                                                                gccgcTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTgg                                                                                                                            >  1:711519/1‑107 (MQ=255)
                                                                gccgcTCGACACACTGCACGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCa                                                                                 >  2:107435/1‑150 (MQ=255)
                                                                 ccgcTCGACACACTGCCCGGCCTGGGGGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAg                                                                                >  1:766624/1‑150 (MQ=255)
                                                                    cTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGgcg                                                                             >  2:923285/1‑150 (MQ=255)
                                                                    cTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGACGCCTCCCGCCAACAGAAACGGCTGCTGACAGATTGCCTCGTGCGGGCTACCGGGTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGgcg                                                                             >  2:1040965/1‑150 (MQ=255)
                                                                            acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                                     >  2:808760/1‑150 (MQ=255)
                                                                              tGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCTGCCAGCAGCCGCTGCTGCTGGcggga                                                                                                                                         <  2:1155599/80‑5 (MQ=37)
                                                                              tGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCTGCCAGCAGCCGCTGCTGCTGGcggga                                                                                                                                         >  1:1155599/1‑76 (MQ=37)
                                                                              tGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGgc                                                                   <  1:1039145/150‑1 (MQ=255)
                                                                              tGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGgc                                                                   <  1:107435/150‑1 (MQ=255)
                                                                                    gCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGcc                                                             >  2:1060226/1‑150 (MQ=255)
                                                                                       tGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAgg                                                          >  2:312798/1‑150 (MQ=255)
                                                                                         ttCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGTCTTTGGCTACCGTTTGGGGTGAGACGCAGTTTAGTTGGCAAGGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTc                                                        >  1:270599/1‑150 (MQ=255)
                                                                                            gACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGa                                                                                                 <  1:596216/106‑1 (MQ=255)
                                                                                            gACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGa                                                                                                 >  2:596216/1‑106 (MQ=255)
                                                                                               cgcgCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGcc                                                  <  1:312798/150‑1 (MQ=255)
                                                                                                   cccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGC‑‑CGCTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGttt                                                                                                                   <  1:641225/78‑1 (MQ=255)
                                                                                                   cccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGC‑‑CGCTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGttt                                                                                                                   >  2:641225/1‑78 (MQ=255)
                                                                                                    ccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGAc                                             >  1:228757/1‑150 (MQ=255)
                                                                                                      gCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCg                                           <  1:1055753/150‑1 (MQ=255)
                                                                                                       cAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCTTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCgg                                          <  2:766624/150‑1 (MQ=255)
                                                                                                              cgcTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAg                                   <  1:303513/150‑1 (MQ=255)
                                                                                                                cTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGc                                 <  1:524463/150‑1 (MQ=255)
                                                                                                                cTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGc                                 <  1:1040965/150‑1 (MQ=255)
                                                                                                                 tCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGAcc                                            >  2:983972/1‑138 (MQ=255)
                                                                                                                 tCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGAcc                                            <  1:983972/138‑1 (MQ=255)
                                                                                                                 tCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGcc                                >  1:798557/1‑150 (MQ=255)
                                                                                                                        cggccgCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTc                         >  2:764181/1‑150 (MQ=255)
                                                                                                                            cgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCa                     <  1:300955/150‑1 (MQ=255)
                                                                                                                            cgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGGGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCa                     <  2:285780/150‑1 (MQ=255)
                                                                                                                                  ccGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAg                                                           <  2:1136118/106‑1 (MQ=255)
                                                                                                                                  ccGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGc               >  1:1142503/1‑150 (MQ=255)
                                                                                                                                  ccGCCAGCAGCAGCGGCTGCTGGCAGCTTGACTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAg                                                           >  1:1136118/1‑106 (MQ=255)
                                                                                                                                   tgCCAGCAGCCGCTGCTGCTGGCGGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCa                                                                                          <  2:359439/73‑1 (MQ=37)
                                                                                                                                   tgCCAGCAGCCGCTGCTGCTGGCGGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCa                                                                                          >  1:359439/2‑74 (MQ=37)
                                                                                                                                     ccAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGcaca            >  1:111555/1‑150 (MQ=255)
                                                                                                                                      cagcagccgCTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGc                                                              <  2:690134/99‑1 (MQ=255)
                                                                                                                                      cagcagccgCTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGc                                                              >  1:690134/1‑99 (MQ=255)
                                                                                                                                      cagcagcagCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGcc                                                             >  1:885948/1‑100 (MQ=255)
                                                                                                                                      cagcagcagCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGcc                                                             <  2:885948/100‑1 (MQ=255)
                                                                                                                                      cagcagcagCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAg           >  1:159477/1‑150 (MQ=255)
                                                                                                                                       agcagccgcTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGcc                            <  1:549624/132‑1 (MQ=255)
                                                                                                                                       agcagccgcTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGcc                            >  2:549624/1‑132 (MQ=255)
                                                                                                                                          agcagcGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGccc       >  2:345278/1‑150 (MQ=255)
                                                                                                                                               cGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCTCAAg  >  2:1054795/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAGGCCCTCAAG  >  NC_002947/1202719‑1203011

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: