Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,886,825 C→G 8.3% 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 8.3% 247.8 / 12.5 84A83P (GCA→CCA) PP_5161lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (41/36);  new base G (4/3);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

GACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGCGGT  >  NC_002947/5886680‑5886974
                                                                                                                                                 |                                                                                                                                                     
gACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgcctgc                                                                                                                                                   >  2:1058780/1‑150 (MQ=255)
    ggCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGa                                                                                                                                               <  2:181182/150‑1 (MQ=255)
     gCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGAc                                                                                                                                              <  2:1458820/150‑1 (MQ=255)
             cGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGa                                                                                                                                      <  2:642938/150‑1 (MQ=255)
              tggCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTAGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGaa                                                                                                                                     >  2:1398952/2‑150 (MQ=255)
                    gTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGa                                                                                                                                               <  1:1122658/134‑1 (MQ=255)
                    gTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGa                                                                                                                                               >  2:1122658/1‑134 (MQ=255)
                     tCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                                       >  1:507468/1‑125 (MQ=255)
                     tCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                                       <  2:507468/125‑1 (MQ=255)
                       acaGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGc                                                                                                                            <  1:154941/150‑1 (MQ=255)
                       acaGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGc                                                                                                                            <  1:787437/150‑1 (MQ=255)
                       acaGCTGATCACGGCGAAACCCTGCTTAAATGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGc                                                                                                                            >  1:1322905/1‑150 (MQ=255)
                         aGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCt                                                                                                                          >  2:216213/1‑150 (MQ=255)
                         aGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCt                                                                                                                          >  2:1093589/1‑150 (MQ=255)
                             gATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAACGTCGCGATTAGTTGCTCGCGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGaaa                                                                                                                      <  1:1494080/150‑1 (MQ=255)
                               tCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGaa                                                                                                                                     >  2:860296/1‑133 (MQ=255)
                               tCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGaa                                                                                                                                     <  1:860296/133‑1 (MQ=255)
                                 aCGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTa                                                                                                                  <  2:857618/150‑1 (MQ=255)
                                  cGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGACAGCTAc                                                                                                                 >  2:1085060/1‑150 (MQ=255)
                                  cGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTAc                                                                                                                 <  1:670435/150‑1 (MQ=255)
                                           ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCcagccagc                                                                                                        >  2:445172/1‑150 (MQ=255)
                                           ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCcagccagc                                                                                                        >  2:865160/1‑150 (MQ=255)
                                            cTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCACGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCt                                                                                                       <  1:1683183/150‑1 (MQ=255)
                                                ttGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATCGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                   >  1:850104/1‑150 (MQ=255)
                                                  gAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGc                                                                                                 <  1:865160/150‑1 (MQ=255)
                                                   aGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGcc                                                                                                <  1:445172/150‑1 (MQ=255)
                                                      gCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGc                                                                                             >  1:1180083/1‑150 (MQ=255)
                                                        ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                     >  1:281754/1‑140 (MQ=255)
                                                        ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                     <  2:281754/140‑1 (MQ=255)
                                                        ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAg                                                                                           >  2:169070/1‑150 (MQ=255)
                                                          aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGccc                                                                                                                                                 >  2:1072330/1‑94 (MQ=255)
                                                          aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGccc                                                                                                                                                 <  1:1072330/94‑1 (MQ=255)
                                                          aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCCCGAGCAGCCAagc                                                                                                                                                       <  1:152933/88‑1 (MQ=255)
                                                              gATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAg                                                                                     <  1:57375/150‑1 (MQ=255)
                                                              gATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAg                                                                                     >  2:57375/1‑150 (MQ=255)
                                                                tCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGc                                                                                   >  2:148789/1‑150 (MQ=255)
                                                                    aTGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCTCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATa                                                                               <  2:1322905/150‑1 (MQ=255)
                                                                     tGTCGGTAGGGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAg                                                                              <  2:231548/150‑1 (MQ=255)
                                                                       tCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGaa                                                                            <  1:783950/150‑1 (MQ=255)
                                                                         ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCGTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAgc                                                                          >  1:1390372/1‑150 (MQ=255)
                                                                         ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAgc                                                                          <  1:216213/150‑1 (MQ=255)
                                                                          gTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGagc                                                                                                                                             <  2:743907/82‑3 (MQ=38)
                                                                          gTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGagc                                                                                                                                             >  1:743907/1‑80 (MQ=39)
                                                                          gTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAgcg                                                                         >  2:1602774/1‑150 (MQ=255)
                                                                          gTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAg                                                                           <  2:1658744/148‑1 (MQ=255)
                                                                          gTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAg                                                                           >  1:1658744/1‑148 (MQ=255)
                                                                           tAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAgcgc                                                                        >  1:231494/1‑150 (MQ=255)
                                                                              agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                     <  1:148789/150‑1 (MQ=255)
                                                                              agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGACTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                                     >  2:268135/1‑150 (MQ=255)
                                                                                gCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTa                                                                   >  2:1074752/1‑150 (MQ=255)
                                                                                  ccATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGt                                                                 >  1:1866541/1‑150 (MQ=255)
                                                                                      gtgGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTg                                                             <  1:419658/150‑1 (MQ=255)
                                                                                         gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                   >  2:1022756/1‑109 (MQ=255)
                                                                                         gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                                   <  1:1022756/109‑1 (MQ=255)
                                                                                         gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTgcgc                                                          >  1:1105316/1‑150 (MQ=255)
                                                                                             tAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgt                                                      <  1:251731/150‑1 (MQ=255)
                                                                                              aaTAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgtg                                                     >  1:266075/1‑150 (MQ=255)
                                                                                                       gcagcaGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtg                                            <  1:169070/150‑1 (MQ=255)
                                                                                                       gcagcaGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtg                                            <  1:215151/150‑1 (MQ=255)
                                                                                                        cagcagGAACAGCCTCGCGCTTGGCTGCTCGGGTAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtgg                                           <  1:274959/150‑1 (MQ=255)
                                                                                                         agcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggc                                          >  1:688670/1‑150 (MQ=255)
                                                                                                           cagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctg                                        >  2:1737496/1‑150 (MQ=255)
                                                                                                            agGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctgg                                       <  2:914249/150‑1 (MQ=255)
                                                                                                             ggAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctggc                                      >  1:196145/1‑150 (MQ=255)
                                                                                                              gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCg                                     >  1:634375/1‑150 (MQ=255)
                                                                                                              gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCg                                     <  2:1134370/150‑1 (MQ=255)
                                                                                                               aaCAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTTCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGa                                    >  2:1012778/1‑150 (MQ=255)
                                                                                                                aCAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAgcgc                                                                        >  2:99580/1‑113 (MQ=255)
                                                                                                                aCAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAgcgc                                                                        <  1:99580/113‑1 (MQ=255)
                                                                                                                  aGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtg                                 <  2:1105316/150‑1 (MQ=255)
                                                                                                                        gcgcTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTATTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAAc                           >  1:1553128/1‑150 (MQ=255)
                                                                                                                         gttcgacctggtCGGGCAGGCATGGCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACg                          >  2:198791/8‑150 (MQ=255)
                                                                                                                            ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgcg                       >  2:1802521/1‑150 (MQ=255)
                                                                                                                            ttGGCTGCTCGGGCAGCCATGCCTGCCAGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgcg                       >  2:741234/1‑150 (MQ=255)
                                                                                                                              ggCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTCCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGcc                     >  1:1307079/1‑150 (MQ=255)
                                                                                                                                    tCGGGCAGGCATGGCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg               >  2:1630511/1‑150 (MQ=255)
                                                                                                                                    tCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctg                                        >  2:566141/1‑125 (MQ=255)
                                                                                                                                    tCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctg                                        <  1:566141/125‑1 (MQ=255)
                                                                                                                                    tCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg               >  2:1516483/1‑150 (MQ=255)
                                                                                                                                      gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAgg                                                                                    >  1:768306/1‑79 (MQ=255)
                                                                                                                                      gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAgg                                                                                    <  2:768306/79‑1 (MQ=255)
                                                                                                                                      gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGTCTGGCGAgtg                                 <  2:477636/130‑1 (MQ=255)
                                                                                                                                      gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtg                                 >  1:477636/1‑130 (MQ=255)
                                                                                                                                           gcgaTGCCTGCCCGACCAGGTCGAACCCCGACGACCGCCAGGTGCGCGGTGAGCTGTGTGCCGGCAGGTTCAGCGCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCa        <  1:1398952/147‑1 (MQ=255)
                                                                                                                                            ccATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg       <  2:1180083/150‑1 (MQ=255)
                                                                                                                                              aTGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtg                                            <  1:1800806/111‑1 (MQ=255)
                                                                                                                                              aTGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtg                                            >  2:1800806/1‑111 (MQ=255)
                                                                                                                                              aTGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGc     >  2:1907106/1‑150 (MQ=255)
                                                                                                                                                 cctgccCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGTGGGTGCGCAGGCGGt  >  2:1045995/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAGGCGGT  >  NC_002947/5886680‑5886974

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: ‑