Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,783,969 G→C 10.6% E127Q (GAG→CAG)  PP_2437 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,783,9690GC10.6% 201.9 / 14.5 66E127Q (GAG→CAG) PP_2437acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base C (4/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

GCCTGAGCAACCTGGAGTACGCACCGCTGGCAGAGATCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTGGGTGATCAACGGCCGCAAGTGGTGGACGTC  >  NC_002947/2783820‑2784111
                                                                                                                                                      |                                                                                                                                              
gCCTGAGCAACCTGGAGTACGCACCGCTGGCAGAGATCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTgg                                                                                                                                                >  1:268181/1‑150 (MQ=255)
 ccTGAGCAACCTGGAGTACGCACCGCTGGCAGAGATCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGa                                                                                                                                               <  2:706733/150‑1 (MQ=255)
   tGAGCAACCTGGAGTGCGAACAGATCGCGAATGTCATGCCCCTCTCGTTGCTCGGGCCGGATGGTTTTCACGTCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGc                                                                                                                                             <  1:21786/150‑1 (MQ=255)
      gCAACCTGGAGTACGCACCGCTGGCAGAGATCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCAc                                                                                                                                          <  2:276099/150‑1 (MQ=255)
         aCCTGGAGTACGCACCGCTGGCAGAGATCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCActgc                                                                                                                                       >  1:364721/1‑150 (MQ=255)
          ccTGGAGTACGCACCGCTGGCAGAGATCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCActgct                                                                                                                                      >  1:618504/1‑150 (MQ=255)
            tGGAGTACGCACCGCTGGCAGAGATCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACtgctgc                                                                                                                                    >  1:956050/1‑150 (MQ=255)
                 tACGCACCGCTGGCAGAGATCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCgcggc                                                                                                                               >  2:834172/1‑150 (MQ=255)
                                  gaTCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTTGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGtt                                                                                                              <  1:834172/150‑1 (MQ=255)
                                    tCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTg                                                                                                            <  2:850426/150‑1 (MQ=255)
                                      aTGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGcc                                                                                                          <  2:877711/150‑1 (MQ=255)
                                          gCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTCCAGCCActgg                                                                                                                                       >  2:282411/1‑116 (MQ=255)
                                          gCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTCCAGCCActgg                                                                                                                                       <  1:282411/117‑2 (MQ=255)
                                           ccGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGAc                                                                                                     >  1:262834/1‑150 (MQ=255)
                                                cGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAAc                                                                                                <  2:956050/150‑1 (MQ=255)
                                                      tCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATc                                                                                                                         >  2:313359/1‑119 (MQ=255)
                                                      tCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATc                                                                                                                         <  1:313359/119‑1 (MQ=255)
                                                      tCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTCCAGCCActgg                                                                                                                                       <  2:237215/105‑2 (MQ=255)
                                                      tCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTCCAGCCActgg                                                                                                                                       >  1:237215/1‑104 (MQ=255)
                                                        gggCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGc                                                                                                           <  2:716362/131‑1 (MQ=255)
                                                        gggCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGc                                                                                                           >  1:716362/1‑131 (MQ=255)
                                                        gggCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTg                                                                                        <  2:268181/150‑1 (MQ=255)
                                                        gggCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTg                                                                                        <  2:364721/150‑1 (MQ=255)
                                                          gCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATCACCGAACCGGACGTGGc                                                                                      <  2:618504/150‑1 (MQ=255)
                                                            ccgAGCAGTTGAACGGCTCCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCgcgg                                                                                                                                <  2:338392/105‑1 (MQ=255)
                                                                gCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCctcctc                                                                                <  2:262834/150‑1 (MQ=255)
                                                                   gTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGa                                                                             >  2:1158002/1‑150 (MQ=255)
                                                                    ttCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGAc                                                                            >  1:301638/1‑150 (MQ=255)
                                                                    ctCAACTGCTCGGGGGCGGACACTGGCAAAATGGAAGTGCTGGTGCGCTAGGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGAc                                                                            >  2:111929/2‑150 (MQ=255)
                                                                        aCTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGcca                                                                        <  2:1127188/150‑1 (MQ=255)
                                                                            cTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGAc                                                                                                     >  2:350009/1‑117 (MQ=255)
                                                                            cTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGAc                                                                                                     <  1:350009/117‑1 (MQ=255)
                                                                             tCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTACGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAAc                                                                   >  2:190054/1‑150 (MQ=255)
                                                                                ggCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCgaga                                                                                                                           >  1:5622/1‑92 (MQ=255)
                                                                                ggCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCgaga                                                                                                                           <  2:5622/92‑1 (MQ=255)
                                                                                ggCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACAt                                                                 >  1:729506/1‑150 (MQ=255)
                                                                                 gcgcCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATg                                                                <  1:1158002/150‑1 (MQ=255)
                                                                                      ggACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCgaga                                                                                                                           <  2:1049241/86‑1 (MQ=255)
                                                                                      ggACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCgaga                                                                                                                           >  1:1049241/1‑86 (MQ=255)
                                                                                      ggACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATtgctgc                                                           <  2:729506/150‑1 (MQ=255)
                                                                                      ggACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGgctgc                                                           >  1:1017696/1‑150 (MQ=255)
                                                                                            tGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGc                                                     >  1:232219/1‑150 (MQ=255)
                                                                                               cAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAAc                                                                                                <  1:1089920/104‑1 (MQ=255)
                                                                                               cAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAAc                                                                                                >  2:1089920/1‑104 (MQ=255)
                                                                                                            cTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCTAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGAc                                     >  2:436040/1‑150 (MQ=255)
                                                                                                             tGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGGCAAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGAcc                                    <  1:980532/150‑1 (MQ=255)
                                                                                                             tGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGAc                                                                            >  1:738621/1‑110 (MQ=255)
                                                                                                             tGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGAc                                                                            <  2:738621/110‑1 (MQ=255)
                                                                                                                  cgcTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTg                                <  1:457632/149‑1 (MQ=255)
                                                                                                                  cgcTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTg                                >  2:457632/1‑149 (MQ=255)
                                                                                                                   gcTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTggg                              <  2:232219/150‑1 (MQ=255)
                                                                                                                   gcTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTggg                              >  1:541328/1‑150 (MQ=255)
                                                                                                                             cGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTGGGTGATCAACgg                    <  1:196033/150‑1 (MQ=255)
                                                                                                                               aGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTGGGTGATCAACGGcc                  <  2:1017696/150‑1 (MQ=255)
                                                                                                                                 gcgcAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTAc                                                        >  2:805633/1‑110 (MQ=255)
                                                                                                                                 gcgcAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTAc                                                        <  1:805633/110‑1 (MQ=255)
                                                                                                                                           gCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACgg                                              <  2:410738/110‑1 (MQ=255)
                                                                                                                                           gCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACgg                                              >  1:410738/1‑110 (MQ=255)
                                                                                                                                           gCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGAc                                     >  1:319217/1‑119 (MQ=255)
                                                                                                                                           gCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGAc                                     <  2:319217/119‑1 (MQ=255)
                                                                                                                                           agcAGTGGCTCCAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGAc                                     >  1:1069950/3‑119 (MQ=255)
                                                                                                                                           agcAGTGGCTCCAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGAc                                     <  2:1069950/117‑1 (MQ=255)
                                                                                                                                            ccAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTGGGTGATCAACGGCCGCAAGTGGTGGAc     >  1:686957/1‑150 (MQ=255)
                                                                                                                                            ccAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTGGGTGATCAACGGCCGCAAGTGGTGGAc     >  1:831445/1‑150 (MQ=255)
                                                                                                                                             cAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTGGGTGATCAACGGCCGCAAGTGGTGGACg    >  2:14722/1‑150 (MQ=255)
                                                                                                                                             cAGTGGCTCCAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTGGGTGATCAAAGGAGTCAAGGGGGGGACg    >  2:634314/1‑150 (MQ=255)
                                                                                                                                               gTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTGGGTGATCAACGGCCGCAAGTGGTGGACGTc  >  1:101698/1‑150 (MQ=255)
                                                                                                                                                      gAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGgct                                                             >  2:717193/1‑84 (MQ=255)
                                                                                                                                                      gAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGgct                                                             <  1:717193/84‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
GCCTGAGCAACCTGGAGTACGCACCGCTGGCAGAGATCATGGGCCGCTCGTTGCTCGGGCCGGAGCCGTTCAACTGCTCGGCGCCGGACACTGGCAACATGGAAGTGCTGGTGCGCTATGGCAGCGAGGCGCAGAAACGCCAGTGGCTGGAGCCACTGCTGCGCGGCGAGATCCGCTCAGCGTTTGCCATGACCGAACCGGACGTGGCCTCCTCGGACGCCACCAACATGGCTGCTACCGCCATACGGGACGGTGACCAGTGGGTGATCAACGGCCGCAAGTGGTGGACGTC  >  NC_002947/2783820‑2784111

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