Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,886,825 C→G 7.4% A83P (GCA→CCA)  PP_5161 ← lysophospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,886,8250CG 7.4% 298.6 / 12.9 95A83P (GCA→CCA) PP_5161lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (44/44);  new base G (5/2);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

GCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGC  >  NC_002947/5886676‑5886971
                                                                                                                                                     |                                                                                                                                                  
gcTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                                    <  1:1608494/150‑1 (MQ=255)
      cAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAGGCCTGccc                                                                                                                                              >  2:1041444/1‑150 (MQ=255)
         gCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGag                                                                                                                                           <  2:751805/150‑2 (MQ=255)
         gCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGa                                                                                                                                            >  1:1871909/1‑149 (MQ=255)
         gCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGa                                                                                                                                            <  2:1871909/149‑1 (MQ=255)
             tGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAgg                                                                                                                                       <  1:492337/150‑1 (MQ=255)
                   ggCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAAc                                                                                                                                 >  2:1772457/1‑150 (MQ=255)
                         tCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCt                                                                                                                           <  2:1697699/150‑1 (MQ=255)
                         tCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCt                                                                                                                           >  1:535547/1‑150 (MQ=255)
                          cacaGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTg                                                                                                                          >  2:1633493/1‑150 (MQ=255)
                          cacaGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCACGTCGAACCCCGCTg                                                                                                                          >  2:779368/1‑150 (MQ=255)
                           acaGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGa                                                                                                                                            >  2:60184/1‑131 (MQ=255)
                           acaGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGa                                                                                                                                            <  1:60184/131‑1 (MQ=255)
                             aGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCt                                                                                                                       <  2:852083/150‑1 (MQ=255)
                               cTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGa                                                                                                                     >  2:526305/1‑150 (MQ=255)
                                tGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGaa                                                                                                                    >  1:1890629/1‑150 (MQ=255)
                                 gATCACGGCGAAACCCTGCTTGAGTGCCAACTCGACCGCAAAGGGCAAGGGGTTGATGTGGTCGTAATAGCCGGGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGaaa                                                                                                                   <  1:1867538/150‑1 (MQ=255)
                                   tCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCTTGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGc                                                                                                                 >  1:1079769/1‑150 (MQ=255)
                                      cGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTAc                                                                                                              <  2:1769388/150‑1 (MQ=255)
                                      cGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTAc                                                                                                              <  1:526305/150‑1 (MQ=255)
                                         cGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTa                                                                                                               >  2:1440395/1‑146 (MQ=255)
                                         cGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTa                                                                                                               <  1:1440395/146‑1 (MQ=255)
                                          gAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACccag                                                                                                          >  2:1171079/1‑150 (MQ=255)
                                             aCCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACccagcca                                                                                                       >  2:210093/1‑150 (MQ=255)
                                               ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCcagccagc                                                                                                     >  1:1736639/1‑150 (MQ=255)
                                               ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCcagccagc                                                                                                     >  1:21275/1‑150 (MQ=255)
                                               ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCcagccagc                                                                                                     >  1:1090649/1‑150 (MQ=255)
                                                  gCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCGGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                  <  2:1225186/150‑1 (MQ=255)
                                                  gCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                  <  1:1633493/150‑1 (MQ=255)
                                                  gCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                  <  2:1337049/150‑1 (MQ=255)
                                                    ttGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                <  1:701124/150‑1 (MQ=255)
                                                    ttGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                <  1:844505/150‑1 (MQ=255)
                                                           cccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACccatcc                                                                                                        >  1:691714/1‑132 (MQ=255)
                                                           cccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACccatcc                                                                                                        <  2:691714/135‑1 (MQ=255)
                                                            ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAg                                                                                        <  2:1148841/150‑1 (MQ=255)
                                                            ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAg                                                                                        >  2:997580/1‑150 (MQ=255)
                                                            ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAg                                                                                        <  2:535547/150‑1 (MQ=255)
                                                             cACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAgg                                                                                       >  1:1560583/1‑150 (MQ=255)
                                                              aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGGGCTTGACTGCTCGGGCAGCCATGCCTGTCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGt                                                                                      >  1:147068/1‑150 (MQ=255)
                                                                tCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCTGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTc                                                                                    >  1:220824/1‑150 (MQ=255)
                                                                tCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCTATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTc                                                                                    <  2:1090649/150‑1 (MQ=255)
                                                                tCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTc                                                                                    >  2:1518609/1‑150 (MQ=255)
                                                                tCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTc                                                                                    >  1:1309139/1‑150 (MQ=255)
                                                                  gATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCt                                                                                                                       <  2:182656/113‑1 (MQ=255)
                                                                  gATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCt                                                                                                                       >  1:182656/1‑113 (MQ=255)
                                                                    tCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGc                                                                                <  2:1104915/150‑1 (MQ=255)
                                                                    tCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGc                                                                                <  2:193553/150‑1 (MQ=255)
                                                                         tGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAg                                                                           <  1:246522/150‑1 (MQ=255)
                                                                           tCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGaa                                                                         <  1:1216577/150‑1 (MQ=255)
                                                                           tCGGTAGAGGCCAATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGaa                                                                         <  2:1031572/150‑1 (MQ=255)
                                                                            cGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAg                                                                        >  1:274189/1‑150 (MQ=255)
                                                                               tAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAgcgc                                                                     <  1:788761/150‑1 (MQ=255)
                                                                                agagGCCCAGGTGGTCGTAATAGCCGTGCAGCAGGACCAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGGTGCCTGAACGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCt                                                                    <  2:279359/150‑1 (MQ=255)
                                                                                 gagGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCTGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTg                                                                   >  1:9813/1‑150 (MQ=255)
                                                                                  gggCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                  <  2:1237094/149‑1 (MQ=255)
                                                                                  agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                  <  1:1238257/150‑1 (MQ=255)
                                                                                  agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                  <  2:1309139/150‑1 (MQ=255)
                                                                                  agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                  <  1:1447291/150‑1 (MQ=255)
                                                                                  agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                  <  2:21275/150‑1 (MQ=255)
                                                                                  agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                  <  2:1592889/150‑1 (MQ=255)
                                                                                      ccATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGt                                                              >  1:848533/1‑150 (MQ=255)
                                                                                        aTGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                                    >  2:293217/1‑62 (MQ=255)
                                                                                        aTGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                                    <  1:293217/62‑1 (MQ=255)
                                                                                           tgGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAGCCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTgc                                                         <  1:1847057/150‑1 (MQ=255)
                                                                                           tgGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTgc                                                         <  1:1216375/150‑1 (MQ=255)
                                                                                             gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtg                                                                                                   <  1:93430/106‑1 (MQ=255)
                                                                                             gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtg                                                                                                   >  2:93430/1‑106 (MQ=255)
                                                                                             gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCAt                                                                             >  2:324645/1‑128 (MQ=255)
                                                                                             gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCAt                                                                             <  1:324645/128‑1 (MQ=255)
                                                                                                  aaTAGCCGTGCAGCAGGAAAAGCGTCGCGCTTAGCTGCTCTGGCAGCCATGCCTGCCCGACCAGGTCGCGACCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgtg                                                  <  1:1772457/150‑1 (MQ=255)
                                                                                                   aTAGCCGTGCAGCAGGAACACCGTCGCCCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATATAAGCGCTGGTAGTCCTGCGCTTGTGc                                                 <  1:227195/150‑1 (MQ=255)
                                                                                                         gTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAg                                           <  2:1495507/150‑1 (MQ=255)
                                                                                                         gTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAg                                           >  1:236319/1‑150 (MQ=255)
                                                                                                             agcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggc                                       >  2:1257872/1‑150 (MQ=255)
                                                                                                                agGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctgg                                    <  2:437591/150‑1 (MQ=255)
                                                                                                                  gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCg                                  >  2:360694/1‑150 (MQ=255)
                                                                                                                  gAACAGCGTCGCGCTTGGCTGCACGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCg                                  >  2:551624/1‑150 (MQ=255)
                                                                                                                   aaCAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGa                                 <  2:274189/150‑1 (MQ=255)
                                                                                                                           cgcgCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGaa                         <  2:1890629/150‑1 (MQ=255)
                                                                                                                           caggcatGGCTGCCCGAGCAGCCAAGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGaa                         >  1:241125/7‑150 (MQ=255)
                                                                                                                               cTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgc                     >  2:1221202/1‑150 (MQ=255)
                                                                                                                               cTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgc                     >  1:1133524/1‑150 (MQ=255)
                                                                                                                                ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgcg                    >  1:459127/1‑150 (MQ=255)
                                                                                                                                ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgcg                    >  1:1253600/1‑150 (MQ=255)
                                                                                                                                  ggCTGCCCGAGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGcc                  >  2:1346336/1‑150 (MQ=255)
                                                                                                                                   gctgctCGGGCAGGCAGGCCTGCCCGACCAGGTCGAACCCCGCTGCCGGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCa                 <  2:1643685/150‑1 (MQ=255)
                                                                                                                                    ctggtCGGGCAGGCATGGCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg                >  1:265160/1‑150 (MQ=255)
                                                                                                                                      ggtCGGGCAGGCATGGCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGc              >  2:1862069/3‑150 (MQ=255)
                                                                                                                                       gtCGGGCAAGCATGGCTGCCCGAGCAGGTCGAACCCCGCTGCCGGAAAGCTACCCAGTCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCg             >  2:1500453/2‑150 (MQ=255)
                                                                                                                                        tCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATGGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg            >  2:447994/1‑150 (MQ=255)
                                                                                                                                          gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGt          <  1:360694/150‑1 (MQ=255)
                                                                                                                                           ggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTg         >  1:1217436/1‑150 (MQ=255)
                                                                                                                                            gCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgc        >  1:1588317/1‑150 (MQ=255)
                                                                                                                                             cAGGCATGGCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgcg       >  1:305059/1‑150 (MQ=255)
                                                                                                                                                ccATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGACAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg    >  2:1430836/1‑150 (MQ=255)
                                                                                                                                                 cATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTCGTGCTCACCTGAAGCCCTCACGGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAgg   >  1:1376546/1‑150 (MQ=255)
                                                                                                                                                 cATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAgg   <  1:997580/150‑1 (MQ=255)
                                                                                                                                                  aTGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGc  >  1:1082965/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGC  >  NC_002947/5886676‑5886971

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