Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,665,249 T→G 5.7% I207L (ATC→CTC)  PP_1459 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,665,2490TG 5.7% 201.4 / 2.0 66I207L (ATC→CTC) PP_1459CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base G (2/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GCGCTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGTC  >  NC_002947/1665104‑1665393
                                                                                                                                                 |                                                                                                                                                   
gcgcTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCg                                                                                                                                                  >  6:205040/1‑149 (MQ=255)
        cAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCg                                                                                                                                          >  7:288860/1‑149 (MQ=255)
          gcgcTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGc                                                                                                                                        >  1:39452/1‑149 (MQ=255)
          gcgcTTCCAGGGTCAGCTCGACCTTGGGCAGCAGGTGGCTGATCAGCTTGGGGGTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGc                                                                                                                                        <  8:147429/149‑1 (MQ=255)
            gcTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCtc                                                                                                                                      <  5:61072/149‑1 (MQ=255)
              ttCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCtcgt                                                                                                                                    >  3:178651/1‑149 (MQ=255)
              ttCCAGGGTCAGCTCGGCCTTGGGCAGCAGCTCGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTAGACGATGCGTTCGATCGGCtcgt                                                                                                                                    >  5:115386/1‑149 (MQ=255)
                 cAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGa                                                                                                                                 <  3:81305/149‑1 (MQ=255)
                  cctggaaGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAg                                                                                                                                <  2:325564/143‑1 (MQ=255)
                  cctggaaGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAg                                                                                                                                <  2:331041/143‑1 (MQ=255)
                   gggTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAgg                                                                                                                               <  1:299186/149‑1 (MQ=255)
                      tCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAgcagc                                                                                                                            <  4:39345/149‑5 (MQ=255)
                         gCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATg                                                                                                                         <  2:39452/149‑1 (MQ=255)
                          cTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGc                                                                                                                        >  3:82575/1‑149 (MQ=255)
                             ggCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGcc                                                                                                                       <  7:205800/147‑1 (MQ=255)
                              gCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGGCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCATTCGATCGGCTCGTCGAGGTTGATGCCATc                                                                                                                    <  6:249239/149‑1 (MQ=255)
                               ccTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATc                                                                                                                    >  5:2342/1‑148 (MQ=255)
                                   gggCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCAc                                                                                                                 <  5:205040/147‑1 (MQ=255)
                                    ggCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTc                                                                                                              <  6:2342/149‑1 (MQ=255)
                                      cagctgGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGt                                                                                                            <  7:298479/149‑1 (MQ=255)
                                           gttcgCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGGTGATGGCGTTGTGGTAGACCGGCCGGCGGGTGTGGCGGCTGACGATCAGCTGCGCGCCGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTgcg‑c                                                                                                      <  1:380121/145‑1 (MQ=255)
                                                 tGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGGCGGTCGATCGGCTCGGCGAGGTTGATGCCATCCACTTCGTAGcg‑cgcgatc                                                                                                >  2:380121/1‑144 (MQ=255)
                                                            gtgtTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGtt                                                                                   <  8:288860/149‑1 (MQ=255)
                                                              gtTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGa                                                                                 >  6:361600/1‑149 (MQ=255)
                                                               ttGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGAc                                                                                >  7:79379/1‑149 (MQ=255)
                                                                 gAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCAc                                                                                                 <  5:135790/132‑1 (MQ=255)
                                                                 gAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCAc                                                                                                 >  6:135790/1‑132 (MQ=255)
                                                                 gAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACg                                                                               <  4:49999/148‑1 (MQ=255)
                                                                   ggATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGt                                                                             >  5:109402/1‑148 (MQ=255)
                                                                   ggATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCGGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTg                                                                            <  4:31509/149‑1 (MQ=255)
                                                                     aTCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGAt                                                                          >  5:290626/1‑149 (MQ=255)
                                                                      tCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATc                                                                         <  6:109402/149‑1 (MQ=255)
                                                                       cGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATc                                                                         <  4:82575/148‑1 (MQ=255)
                                                                       cGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATc                                                                         <  6:115386/148‑1 (MQ=255)
                                                                          ttCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGt                                                                     >  3:290759/1‑149 (MQ=255)
                                                                          ttCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGt                                                                     >  2:138995/1‑149 (MQ=255)
                                                                                cTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGc                                                                >  8:266265/1‑148 (MQ=255)
                                                                                        tGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAgg                                                       >  5:46937/1‑149 (MQ=255)
                                                                                         gTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGc                                                      >  2:108132/1‑149 (MQ=255)
                                                                                          tCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGcc                                                     >  7:74148/1‑149 (MQ=255)
                                                                                          tCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGcc                                                     >  3:213642/1‑149 (MQ=255)
                                                                                            gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGt                                                   >  6:219580/1‑149 (MQ=255)
                                                                                            gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGATAGTAGAGGTTGATGCCATCCACTTAGTTGCG‑CGGCAACAGGAT‑G‑TCGTTGACGGTGATCTTGTCAAGCGACTGCAGGCCGt                                                   >  7:28448/1‑149 (MQ=255)
                                                                                             ttGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTc                                                  >  8:34276/1‑149 (MQ=255)
                                                                                             ttGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTc                                                  <  4:290759/149‑1 (MQ=255)
                                                                                             ttGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTc                                                  <  3:44085/149‑1 (MQ=255)
                                                                                             ttGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGAG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTc                                                  >  7:17039/1‑149 (MQ=255)
                                                                                                tgGTAGACCGGCAGGCGTGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTTATGCCATCCACTTCGTTGCGTCGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCg                                                 <  1:283091/148‑1 (MQ=255)
                                                                                                 ggTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAg                                               >  5:243413/1‑148 (MQ=255)
                                                                                                      aCCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGa                                           >  4:265885/1‑147 (MQ=255)
                                                                                                      aCCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGa                                           <  3:265885/147‑1 (MQ=255)
                                                                                                      aCCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGAcc                                         >  5:59374/1‑149 (MQ=255)
                                                                                                        cGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCt                                       >  5:4431/1‑149 (MQ=255)
                                                                                                        cGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCGCGACGATGCGTGCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCt                                       <  6:270167/149‑1 (MQ=255)
                                                                                                         ggCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTg                                      <  6:59374/149‑1 (MQ=255)
                                                                                                             ggCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGa                                  <  8:74148/149‑1 (MQ=255)
                                                                                                                        cggcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTg                       <  5:219580/149‑1 (MQ=255)
                                                                                                                                 aTCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑GCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATc              <  6:46937/149‑1 (MQ=255)
                                                                                                                                     gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCCTCCTCa          >  5:379317/1‑149 (MQ=255)
                                                                                                                                     gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTc           >  2:344132/1‑148 (MQ=255)
                                                                                                                                         ctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCCGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCg       >  2:263554/1‑148 (MQ=255)
                                                                                                                                         ctCGACGAGCCGATCGAACGCATCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCg       >  7:355005/1‑148 (MQ=255)
                                                                                                                                          tCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGATGCGCTTCGTGCTGTTCATACTCATCg       >  3:219018/1‑147 (MQ=255)
                                                                                                                                          tCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCg       >  6:62339/1‑147 (MQ=255)
                                                                                                                                            gacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGaa     <  8:59385/147‑1 (MQ=255)
                                                                                                                                               gaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGATAGTTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAgtc  >  1:297805/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGCTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCG‑CGGCACCAGGAT‑G‑TCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGTC  >  NC_002947/1665104‑1665393

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