Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,602,094 A→G 5.2% V480A (GTC→GCC)  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,0940AG 5.2% 342.1 / 4.4 117V480A (GTC→GCC) PP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base A (51/57);  new base G (2/4);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 119 total aligned reads displayed.

GTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGG  >  NC_002947/5601945‑5602241
                                                                                                                                                     |                                                                                                                                                   
gTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACTGCGCAAGTGCCTCGGGATATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGa                                                                                                                                                     <  1:1500703/150‑1 (MQ=255)
gTAGAGTGCCACACTCCTAGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGa                                                                                                                                                     >  2:1419629/1‑150 (MQ=255)
 tAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGAc                                                                                                                                                    <  1:1564984/150‑1 (MQ=255)
      tGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACCGTGCAAGCGCTTCGGAAGATCTCCTCGGAGCTCAGACGCACCGCTCCGACCGGCTTGCCCTggg                                                                                                                                               >  2:1232455/1‑149 (MQ=255)
          aCCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAg                                                                                                                                           <  1:32220/150‑1 (MQ=255)
          aCCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAg                                                                                                                                           >  2:52833/1‑150 (MQ=255)
          aCCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCAAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGGCCAGGTCGAg                                                                                                                                           >  1:1936654/1‑150 (MQ=255)
            cctcctGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGCCCTGGcc                                                                                                                                              <  2:1612221/145‑1 (MQ=255)
            cctcctGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGcc                                                                                                                                              >  1:1612221/1‑145 (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)
                   ggCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTg                                                                                                                                  <  2:1779920/150‑1 (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)
                                  agcGATGGTATCGGCGAGAAACCTCACTGACAGTACTTCCCAGTCCTCAGCATTAACTGTGAAAGCGCCTCGGGAAATGAGCTAGGTAGCCTGACGAACCTCTCCTACCGGCTCGCCGTGGCCGGGCTTGATCTGCCTAATGCCATGTTc                                                                                                                   <  2:364797/148‑1 (MQ=255)
                                    cGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCt                                                                                                                 <  1:1419629/150‑1 (MQ=255)
                                    cGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCt                                                                                                                 <  2:939430/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)
                                                   gaAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCATCTTGCAGGTCTGCTACCGGCTCGACCTTGCCGAGCTTGCTCTGCTTTTTGCTCGGCCCCGGCAATGCCCTTAGCt                                                                                                  >  1:1252137/1‑150 (MQ=255)
                                                      aCCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCAc                                                                                               <  1:100840/150‑1 (MQ=255)
                                                       ccTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGGCCAGGTCGAGccg                                                                                                                                        <  2:537062/108‑3 (MQ=255)
                                                       ccTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGGCCAGGTCGAGccg                                                                                                                                        >  1:537062/1‑106 (MQ=255)
                                                       ccTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCAGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGa                                                                                                                                      >  2:1032651/1‑110 (MQ=255)
                                                       ccTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCAGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGa                                                                                                                                      <  1:1032651/110‑1 (MQ=255)
                                                         tCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCAt                                                                                                                       >  1:1843917/1‑123 (MQ=255)
                                                         tCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCAt                                                                                                                       <  2:1843917/123‑1 (MQ=255)
                                                           aCCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCa                                                                                          >  2:1581008/1‑150 (MQ=255)
                                                               ccAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGAATAGGGAAGTGAACTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCgg                                                                                      >  1:2102767/1‑150 (MQ=255)
                                                                 aGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCa                                                                                    <  1:1602594/150‑1 (MQ=255)
                                                                  gTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCAc                                                                                   >  2:1993172/1‑150 (MQ=255)
                                                                      ttCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGccc                                                                               >  2:851678/1‑150 (MQ=255)
                                                                      ttCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGccc                                                                               <  1:278312/150‑1 (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                                                                                                               <  2:712795/106‑1 (MQ=255)
                                                                                  aCCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGc                                                                                                               >  1:712795/1‑106 (MQ=255)
                                                                                  aCCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGcc                                                                   >  2:1387088/1‑150 (MQ=255)
                                                                                    cATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCAt                                                                 <  1:52833/150‑1 (MQ=255)
                                                                                     aTTAACAGTGAAAGCGCCTCGGGAAATGACCGCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATc                                                                >  2:1485473/1‑150 (MQ=255)
                                                                                       tAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGa                                                              <  2:132686/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)
                                                                                             tGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATg                                                                                                                      <  2:255618/88‑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)
                                                                                                                           ccTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGAt                                                                                                                                     >  1:466071/1‑43 (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)
                                                                                                                                   aCCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTGAACGTtgatg                  >  1:225318/1‑150 (MQ=255)
                                                                                                                                   aCCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTtgatg                  >  2:1974625/1‑150 (MQ=255)
                                                                                                                                    ccTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTtgatga                 >  2:1252124/1‑150 (MQ=255)
                                                                                                                                    ccTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTtgatga                 >  2:1252127/1‑150 (MQ=255)
                                                                                                                                     ctctCCTACCGGCTCGACCTGGCCGTGCTTGATCTGGCTAATGCCCTGCAGATACAATGCCTTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCGTTCGGCGTCTTCAACGTTGATGAc                <  2:65783/150‑1 (MQ=255)
                                                                                                                                     ctctCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGAc                >  2:2257085/1‑150 (MQ=255)
                                                                                                                                     ctctCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGAc                >  1:1305030/1‑150 (MQ=255)
                                                                                                                                     ctctCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGAc                >  2:2296389/1‑150 (MQ=255)
                                                                                                                                       ctcCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGAcc                                                            >  2:596605/1‑104 (MQ=255)
                                                                                                                                       ctcCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGAcc                                                            <  1:596605/104‑1 (MQ=255)
                                                                                                                                       ctcCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGc              >  1:1648736/1‑150 (MQ=255)
                                                                                                                                        tcCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGc                                                                                                                          >  2:1870135/1‑41 (MQ=255)
                                                                                                                                        tcCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGc                                                                                                                          <  1:1870135/41‑1 (MQ=255)
                                                                                                                                          cTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAt           >  2:292714/1‑150 (MQ=255)
                                                                                                                                          cTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAt           >  2:1464426/1‑150 (MQ=255)
                                                                                                                                          cTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCAACAGCATCGGCACGCCCGTCGTGGTTGCCATCTACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAt           >  2:2079057/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         <  1:1599745/148‑1 (MQ=255)
                                                                                                                                              cGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCATCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGc       >  1:592997/1‑150 (MQ=255)
                                                                                                                                                   cGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTgg  <  1:292714/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGCCAGTACTTCCCGGTCCTCACCATTAACAGTGAAAGCGCCTCGGGAAATGACCTCGTTACCCTGACGCACCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGG  >  NC_002947/5601945‑5602241

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