Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,202,860 A→C 7.8% Q228P (CAG→CCG)  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,8600AC 7.8% 269.1 / 9.4 100Q228P (CAG→CCG) xcpYtype II secretion pathway protein XcpY
Reads supporting (aligned to +/- strand):  ref base A (48/44);  new base C (3/5);  total (51/49)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 102 total aligned reads displayed.

TGGCGTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAG  >  NC_002947/1202711‑1203002
                                                                                                                                                     |                                                                                                                                              
tGGCGTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCcagcagca                                                                                                                                                <  2:766720/150‑1 (MQ=255)
   cGTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCTACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCCGCCGCCTCCCGCCAGCAGCAGCg                                                                                                                                             <  2:945825/150‑1 (MQ=255)
    gTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCgg                                                                                                                                            <  1:491178/150‑1 (MQ=255)
       ccGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCAACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctg                                                                                                                                         >  1:497980/1‑150 (MQ=255)
        cGCAGACACCGTAGGCAGTGCTGCACAGGACGTCGCTGCGCGGCAGTTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctgc                                                                                                                                        >  2:884882/1‑150 (MQ=255)
         gCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctgct                                                                                                                                       >  1:10711/1‑150 (MQ=255)
         gCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctgct                                                                                                                                       >  1:177346/1‑150 (MQ=255)
         ccAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctgct                                                                                                                                       <  2:662915/149‑1 (MQ=255)
           aGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGCAAGATTATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTgg                                                                                                                                     <  1:1534218/150‑1 (MQ=255)
           aGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCCGCCCGCAAGGTGCGCTCGTTGCCGGCCCCCTCCCGCCAGCAGCCGCCGCTGCTGGc                                                                                                                                    <  1:1570151/150‑1 (MQ=255)
            gCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCCGCTGCTGCCgga                                                                                                                                    <  1:65954/150‑2 (MQ=255)
            gCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGc                                                                                                                                    >  2:1239555/1‑150 (MQ=255)
            gCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGc                                                                                                                                    <  2:124982/150‑1 (MQ=255)
            gCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCACGCCAGCAGCAGCGGCTGCTGGc                                                                                                                                    <  1:1416677/150‑1 (MQ=255)
              cACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCATCCGCTGCTGCTGGcgg                                                                                                                                  <  2:597332/150‑3 (MQ=255)
               tccGTGGGCAGTGCTGCAAAGGAAGCCGCTAACCGGAAATTGGACAAACGTATTTGCGCCGCAAGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGc                                                                                                                                 <  2:863955/149‑1 (MQ=255)
                 cGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCtt                                                                                                                               <  2:220660/150‑1 (MQ=255)
                      gCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTg                                                                                                                          >  2:806666/1‑150 (MQ=255)
                       cAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGGCCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTgg                                                                                                                         <  2:617644/150‑1 (MQ=255)
                        aGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGt                                                                                                                        >  1:786606/1‑150 (MQ=255)
                        aGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGAGCGACTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGt                                                                                                                        <  1:117045/150‑1 (MQ=255)
                                acaGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTa                                                                                                                >  2:524049/1‑150 (MQ=255)
                                acaGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAACGCGCGCCCGTTGCCGGCCGCCCCCCGCCAGCAACAGCGGCTGCTGGCACCTTGCCTACTGCTGGCTa                                                                                                                >  2:585945/1‑150 (MQ=255)
                                    gACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGt                                                                                                            >  1:1239976/1‑150 (MQ=255)
                                      cGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGttt                                                                                                          <  2:10711/150‑1 (MQ=255)
                                      cGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGttt                                                                                                          <  1:131298/150‑1 (MQ=255)
                                        tCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTgg                                                                                                        <  2:803476/150‑1 (MQ=255)
                                        tCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTgg                                                                                                        <  1:165029/150‑1 (MQ=255)
                                         cGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTggg                                                                                                       <  2:1520770/150‑1 (MQ=255)
                                            tGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTg                                                                                                    >  2:1700916/1‑150 (MQ=255)
                                               gCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCt                                                                                                                           >  2:843237/1‑124 (MQ=255)
                                               gCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCt                                                                                                                           <  1:843237/124‑1 (MQ=255)
                                                 ggCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGCCAGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtg                                                                                               >  2:59719/1‑150 (MQ=255)
                                                  gCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtgt                                                                                              <  2:241189/150‑1 (MQ=255)
                                                   cAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtgtg                                                                                             <  1:270111/150‑1 (MQ=255)
                                                   cAATTGCCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtgtg                                                                                             <  1:693645/150‑1 (MQ=255)
                                                      ttGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                          <  2:1239976/150‑1 (MQ=255)
                                                         gCCAACCGTATTTGCGCCGCTCGACACACNGCCCGGCCTGTTCGNCCGCNCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCCCCCGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTGTGGGGTGGCCTGTGGTTGAg                                                                                       >  1:322047/1‑150 (MQ=255)
                                                          ccAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGt                                                                                      >  2:648949/1‑150 (MQ=255)
                                                          ccAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGt                                                                                      >  1:536403/1‑150 (MQ=255)
                                                            aaCCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCa                                                                                    >  2:11773/1‑150 (MQ=255)
                                                             aCCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAg                                                                                   <  2:1282173/150‑1 (MQ=255)
                                                              ccGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGc                                                                                  >  1:1105671/1‑150 (MQ=255)
                                                                  aTTGGCGCCGCTCGACACACTGCCCGGCCTGTCCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATg                                                                              <  2:322047/150‑1 (MQ=255)
                                                                   tttGCGCCGGTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATgg                                                                             <  2:761763/150‑1 (MQ=255)
                                                                   tttGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCtt                                                                                                                               >  1:1575532/1‑100 (MQ=255)
                                                                   tttGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCtt                                                                                                                               <  2:1575532/100‑1 (MQ=255)
                                                                   tttGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTgg                                                                                                        <  2:735927/123‑1 (MQ=255)
                                                                   tttGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTgg                                                                                                        >  1:735927/1‑123 (MQ=255)
                                                                               gACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCaa                                                                 >  2:1880273/1‑150 (MQ=255)
                                                                                    acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                            >  1:637918/1‑150 (MQ=255)
                                                                                    acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                            >  1:784904/1‑150 (MQ=255)
                                                                                    acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                            >  1:174800/1‑150 (MQ=255)
                                                                                    acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                            >  2:737357/1‑150 (MQ=255)
                                                                                    acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCGTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGGCAGGCGCAACTCTg                                                            >  2:75748/1‑150 (MQ=255)
                                                                                     cTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTgg                                                           >  2:30877/1‑150 (MQ=255)
                                                                                      tGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGgc                                                          >  1:1712909/1‑150 (MQ=255)
                                                                                         ccGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCCGCTGCTGCTGGc                                                                                                                                    >  1:832884/1‑73 (MQ=39)
                                                                                         ccGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCCGCTGCTGCTGGc                                                                                                                                    <  2:832884/73‑1 (MQ=39)
                                                                                          cGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                          <  2:10662/114‑1 (MQ=255)
                                                                                          cGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                          >  1:10662/1‑114 (MQ=255)
                                                                                              cTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCGGCTTGCCTGGTGCTGGCTACCGTACACGTTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAg                                                  >  1:497621/1‑150 (MQ=255)
                                                                                               tGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                            >  1:510311/1‑139 (MQ=255)
                                                                                               tGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                            <  2:510311/139‑1 (MQ=255)
                                                                                               tGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAgg                                                 >  2:1804172/1‑150 (MQ=255)
                                                                                                         cgcCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGt                                       <  1:694119/150‑1 (MQ=255)
                                                                                                            ccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTc                                                                         >  2:763995/1‑113 (MQ=255)
                                                                                                            ccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTc                                                                         <  1:763995/113‑1 (MQ=255)
                                                                                                            ccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGAc                                    >  2:225452/1‑150 (MQ=255)
                                                                                                            ccGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGGTGAGTCAGCAATGGCGTCAGGCGCCACTCTGGCGCAGCCAGGTGATTGCCGTGAc                                    >  1:217547/1‑150 (MQ=255)
                                                                                                              gCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCg                                  >  1:1227780/1‑150 (MQ=255)
                                                                                                                  ggTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGa                              <  1:1804172/150‑1 (MQ=255)
                                                                                                                  ggTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGa                              <  2:177346/150‑1 (MQ=255)
                                                                                                                     gcgcTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCa                           <  2:786606/150‑1 (MQ=255)
                                                                                                                      cgcTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAg                          <  1:806666/150‑1 (MQ=255)
                                                                                                                      cgcTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAg                          <  1:1510545/150‑1 (MQ=255)
                                                                                                                      cgcTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAg                          <  2:1435659/150‑1 (MQ=255)
                                                                                                                      cagtCGGTGCCGTGCACCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAg                          <  2:1047884/147‑1 (MQ=255)
                                                                                                                       gcTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAgg                         >  2:700901/1‑150 (MQ=255)
                                                                                                                        cTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGc                        <  1:700901/150‑1 (MQ=255)
                                                                                                                         tCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                            >  1:167767/1‑113 (MQ=255)
                                                                                                                         tCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                            <  2:167767/113‑1 (MQ=255)
                                                                                                                         tCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAggag                       >  1:1047884/1‑148 (MQ=255)
                                                                                                                             tGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGcc                   <  1:648949/150‑1 (MQ=255)
                                                                                                                                cggccgCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTc                >  1:323623/1‑150 (MQ=255)
                                                                                                                                cggccgCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTc                >  2:984070/1‑150 (MQ=255)
                                                                                                                                 ggccgCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCg               <  1:11773/150‑1 (MQ=255)
                                                                                                                                 ggccgCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGCGACCggg                                >  1:268891/1‑133 (MQ=255)
                                                                                                                                 ggccgCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGCGACCggg                                <  2:268891/133‑1 (MQ=255)
                                                                                                                                  gccgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACGCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCgg              <  1:290856/150‑1 (MQ=255)
                                                                                                                                   ccgccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGc             >  1:1358089/1‑150 (MQ=255)
                                                                                                                                       cTCGCGCCTGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTCGCTACCGTTTGGGGTGGCCTGTGGGTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGCGCAGGCCAGCCCTCGGCAcgc         >  2:27584/1‑150 (MQ=255)
                                                                                                                                        gcccGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGAAATGGCGTCAGGCGCAACTCTGGCTCAGCCAGGTCATTGCCTTGACCGGGGAGCAGGCCAGCCCTCGGCAcgcg        >  1:1716092/2‑150 (MQ=255)
                                                                                                                                          ctgCCAGCAGCCGCTGCTGCTGGCGGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCg               >  2:263817/3‑141 (MQ=255)
                                                                                                                                          ctgCCAGCAGCCGCTGCTGCTGGCGGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCg               <  1:263817/139‑1 (MQ=255)
                                                                                                                                          ccGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAAt                                                                               <  2:1003576/77‑1 (MQ=255)
                                                                                                                                          ccGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAAt                                                                               >  1:1003576/1‑77 (MQ=255)
                                                                                                                                          ccGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGc      >  2:1232450/1‑150 (MQ=255)
                                                                                                                                           tgCCAGCAGCCGCTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGca     >  1:591494/2‑150 (MQ=255)
                                                                                                                                              cagcagcagCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAg  >  2:1588338/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TGGCGTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAG  >  NC_002947/1202711‑1203002

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