Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,602,097 T→A 5.1% Q479L (CAG→CTG)  PP_4924 ← subtilase family serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,602,0970TA 5.1% 378.0 / 5.7 117Q479L (CAG→CTG) PP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base T (51/59);  new base A (2/4);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

TGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTC  >  NC_002947/5601951‑5602245
                                                                                                                                                  |                                                                                                                                                    
tGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACCGTGCAAGCGCTTCGGAAGATCTCCTCGGAGCTCAGACGCACCGCTCCGACCGGCTTGCCCTggg                                                                                                                                                   >  2:1232455/1‑149 (MQ=255)
   cACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGa                                                                                                                                                >  1:2042546/1‑150 (MQ=255)
    aCCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAg                                                                                                                                               >  2:52833/1‑150 (MQ=255)
    aCCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAg                                                                                                                                               <  1:32220/150‑1 (MQ=255)
    aCCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCAAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGGCCAGGTCGAg                                                                                                                                               >  1:1936654/1‑150 (MQ=255)
      cctcctGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGCCCTGGcc                                                                                                                                                  <  2:1612221/145‑1 (MQ=255)
      cctcctGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCt                                                                                                                                             >  2:2229148/1‑150 (MQ=255)
        tcctGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTg                                                                                                                                           >  1:355579/1‑150 (MQ=255)
        ccctGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCCCGGCCAGGTCGAGccgg                                                                                                                                           <  2:2279208/149‑4 (MQ=255)
           tGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGGCCAGGTCGAGccggtag                                                                                                                                        <  1:2356323/150‑7 (MQ=255)
           tGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATc                                                                                                                                        <  1:107857/150‑1 (MQ=255)
           tGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCAAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGGCCAGGTCGAGccggtag                                                                                                                                        <  2:1936654/150‑7 (MQ=255)
              gCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTgg                                                                                                                                     <  1:164643/150‑1 (MQ=255)
                   gATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAAt                                                                                                                                <  2:1572008/150‑1 (MQ=255)
                    aTCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATg                                                                                                                               <  1:1335547/150‑1 (MQ=255)
                       cGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCa                                                                                                                            <  1:2229148/150‑1 (MQ=255)
                         aTTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGTGCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATg                                                                                                                          >  2:467875/1‑150 (MQ=255)
                              cGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCt                                                                                                                     <  2:939430/150‑1 (MQ=255)
                              cGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCt                                                                                                                     <  1:1419629/150‑1 (MQ=255)
                                   gTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAAt                                                                                                                <  2:449985/150‑1 (MQ=255)
                                      tCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGTCCGATCTTTATCTTGCTAATGCCATGCTCCTGCAATGcc                                                                                                             >  2:1111393/1‑150 (MQ=255)
                                        ggagagaAACCTCCCCGCCACTACTTCCCGGGCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTTGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCt                                                                                                           <  1:37946/147‑1 (MQ=255)
                                         gCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCtt                                                                                                          <  1:1781062/150‑1 (MQ=255)
                                          cGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGCGCTTGCTCTGGCTAATGCCCTGCTCCTGCAATGCCCTTa                                                                                                         <  2:1241133/150‑1 (MQ=255)
                                          cGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTa                                                                                                         <  2:135741/150‑1 (MQ=255)
                                          cGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTa                                                                                                         <  1:300503/150‑1 (MQ=255)
                                             gaAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCATCTTGCAGGTCTGCTACCGGCTCGACCTTGCCGAGCTTGCTCTGCTTTTTGCTCGGCCCCGGCAATGCCCTTAGCt                                                                                                      >  1:1252137/1‑150 (MQ=255)
                                                aCCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCAc                                                                                                   <  1:100840/150‑1 (MQ=255)
                                                 ccTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGGCCAGGTCGAGccg                                                                                                                                            >  1:537062/1‑106 (MQ=255)
                                                 ccTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGGCCAGGTCGAGccg                                                                                                                                            <  2:537062/108‑3 (MQ=255)
                                                 ccTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCAGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGa                                                                                                                                          >  2:1032651/1‑110 (MQ=255)
                                                 ccTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCAGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGa                                                                                                                                          <  1:1032651/110‑1 (MQ=255)
                                                  cTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAGAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCGGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCg                                                                                                 >  1:2302926/1‑150 (MQ=255)
                                                  cTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTg                                                                                                                                      >  2:1421146/1‑113 (MQ=255)
                                                   tCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCAt                                                                                                                           <  2:1843917/123‑1 (MQ=255)
                                                   tCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCAt                                                                                                                           >  1:1843917/1‑123 (MQ=255)
                                                     aCCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCa                                                                                              >  2:1581008/1‑150 (MQ=255)
                                                        gCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCg                                                                                           <  1:467875/150‑1 (MQ=255)
                                                         ccAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGAATAGGGAAGTGAACTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCgg                                                                                          >  1:2102767/1‑150 (MQ=255)
                                                           aGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCa                                                                                        <  1:1602594/150‑1 (MQ=255)
                                                            gTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCAc                                                                                       >  2:1993172/1‑150 (MQ=255)
                                                             tACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACg                                                                                      >  2:278312/1‑150 (MQ=255)
                                                                ttCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGccc                                                                                   <  1:278312/150‑1 (MQ=255)
                                                                ttCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGccc                                                                                   >  2:851678/1‑150 (MQ=255)
                                                                 tCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTgg                                                                                                                                     <  1:1230641/99‑1 (MQ=255)
                                                                 tCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTgg                                                                                                                                     >  2:1230641/1‑99 (MQ=255)
                                                                   ccGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCcgtc                                                                                <  2:2102767/150‑1 (MQ=255)
                                                                       tCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTgg                                                                            <  1:1993172/150‑1 (MQ=255)
                                                                          tCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTg                                                                             <  1:1126159/146‑1 (MQ=255)
                                                                          tCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTg                                                                             >  2:1126159/1‑146 (MQ=255)
                                                                            aCCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGc                                                                                                                   >  1:712795/1‑106 (MQ=255)
                                                                            aCCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGc                                                                                                                   <  2:712795/106‑1 (MQ=255)
                                                                            aCCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGcc                                                                       >  2:1387088/1‑150 (MQ=255)
                                                                             ccATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCa                                                                      <  2:1252137/150‑1 (MQ=255)
                                                                              cATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCAt                                                                     <  1:52833/150‑1 (MQ=255)
                                                                               aTTAACAGTGAAAGCGCCTCGGGAAATGACCGCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATc                                                                    >  2:1485473/1‑150 (MQ=255)
                                                                                  aaCAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGAc                                                                 <  1:1581008/150‑1 (MQ=255)
                                                                                     aGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACTCCCGTCGTGGTTGCCATCGAc                                                                 >  2:265921/1‑147 (MQ=255)
                                                                                     aGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGAc                                                                 <  1:265921/147‑1 (MQ=255)
                                                                                      gTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCcac                                                             <  2:2042546/150‑1 (MQ=255)
                                                                                       tGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATg                                                                                                                          >  1:255618/1‑88 (MQ=255)
                                                                                        gAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCacac                                                           <  1:1758394/150‑1 (MQ=255)
                                                                                        gAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCAAGCCCGTCGTGGTTGCCATCGACCCacac                                                           >  1:1732095/1‑150 (MQ=255)
                                                                                           aGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAGTCTCCGGTTTCTACAATGCCCTTAGCTCAACGGCATCGCCACGCCCGTCTTAGTTGCCATCGACCCACACCtt                                                        >  1:418864/1‑150 (MQ=255)
                                                                                              ctcTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCaa                                                     <  1:871760/148‑1 (MQ=255)
                                                                                                cTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTGCCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAAtt                                                   >  2:1393720/1‑150 (MQ=255)
                                                                                                    ggAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTgcgc                                               <  2:2302926/150‑1 (MQ=255)
                                                                                                       aaTGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCa                                            <  1:1485473/150‑1 (MQ=255)
                                                                                                        aTGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAg                                           >  1:1534089/1‑150 (MQ=255)
                                                                                                             cTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACTGGCATCGGCACGCCCGTCGTGGTTGCCATCGAccc                                                               >  1:1423827/1‑125 (MQ=255)
                                                                                                             cTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACTGGCATCGGCACGCCCGTCGTGGTTGCCATCGAccc                                                               <  2:1423827/125‑1 (MQ=255)
                                                                                                             cTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGc                                      <  1:1111393/150‑1 (MQ=255)
                                                                                                                gTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTc                                   <  2:1629792/150‑1 (MQ=255)
                                                                                                                     ccTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGCTCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCaa                              >  2:649915/1‑150 (MQ=255)
                                                                                                                     ccTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGAt                                                                                                                                         <  2:466071/43‑1 (MQ=255)
                                                                                                                     ccTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTAAGCCCGGACGGCGTCTTCaa                              >  2:1340439/1‑150 (MQ=255)
                                                                                                                         acggtccTCTCCTCCCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAACTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACCCCTTCAATTGCGCCCAGACGGCGTCTTCAACGtt                          <  2:284647/145‑1 (MQ=255)
                                                                                                                          cGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATACGATCCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTtg                         <  1:1353545/150‑1 (MQ=255)
                                                                                                                           gCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTtga                        >  1:2212020/1‑150 (MQ=255)
                                                                                                                            cACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTtgat                       <  2:355579/150‑1 (MQ=255)
                                                                                                                             aCCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTtgatg                      >  2:1974625/1‑150 (MQ=255)
                                                                                                                              ccTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTtgatga                     >  2:1252127/1‑150 (MQ=255)
                                                                                                                              ccTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTtgatga                     >  2:1252124/1‑150 (MQ=255)
                                                                                                                               ctctCCTACCGGCTCGACCTGGCCGTGCTTGATCTGGCTAATGCCCTGCAGATACAATGCCTTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCGTTCGGCGTCTTCAACGTTGATGAc                    <  2:65783/150‑1 (MQ=255)
                                                                                                                               ctctCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGAc                    >  2:2296389/1‑150 (MQ=255)
                                                                                                                               ctctCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGAc                    >  2:2257085/1‑150 (MQ=255)
                                                                                                                               ctctCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGAc                    >  1:1305030/1‑150 (MQ=255)
                                                                                                                                 ctcCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGAcc                                                                <  1:596605/104‑1 (MQ=255)
                                                                                                                                 ctcCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGAcc                                                                >  2:596605/1‑104 (MQ=255)
                                                                                                                                 ctcCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGc                  >  1:1648736/1‑150 (MQ=255)
                                                                                                                                  tcCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGc                                                                                                                              >  2:1870135/1‑41 (MQ=255)
                                                                                                                                    cTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAt               >  2:292714/1‑150 (MQ=255)
                                                                                                                                     tACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAAc                             >  2:312482/1‑135 (MQ=255)
                                                                                                                                     tACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAAc                             <  1:312482/135‑1 (MQ=255)
                                                                                                                                     tACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATc              <  1:2257085/150‑1 (MQ=255)
                                                                                                                                      aCCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCGATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCg             >  1:562645/1‑150 (MQ=255)
                                                                                                                                        cGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCg             >  2:1599745/1‑148 (MQ=255)
                                                                                                                                        cGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCg             <  1:1599745/148‑1 (MQ=255)
                                                                                                                                        cGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCATCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGc           >  1:592997/1‑150 (MQ=255)
                                                                                                                                                 cTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTc  <  1:1340439/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTC  >  NC_002947/5601951‑5602245

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