Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,886,820 C→G 7.9% W84C (TGG→TGC PP_5161 ← lysophospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,886,8200CG 7.9% 273.4 / 12.9 88W84C (TGG→TGCPP_5161lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (44/37);  new base G (4/3);  total (48/40)
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.02e-01

TCGCCGCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAG  >  NC_002947/5886671‑5886969
                                                                                                                                                     |                                                                                                                                                     
tCGCCGCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTTCTCGGGCAGc                                                                                                                                                       >  2:215151/1‑150 (MQ=255)
  gccgcTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCa                                                                                                                                                     >  1:1134370/1‑150 (MQ=255)
   ccgcTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt                                                                                                                                                    <  1:111890/150‑1 (MQ=255)
         gACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgcctgc                                                                                                                                              >  2:1058780/1‑150 (MQ=255)
            aGGTCGTGGCCGGGCAGTTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt                                                                                                                                                    <  2:939783/141‑1 (MQ=255)
            aGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt                                                                                                                                                    >  1:939783/1‑141 (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:787437/150‑1 (MQ=255)
                                acaGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGc                                                                                                                       <  1:154941/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                                                                                                                                <  1:860296/133‑1 (MQ=255)
                                        tCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGaa                                                                                                                                >  2:860296/1‑133 (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                                                                                                <  2:281754/140‑1 (MQ=255)
                                                                 ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                                >  1:281754/1‑140 (MQ=255)
                                                                 ccACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAg                                                                                      >  2:169070/1‑150 (MQ=255)
                                                                   aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGccc                                                                                                                                            <  1:1072330/94‑1 (MQ=255)
                                                                   aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGccc                                                                                                                                            >  2:1072330/1‑94 (MQ=255)
                                                                   aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCCCGAGCAGCCaagc                                                                                                                                                  >  2:152933/1‑85 (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                                                                      >  1:1658744/1‑148 (MQ=255)
                                                                                   gTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAg                                                                      <  2:1658744/148‑1 (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                                                                                              <  1:1022756/109‑1 (MQ=255)
                                                                                                  gTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                              >  2:1022756/1‑109 (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:215151/150‑1 (MQ=255)
                                                                                                                gcagcaGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtg                                       <  1:169070/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                                <  2:1134370/150‑1 (MQ=255)
                                                                                                                       gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCg                                >  1:634375/1‑150 (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                                                                               <  2:768306/79‑1 (MQ=255)
                                                                                                                                               gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAgg                                                                               >  1:768306/1‑79 (MQ=255)
                                                                                                                                               gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGTCTGGCGAgtg                            <  2:477636/130‑1 (MQ=255)
                                                                                                                                               gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtg                            >  1:477636/1‑130 (MQ=255)
                                                                                                                                                     ccATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAg  <  2:1180083/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCCGCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCGCAG  >  NC_002947/5886671‑5886969

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