Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 152,158 A→T 9.4% intergenic (‑52/+42) PP_0143 ← / ← PP_0144 membrane protein/insulinase family metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947152,1580AT 9.4% 259.6 / 18.1 98intergenic (‑52/+42)PP_0143/PP_0144membrane protein/insulinase family metalloprotease
Reads supporting (aligned to +/- strand):  ref base A (46/42);  new base T (4/5);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

CGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTA  >  NC_002947/152013‑152304
                                                                                                                                                 |                                                                                                                                                  
cGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGccc                                                                                                                                                >  1:1557502/1‑150 (MQ=255)
 gATACCGGAGAGTAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTTAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcg                                                                                                                                               >  2:1240733/1‑150 (MQ=255)
 gATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGAAAGGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAACTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcg                                                                                                                                               >  1:687124/1‑150 (MQ=255)
  aTACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcgc                                                                                                                                              <  2:2105044/150‑1 (MQ=255)
  aTACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcgc                                                                                                                                              <  2:17965/150‑1 (MQ=255)
   tACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcgcg                                                                                                                                             >  1:276054/1‑150 (MQ=255)
     ccGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                           <  1:1084500/150‑1 (MQ=255)
      cGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaag                                                                                                                                          <  2:1557502/150‑2 (MQ=255)
      cGGAGAGGAAACCGCCGAGCGCAATGATGCCAGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGaaa                                                                                                                                          >  2:744491/1‑150 (MQ=255)
        gagaGTAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTTAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAgg                                                                                                                                        <  1:1240733/150‑1 (MQ=255)
        gagaGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaagag                                                                                                                                        <  1:62873/150‑4 (MQ=255)
        gagaGGAAACCGCCGAGCGCAATGATGCCCCGGTAGCGCTTGAGTAGTTCGACGATGCGCTAGCCGAAATGCGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTACCCTCCTCTTCGCGGGCAAGCCCGCGAAAgg                                                                                                                                        <  1:859000/150‑1 (MQ=255)
               aaCCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTa                                                                                                                                 >  2:1504738/1‑150 (MQ=255)
                   gccgAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTgg                                                                                                                             <  1:1132032/150‑1 (MQ=255)
                        gcgcAATGGTGGCCGGGTAGCGCTTGAGTAGTTCGATGAGGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCaaaa                                                                                                                        <  2:2471152/150‑1 (MQ=255)
                        gcgcAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCaaaa                                                                                                                        <  2:468404/150‑1 (MQ=255)
                          gcAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGAGAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTGGGTAATGGCAAAAgg                                                                                                                      >  1:1448453/1‑150 (MQ=255)
                           cAATGATGCCAGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGa                                                                                                                     <  1:744491/150‑1 (MQ=255)
                           cAATGAGGCCCGGGTAGCGGTTCAGTAGTTCGATGATGCGCTGGACGATATGGGGGTTTTGGGGCAGTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCGCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGa                                                                                                                     <  2:715674/150‑1 (MQ=255)
                              tgatgCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAc                                                                                                                  <  2:2896275/150‑1 (MQ=255)
                                   cccGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCaaaa                                                                                                                        <  2:1061460/139‑1 (MQ=255)
                                   cccGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCaaaa                                                                                                                        >  1:1061460/1‑139 (MQ=255)
                                    ccGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGa                                                                                                            >  1:380686/1‑150 (MQ=255)
                                        gTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTa                                                                                                        <  2:687124/150‑1 (MQ=255)
                                         tAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTTTGTCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTAc                                                                                                       <  1:333858/150‑1 (MQ=255)
                                         tAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTAc                                                                                                       <  1:354427/150‑1 (MQ=255)
                                          aGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTAcc                                                                                                      <  1:1799029/150‑1 (MQ=255)
                                             gcTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCt                                                                                                   >  1:2787360/1‑150 (MQ=255)
                                             gcTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCt                                                                                                   >  2:1678894/1‑150 (MQ=255)
                                                 gAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAg                                                                                                                       <  1:2194033/126‑1 (MQ=255)
                                                 gAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAg                                                                                                                       >  2:2194033/1‑126 (MQ=255)
                                                   gtagtTCGAGGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtg                                                                                             >  1:2431842/1‑150 (MQ=255)
                                                    tagtTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtgt                                                                                            <  2:2424836/150‑1 (MQ=255)
                                                     agtTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTa                                                                                           <  2:380686/150‑1 (MQ=255)
                                                      gtTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaaga                                                                                                                                         >  2:2323848/1‑101 (MQ=255)
                                                      gtTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGaaga                                                                                                                                         <  1:2323848/103‑3 (MQ=255)
                                                       ttCGATGCAGCGCTGCACCGTCTGCGGGATCCGGGGCCTTGAGGGCTACGAGAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGTGCAAGCCGGCGAAAGGGGCGTTACTGGAAAAAAGATACCCAAGATTTACCGCTGGCGAGTAGc                                                                                         >  2:66013/1‑150 (MQ=255)
                                                       ttCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGc                                                                                         >  2:828415/1‑150 (MQ=255)
                                                         cGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcac                                                                                       <  2:2816974/150‑1 (MQ=255)
                                                          gatgatGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTgg                                                                                                                             >  2:270773/1‑111 (MQ=255)
                                                          gatgatGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTgg                                                                                                                             <  1:270773/111‑1 (MQ=255)
                                                           atgatgCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcacca                                                                                     >  1:1518219/1‑150 (MQ=255)
                                                             gatgCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACt                                                                                   >  1:494172/1‑150 (MQ=255)
                                                              atgCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTg                                                                                  <  2:2787360/150‑1 (MQ=255)
                                                              atgCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTg                                                                                  <  1:828415/150‑1 (MQ=255)
                                                                 cgcTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGt                                                                                                                                      >  2:2979184/1‑95 (MQ=255)
                                                                 cgcTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGt                                                                                                                                      <  1:2979184/95‑1 (MQ=255)
                                                                  gcTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGcc                                                                              >  1:1681437/1‑150 (MQ=255)
                                                                    tGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTATGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGt                                                                            >  2:2532869/1‑150 (MQ=255)
                                                                      gACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtgt                                                                                            <  2:924529/132‑1 (MQ=255)
                                                                      gACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtgt                                                                                            >  1:924529/1‑132 (MQ=255)
                                                                      gACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTc                                                                          >  1:134708/1‑150 (MQ=255)
                                                                              gggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGc                                                                               >  2:2405526/1‑137 (MQ=255)
                                                                              gggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCAt                                                                  <  1:261717/150‑1 (MQ=255)
                                                                              gggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCAt                                                                  <  2:361307/150‑1 (MQ=255)
                                                                              gggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGCATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGc                                                                               <  1:2405526/137‑1 (MQ=255)
                                                                              gggggATATGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGCGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATTCCGCCGGTCTACCGCTGCCGTGTAGTACCACTGGGCCGTTCAGGTCCAt                                                                  <  2:836459/150‑1 (MQ=255)
                                                                              ggggaATTTGGGGCATTTGTGGCTACGACAGAAAAGTGAGTTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCTGTACTGGCAAAAGTATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCAt                                                                  <  1:718853/150‑1 (MQ=255)
                                                                               ggggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATg                                                                 >  2:180205/1‑150 (MQ=255)
                                                                                gggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACt                                                                                   >  1:1251237/1‑131 (MQ=255)
                                                                                gggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACt                                                                                   <  2:1251237/131‑1 (MQ=255)
                                                                                gggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGtt                                                                           >  2:1866595/1‑139 (MQ=255)
                                                                                gggATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGtt                                                                           <  1:1866595/139‑1 (MQ=255)
                                                                                       ggggCATAGGGGGGTACGAGAGAAAAGTGAGGTGGCTGGGCCAGCGTCTTCGTGGGCGAGGCCGGGGGAGGGTGGGTACGGGCAAAAGGATAACGGCGACTTACCGCTGCCGTGTGGGGACACTTGGCCGTTTATGGCCGTgcggcggcc                                                         >  2:584369/1‑150 (MQ=255)
                                                                                          gCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAgc                                                      <  2:1448453/150‑1 (MQ=255)
                                                                                          gCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATACACGGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAgc                                                      >  1:967571/1‑150 (MQ=255)
                                                                                              tGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAg                                                                        <  1:36565/128‑1 (MQ=255)
                                                                                              tGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAg                                                                        >  2:36565/1‑128 (MQ=255)
                                                                                                 gggCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAg                                               >  2:2885878/1‑150 (MQ=255)
                                                                                                     tACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGcccc                                           <  2:134708/150‑1 (MQ=255)
                                                                                                        gACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCccgcc                                        >  2:736182/1‑150 (MQ=255)
                                                                                                          aaGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAGGGGTCCGTACTGGCAAAAAGATACCGTCGATTTACCGCTGCCGTCGAGGACCACTGGGCCGTTCACGTCCAGGCAGCCGCCAGCGCAACAGCCCcgccgc                                      <  2:276054/149‑1 (MQ=255)
                                                                                                           aGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCgccccc                                     >  2:610681/1‑150 (MQ=255)
                                                                                                             aaaggggAGGTGGCTGTGCCAGCGTCTTAGCGGGCAGGCCCGCAAAATGGGCGGTACTGGCGAAAGGATACCGTCGATTTACAGCGGCCGTGTCGGACCCCTGGGACGTTCAGGTCCATCCCGCCGCAAGCGCAACAGCCCCGCCGCCAg                                   <  2:494172/144‑1 (MQ=255)
                                                                                                                 gTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAcc                                                                                                                 >  2:1402389/1‑68 (MQ=255)
                                                                                                                 gTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAcc                                                                                                                 <  1:1402389/68‑1 (MQ=255)
                                                                                                                 gTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGTCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGccc                               >  1:103201/1‑150 (MQ=255)
                                                                                                                 gTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGccc                               <  1:1504738/150‑1 (MQ=255)
                                                                                                                 gTGAGGTGGCGGTGCCAGCCTCTTCGCGGTCAAGCCCGCGAAAGGGTCGGTACTGTCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCAAGCGCAACAGCCGCGCCGCCAGGccc                               <  1:66013/150‑1 (MQ=255)
                                                                                                                      gTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGccc                          >  1:2315144/1‑150 (MQ=255)
                                                                                                                      gTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGccc                          >  1:2766369/1‑150 (MQ=255)
                                                                                                                        ggCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGccccc                        >  2:1257209/1‑149 (MQ=255)
                                                                                                                           tgtgCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGa                     >  2:399040/1‑150 (MQ=255)
                                                                                                                               ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCaa                 >  1:2813868/1‑150 (MQ=255)
                                                                                                                                cAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGATCAAg                >  2:2844645/1‑150 (MQ=255)
                                                                                                                                cAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAg                >  1:2020544/1‑150 (MQ=255)
                                                                                                                                  gCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGc              >  1:364360/1‑150 (MQ=255)
                                                                                                                                  accctttCGCGGGCTTGCCCGCGAAGGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGc              >  1:175861/6‑150 (MQ=255)
                                                                                                                                    ctctTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCa            <  1:1678894/150‑1 (MQ=255)
                                                                                                                                     tctTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCAt           <  1:2530713/150‑1 (MQ=255)
                                                                                                                                     ctttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTa                                                                                           >  1:2757967/3‑70 (MQ=38)
                                                                                                                                     ctttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTa                                                                                           <  2:2757967/68‑1 (MQ=38)
                                                                                                                                          gcgGGCAAGGCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATgccg                                                              >  1:2881440/1‑94 (MQ=255)
                                                                                                                                          gcgGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATgccg                                                              <  2:2881440/94‑1 (MQ=255)
                                                                                                                                             ggCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGt   >  1:1988213/1‑150 (MQ=255)
                                                                                                                                              gCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTa  >  2:2153674/1‑150 (MQ=255)
                                                                                                                                              gCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTa  <  1:736182/150‑1 (MQ=255)
                                                                                                                                              gCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATAACGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTa  >  1:2069211/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTA  >  NC_002947/152013‑152304

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