Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,665,253 T→A 12.2% E205D (GAA→GAT)  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,2530TA12.2% 157.4 / 11.7 64E205D (GAA→GAT) PP_1459CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base A (3/4);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

CGCTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGTCGTCGGCC  >  NC_002947/1665105‑1665400
                                                                                                                                                    |                                                                                                                                                   
cgcTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGt                                                                                                                                                     <  4:231927/149‑1 (MQ=255)
cgcTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGt                                                                                                                                                     <  2:274374/149‑1 (MQ=255)
        aGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACg                                                                                                                                              >  4:291057/1‑148 (MQ=255)
               ccAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCtcgtc                                                                                                                                       >  4:104763/1‑148 (MQ=255)
                cAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGa                                                                                                                                     <  3:291057/149‑1 (MQ=255)
                cAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGa                                                                                                                                     <  6:146076/149‑1 (MQ=255)
                 aGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCATCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGCCGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAg                                                                                                                                    >  1:32161/1‑149 (MQ=255)
                  gggTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAgc                                                                                                                                   <  3:157555/149‑2 (MQ=255)
                   ggTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGt                                                                                                                                  <  5:232181/149‑1 (MQ=255)
                       aGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGAt                                                                                                                              <  2:50728/149‑1 (MQ=255)
                       aGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAgcagctg                                                                                                                              <  6:67885/149‑7 (MQ=38)
                         cTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGc                                                                                                                            >  5:282505/1‑149 (MQ=255)
                             gCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATc                                                                                                                        <  6:282505/149‑1 (MQ=255)
                              ccTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATcc                                                                                                                       <  6:44849/149‑1 (MQ=255)
                                          ggTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTgcgc                                                                                                           <  3:104763/149‑1 (MQ=255)
                                               cTGATCAGCTTGGTGTGGAGGATCGCTTCGACCTGGTTGATGGCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCAc                                                                                                      <  7:45490/149‑1 (MQ=255)
                                                 gATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCAc                                                                                                      >  2:54419/1‑147 (MQ=255)
                                                       cTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGt                                                                                              >  7:276453/1‑149 (MQ=255)
                                                          ggTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGt                                                                                           >  7:252983/1‑149 (MQ=255)
                                                           gtgtTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACTATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGtt                                                                                          <  8:276453/149‑1 (MQ=255)
                                                                gAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACgg                                                                                     >  6:174856/1‑149 (MQ=255)
                                                                 aGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGt                                                                                    >  7:174271/1‑149 (MQ=255)
                                                                       tcTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCtt                                                                              <  7:195385/148‑1 (MQ=255)
                                                                         ttCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCtt                                                                              >  8:195385/1‑147 (MQ=255)
                                                                             aCCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAg                                                                        >  2:230446/1‑149 (MQ=255)
                                                                              ccTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                                       <  3:284223/149‑1 (MQ=255)
                                                                               cTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                                       <  8:21979/148‑1 (MQ=255)
                                                                               cTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                                       >  7:21979/1‑148 (MQ=255)
                                                                               cTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCg                                                                      >  4:229311/1‑149 (MQ=255)
                                                                                  gTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGAt                                                                                 >  5:133374/1‑135 (MQ=255)
                                                                                  gTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGAt                                                                                 <  6:133374/135‑1 (MQ=255)
                                                                                        gTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTg                                                                                   <  5:147389/127‑1 (MQ=255)
                                                                                        gTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTg                                                                                   >  6:147389/1‑127 (MQ=255)
                                                                                          cGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGcc                                                            <  8:174271/148‑1 (MQ=255)
                                                                                           gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGt                                                          <  8:252983/149‑1 (MQ=255)
                                                                                              gtgGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGt                                                                            >  4:47644/1‑128 (MQ=255)
                                                                                              gtgGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGt                                                                            <  3:47644/128‑1 (MQ=255)
                                                                                                 gTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAg                                                      >  7:286817/1‑147 (MQ=255)
                                                                                                       cGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCAt                                                                                                                         >  8:9693/1‑74 (MQ=255)
                                                                                                       cGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCAt                                                                                                                         <  7:9693/74‑1 (MQ=255)
                                                                                                         gCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTgg                                            >  7:262231/1‑149 (MQ=255)
                                                                                                              cGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTTGCGAgg                                       <  5:174856/149‑1 (MQ=255)
                                                                                                                ggTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgcg                                     <  1:230446/149‑1 (MQ=255)
                                                                                                                 gtgtGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgcgc                                    >  6:35232/1‑149 (MQ=255)
                                                                                                                    tgGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTc                                 <  1:54419/149‑1 (MQ=255)
                                                                                                                    tgGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTc                                 <  5:253692/149‑1 (MQ=255)
                                                                                                                         gcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGc                                           >  6:217222/1‑134 (MQ=255)
                                                                                                                         gcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGc                                           <  5:217222/134‑1 (MQ=255)
                                                                                                                          cTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCTTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTg                           <  8:286817/149‑1 (MQ=255)
                                                                                                                            gACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGt                          >  3:107649/1‑148 (MQ=255)
                                                                                                                                 tCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATcc                    >  2:160388/1‑149 (MQ=255)
                                                                                                                                    gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTg                              <  8:26785/136‑1 (MQ=255)
                                                                                                                                    gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTg                              >  7:26785/1‑136 (MQ=255)
                                                                                                                                    gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCa                 >  8:298609/1‑149 (MQ=255)
                                                                                                                                        ctCGACGATGAGTTCGATCTGCTCGTCGAGGTTGATGCCATCCTCTTCGTTGCGCCTCACCAGGATGGCGTTGTCGGTTATCTTGTAGAGCGTGTGCAGGCCGTAGAGCATACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATc               <  2:32161/147‑1 (MQ=255)
                                                                                                                                        ctCGACGAGCCGATCGAACGCATCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCg              >  1:12096/1‑148 (MQ=255)
                                                                                                                                         tCGACGAGCCGCTCGAACGCATCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGAACGAGATCTAGTCCAGCGACTCCAGGCCGTCTAGCAGACCCTGGCGAGGAGCTTCGTGCTGTTCATCGTCAACGa             >  2:252867/1‑148 (MQ=255)
                                                                                                                                         tCGACGAGCCGATCGAACGCATCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCCTCCTCATCg              >  8:299491/1‑147 (MQ=255)
                                                                                                                                          cgacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAg           >  8:54545/1‑149 (MQ=255)
                                                                                                                                          cgacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAg           >  5:66377/1‑149 (MQ=255)
                                                                                                                                          cgacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAg           >  6:288778/1‑149 (MQ=255)
                                                                                                                                           gacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAg           <  5:288778/148‑1 (MQ=255)
                                                                                                                                               aTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAgtcgtc      <  1:124547/149‑1 (MQ=255)
                                                                                                                                                tGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgc                                      <  6:305635/116‑1 (MQ=255)
                                                                                                                                                tGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgc                                      >  5:305635/1‑116 (MQ=255)
                                                                                                                                                   gTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGTCGTCGGcc  >  2:11228/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGTCGTCGGCC  >  NC_002947/1665105‑1665400

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