Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,873,945 G→C 5.5% P198P (CCG→CCC PP_5147 → haloacid dehalogenase‑like family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,873,9450GC 5.5% 229.5 / 3.8 73P198P (CCG→CCCPP_5147haloacid dehalogenase‑like family hydrolase
Reads supporting (aligned to +/- strand):  ref base G (37/32);  new base C (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CCGTGTTTCCGTGAAGGCAAAGTGACGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGGGCTGCGCTGTACCTGTAG  >  NC_002947/5873799‑5874092
                                                                                                                                                  |                                                                                                                                                    
ccGTGTTTCCGTGAAGGCAAAGTGACGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCC                                                                                                                                                   >  2:650996‑M1/1‑150 (MQ=255)
 cGTGTTTCCGTGAAGGCAAAGTGACGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCTGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGccc                                                                                                                                                  >  2:1028699/1‑150 (MQ=255)
          gTGAAGGCAAAGTGACGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGc                                                                                                                                         <  1:1267764/150‑1 (MQ=255)
          gTGAAGGCAAAGTGACGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTATACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACTCGGCCCCGGACCTGc                                                                                                                                         >  2:177202/1‑150 (MQ=255)
                   aaGTGACGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAgg                                                                                                                               >  1:650996‑M2/1‑148 (MQ=255)
                    aGTGACGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGAAAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCtcggcccgcaggtccggggc                                                                                                                               <  2:751848‑M1/150‑21 (MQ=255)
                      tGACGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCgg                                                                                                                            <  1:436376/150‑1 (MQ=255)
                        aCGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCtcggcccgcaggtccggggccggg                                                                                                                           <  1:1268501‑M1/150‑25 (MQ=255)
                                  aGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCgtgggt                                                                                                                >  2:1070357/1‑150 (MQ=255)
                                           tGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGt                                                                                                         >  1:521616/1‑148 (MQ=255)
                                           tGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGt                                                                                                         <  2:521616/148‑1 (MQ=255)
                                             cTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGAtc                                                                                                     <  1:1079468/150‑1 (MQ=255)
                                                        cTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCg                                                                                                           <  2:648202/133‑1 (MQ=255)
                                                        cTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCg                                                                                                           >  1:648202/1‑133 (MQ=255)
                                                         tttGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGc                                                                                                                                         <  2:787167/103‑1 (MQ=255)
                                                         tttGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGc                                                                                                                                         >  1:787167/1‑103 (MQ=255)
                                                         tttGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGAt                                                                                         >  1:453894/1‑150 (MQ=255)
                                                          ttGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGAtt                                                                                        <  1:510957/150‑1 (MQ=255)
                                                             aTCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAg                                                                                     <  2:453894/150‑1 (MQ=255)
                                                             aTCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAg                                                                                     <  1:1070357/150‑1 (MQ=255)
                                                              tCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAgg                                                                                    >  2:848602/1‑150 (MQ=255)
                                                                ttGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCCGCCtcggcccgcaggtccggggccgg                                                                                                                            <  2:1346844‑M1/109‑24 (MQ=255)
                                                                ttGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCCGCCtcggcccgcaggtccggggccgg                                                                                                                            >  1:1346844‑M1/1‑86 (MQ=255)
                                                                    ggACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTcc                                                                              >  2:776381/1‑150 (MQ=255)
                                                                              ttCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTc                                                                    >  2:858139/1‑150 (MQ=255)
                                                                                     gCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCtttt                                                                >  1:979882/1‑147 (MQ=255)
                                                                                     gCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCtttt                                                                <  2:979882/147‑1 (MQ=255)
                                                                                      cGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGGTTTCCCGTGCCGGTCTTTTcgcg                                                            >  2:1151793/1‑150 (MQ=255)
                                                                                       gACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGAGTGCGCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCgg                                                           <  1:812954/150‑1 (MQ=255)
                                                                                        aCTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCtcggcccgcaggtccggggccgg                                                                                                                            >  2:1022962‑M1/1‑62 (MQ=255)
                                                                                        aCTCAATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCtcggcccgcaggtccggggccgg                                                                                                                            <  1:1022962‑M1/85‑24 (MQ=255)
                                                                                          tCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCa                                                        >  1:986330/1‑150 (MQ=255)
                                                                                            gATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCTCGGCCCGCAGGTCCGGggccgggtccactgc                                                                                                                  >  2:1124228‑M2/1‑75 (MQ=255)
                                                                                            gATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCTCGGCCCGCAGGTCCGGggccgggtccactgc                                                                                                                  <  1:1124228‑M2/90‑16 (MQ=255)
                                                                                            gATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGCCTATTCGCgggtgag                                                      >  2:109987/1‑146 (MQ=255)
                                                                                             atgaatgaCCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGa                                                                                                                                 <  1:689694/74‑1 (MQ=255)
                                                                                             atgaatgaCCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGa                                                                                                                                 >  2:689694/1‑74 (MQ=255)
                                                                                             atgaatgaCCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAggc                                                                                                                              >  1:283867‑M2/1‑74 (MQ=255)
                                                                                             atgaatgaCCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAggc                                                                                                                              <  2:283867‑M2/77‑4 (MQ=255)
                                                                                                    atcttccGTTGTTGCACCGGGTGACGCATGCAGTGGAAGTGGAACCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCCGTCTTTTCGCGGGCACGCCCGCTcc                                              <  2:168366/145‑1 (MQ=255)
                                                                                                     ccTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTccc                                             >  2:1238809/1‑150 (MQ=255)
                                                                                                      cTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTGTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCca                                            >  1:1008988/1‑150 (MQ=255)
                                                                                                      cTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCCGCCtcggcccgcaggtccggggccgg                                                                                                                            <  2:31677‑M1/71‑24 (MQ=255)
                                                                                                      cTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCCGCCtcggcccgcaggtccggggccgg                                                                                                                            >  1:31677‑M1/1‑48 (MQ=255)
                                                                                                       tGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCcac                                           >  1:699169/1‑150 (MQ=255)
                                                                                                           gttgttGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGa                                                                                      >  1:496628/1‑103 (MQ=255)
                                                                                                           gttgttGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGa                                                                                      <  2:496628/103‑1 (MQ=255)
                                                                                                                tgCAGCGGGTGACGCATGCAGTGGCAGTCGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATcac                                  <  1:407145/150‑1 (MQ=255)
                                                                                                                   aGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCACTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATcaccac                               >  1:252400/1‑150 (MQ=255)
                                                                                                                             gCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGAGAGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGACCACACCAACAAGAATCGCCACAGGCCTTGAg                     >  1:385838/1‑150 (MQ=255)
                                                                                                                              cATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAgg                    <  1:858139/150‑1 (MQ=255)
                                                                                                                                tGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAg                                                                                     <  2:1199039/83‑1 (MQ=255)
                                                                                                                                tGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAg                                                                                     >  1:1199039/1‑83 (MQ=255)
                                                                                                                                 gcagtggcagtgGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGGGCt                 >  1:1401084/1‑150 (MQ=255)
                                                                                                                                   agtggcaggggACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGGGCTgc               <  1:177202/150‑1 (MQ=255)
                                                                                                                                       gcagtggACCCGGCCccggacctgcgggccgaggccgggtccactgccactgcatgcgtcacccgctgcaacaacggcaggtcattcatc                                                                        <  1:1124228‑M1/90‑76 (MQ=255)
                                                                                                                                       gcagtggACCCGGCCccggacctgcgggccgaggccgggtccactgccactgcatgcgtcacccgctgcaacaacggcaggtcattcatc                                                                        >  2:1124228‑M1/1‑15 (MQ=255)
                                                                                                                                        cagtggACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCgg                                                           <  1:585877/101‑1 (MQ=255)
                                                                                                                                        cagtggACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCgg                                                           >  2:585877/1‑101 (MQ=255)
                                                                                                                                         agtggACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGGGCTGCGCTGTa         >  1:626064/1‑150 (MQ=255)
                                                                                                                                             gACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGGGCTGCGCTGTACCTg     >  2:790820/1‑150 (MQ=255)
                                                                                                                                               cccGGCCCCGGACCTGCGGGCCGAggccgggtccactgccactgcatgcgtcacccgctgcaacaacggcaggtcattcatcgagtcgctgtagaagtagctgtcctcaagatcaaagccgttctccagcatccagcgctccagccgcgt   <  1:1268501‑M2/150‑127 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAggcggggtccactgccactgcatgcgtcacccgctgcaacaacggcaggtcattcatcgagtcgctgtagaagtagctgtcctcaa                                           >  1:1346844‑M2/1‑23 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAggcggggtccactgccactgcatgcgtcacccgctgcaacaacggcaggtcattcatcgagtcgctgtagaagtagctgtcctcaa                                           <  2:1346844‑M2/109‑87 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAggcggggtccactgccactgcatgcgtcacccgctgcaacaacggcag                                                                                 <  2:31677‑M2/71‑49 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAggcggggtccactgccactgcatgcgtcacccgctgcaacaacggcag                                                                                 >  1:31677‑M2/1‑23 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAggccgggtccactgccactgcatgcgtcacccgctgcaacaacggcaggtcattcattgagt                                                                   <  1:1022962‑M2/85‑63 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAggccgggtccactgccactgcatgcgtcacccgctgcaacaacggcaggtcattcatcgagt                                                                   >  2:1022962‑M2/1‑23 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCgtg                                                                                                                   >  1:824070/1‑37 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCgtg                                                                                                                   <  2:824070/37‑1 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGt                                                                                                         <  1:88505/47‑1 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGt                                                                                                         >  2:88505/1‑47 (MQ=255)
                                                                                                                                                ccGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCACTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGGGCTGCGCTGTACCTGTAg  <  2:252400/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCGTGTTTCCGTGAAGGCAAAGTGACGCGGCTGGAGCGCTGGATGCTGGAGAACGGCTTTGATCTTGAGGACAGCTACTTCTACAGCGACTCGATGAATGACCTGCCGTTGTTGCAGCGGGTGACGCATGCAGTGGCAGTGGACCCGGCCCCGGACCTGCGGGCCGAGGCGGATAAGCGTGGGTGGCCGGTGATCTCGCTGCGGGATTGAGGTTTTCCCGTGCCGGTCTTTTCGCGGGCACGCCCGCTCCCACAGGAATCACCACAGGCCTTGAGGGCTGCGCTGTACCTGTAG  >  NC_002947/5873799‑5874092

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: