Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,586,523 T→C 6.2% T255A (ACC→GCC)  PP_4914 ← phytase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,586,5230TC 6.2% 189.8 / 2.1 64T255A (ACC→GCC) PP_4914phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (2/2);  total (30/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 = 9.98e-01

ATCGATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCT  >  NC_002947/5586377‑5586669
                                                                                                                                                  |                                                                                                                                                  
aTCGATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTaa                                                                                                                                                  <  7:287091/149‑1 (MQ=255)
aTCGATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTa                                                                                                                                                   >  1:24488/1‑148 (MQ=255)
aTCGATCACCAGTGCTGCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTaa                                                                                                                                                  <  8:314316/149‑1 (MQ=255)
   gATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAAc                                                                                                                                                <  2:24488/148‑1 (MQ=255)
                 cGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCTACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCGaa                                                                                                                                  <  8:247016/148‑1 (MQ=255)
                    ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccag                                                                                                                              <  7:121995/149‑1 (MQ=255)
                      aggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                            <  4:361764/149‑1 (MQ=255)
                      aggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                            <  3:179190/149‑1 (MQ=255)
                         ccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggccag                                                                                                                         <  3:50647/149‑1 (MQ=255)
                                 cTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGa                                                                                                                 <  7:306747/149‑1 (MQ=255)
                                   gcgcGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCGAAc                                                                                                                                 <  6:197597/131‑1 (MQ=255)
                                   gcgcGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCGAAc                                                                                                                                 >  5:197597/1‑131 (MQ=255)
                                      cgTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTcg                                                                                                            >  7:49853/1‑149 (MQ=255)
                                      cgTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTcg                                                                                                            >  5:242266/1‑149 (MQ=255)
                                          gcgcCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGaa                                                                                                                <  6:230943/141‑1 (MQ=255)
                                          gcgcCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGaa                                                                                                                >  5:230943/1‑141 (MQ=255)
                                            gcCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGa                                                                                                       >  2:206033/1‑148 (MQ=255)
                                               gTGACTCGACCAGCGCATCGGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAgg                                                                                                   >  2:366614/1‑149 (MQ=255)
                                                   cTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTc                                                                                                                                     >  3:254880/1‑111 (MQ=255)
                                                   cTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTc                                                                                                                                     <  4:254880/111‑1 (MQ=255)
                                                         cAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGaa                                                                                                                                  <  8:257229/108‑1 (MQ=255)
                                                         cAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGaa                                                                                                                                  >  7:257229/1‑108 (MQ=255)
                                                          aGCACATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCAtt                                                                                        <  3:107815/149‑1 (MQ=255)
                                                            cgcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTgg                                                                                      <  6:242266/149‑1 (MQ=255)
                                                             gcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                     <  7:160120/149‑1 (MQ=255)
                                                               aTCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCCATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACCGGGCCCTGGCATTGGGCa                                                                                   >  1:362802/1‑149 (MQ=255)
                                                                 cagcagCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATc                                                                                 >  6:23459/1‑149 (MQ=255)
                                                                 cagcagCAAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATc                                                                                 >  8:316397/1‑149 (MQ=255)
                                                                           cGACCAGCAACGCTCGACCTTGCCACTGTGGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGtt                                                                       <  5:47633/149‑1 (MQ=255)
                                                                               cAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                     <  1:35952/131‑1 (MQ=255)
                                                                               cAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                     >  2:35952/1‑131 (MQ=255)
                                                                                aGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTg                                                                      >  6:286193/1‑145 (MQ=255)
                                                                                aGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTg                                                                      <  5:286193/145‑1 (MQ=255)
                                                                                 gCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAttt                                                                 >  6:54106/1‑149 (MQ=255)
                                                                                    aCGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGcgcg                                                                           <  6:196961/136‑1 (MQ=255)
                                                                                    aCGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGcgcg                                                                           >  5:196961/1‑136 (MQ=255)
                                                                                        tCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGtgcg                                                                           <  7:147089/132‑1 (MQ=255)
                                                                                        tCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGtgcg                                                                           >  8:147089/1‑132 (MQ=255)
                                                                                               tGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCg                                                   >  3:102834/1‑149 (MQ=255)
                                                                                                 ccACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGc                                                  <  1:206033/148‑1 (MQ=255)
                                                                                                   aCTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTca                                               <  4:102834/149‑1 (MQ=255)
                                                                                                     tGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGcc                                                                                                >  5:377482/1‑98 (MQ=255)
                                                                                                     tGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGcc                                                                                                <  6:377482/98‑1 (MQ=255)
                                                                                                     tGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagc                                             <  7:316397/149‑1 (MQ=255)
                                                                                                           tatCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACg                                       >  3:14842/1‑149 (MQ=255)
                                                                                                              cACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagc                                             >  1:197750/1‑140 (MQ=255)
                                                                                                              cACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagc                                             <  2:197750/140‑1 (MQ=255)
                                                                                                               aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCaa                                    <  4:356790/148‑1 (MQ=255)
                                                                                                               aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCaa                                    >  3:356790/1‑148 (MQ=255)
                                                                                                               aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAAc                                   <  1:366614/149‑1 (MQ=255)
                                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGcc                                                                                                >  8:292996/1‑82 (MQ=255)
                                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGcc                                                                                                <  7:292996/82‑1 (MQ=255)
                                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTca                                               >  6:114176/1‑131 (MQ=255)
                                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTca                                               <  5:114176/131‑1 (MQ=255)
                                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGcc                             >  8:9322/1‑149 (MQ=255)
                                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGcc                             >  8:3765/1‑149 (MQ=255)
                                                                                                                                 ggCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAg                  <  8:49853/148‑1 (MQ=255)
                                                                                                                                  gCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGAc                <  5:23459/149‑1 (MQ=255)
                                                                                                                                   cGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTgcg                                                                             <  3:51982/87‑1 (MQ=255)
                                                                                                                                   cGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTgcg                                                                             >  4:51982/1‑87 (MQ=255)
                                                                                                                                        cGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGAct           <  5:54106/148‑1 (MQ=255)
                                                                                                                                           tCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActgctg       >  1:241091/1‑149 (MQ=255)
                                                                                                                                            cGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTAAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTgg      >  2:331921/1‑149 (MQ=255)
                                                                                                                                                ggTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCt  >  5:224629/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATCGATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCT  >  NC_002947/5586377‑5586669

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