Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 794,511 A→C 5.9% E218A (GAG→GCG)  PP_0680 → ATP‑dependent protease‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947794,5110AC 5.9% 215.6 / 3.4 68E218A (GAG→GCG) PP_0680ATP‑dependent protease‑like protein
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base C (2/2);  total (36/32)
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.84e-01

GGAAAAGGACGTGGCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGAAAAGTACGCCGAAAACGCGGCGGTG  >  NC_002947/794365‑794659
                                                                                                                                                  |                                                                                                                                                    
ggAAAAGGACGTGGCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGCCATTGCCCTGCTTGAGGAgc                                                                                                                                                    >  8:345054/1‑149 (MQ=255)
 gAAAAGGACGTGGCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAgc                                                                                                                                                    <  5:229255/148‑1 (MQ=255)
  aaaaGGACGTGGCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAgcgc                                                                                                                                                  <  1:24221/149‑1 (MQ=255)
   aaaGGACGTGGCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGcc                                                                                                                                                 >  2:120847/1‑149 (MQ=255)
       gACGTGGCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCaa                                                                                                                                             >  2:3117/1‑149 (MQ=255)
         cGTGGCCTTGTACCGCGATGCCAGCAACGTCGGCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACg                                                                                                                                           <  5:110977/149‑1 (MQ=255)
           tGGCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGaa                                                                                                                                         <  4:319197/149‑1 (MQ=255)
             gCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGCGCTCCTCaagcagg                                                                                                                                        <  7:147163/148‑6 (MQ=255)
             gCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAag                                                                                                                                        <  1:41339/148‑1 (MQ=255)
               cTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGTGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGc                                                                                                                                     >  4:8898/1‑149 (MQ=255)
               cTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGCGCTCCTCaagcagg                                                                                                                                        <  2:330345/146‑6 (MQ=255)
               cTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGCGCTCCTCaagcagg                                                                                                                                        >  1:330345/1‑141 (MQ=255)
                   tACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTgg                                                                                                                                  <  5:381472/148‑1 (MQ=255)
                   tACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGc                                                                                                                                 <  3:262034/149‑1 (MQ=255)
                     ccGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTgg                                                                                                                                  <  5:154507/146‑1 (MQ=255)
                     ccGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTgg                                                                                                                                  >  6:154507/1‑146 (MQ=255)
                      cgcgATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAg                                                                                                                              >  2:361061/1‑149 (MQ=255)
                          aTGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGATGTATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGGGAGCAATTCCACGAGGTCATTGCCCTGCTTGAGGAGCGGCTCAACGAAGAGCTGGCCAGCGTg                                                                                                                          <  2:200653/149‑1 (MQ=255)
                              cAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCg                                                                                                                       >  6:190565/1‑148 (MQ=255)
                               aGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGc                                                                                                                      >  1:227180/1‑148 (MQ=255)
                                gCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCaa                                                                                                                    <  8:119070/149‑1 (MQ=255)
                                 cAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAAt                                                                                                                   <  5:190565/149‑1 (MQ=255)
                                 cAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAAt                                                                                                                   <  7:79754/149‑1 (MQ=255)
                                 cAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAAt                                                                                                                   <  1:361061/149‑1 (MQ=255)
                                  aaCGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGTGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATg                                                                                                                  <  3:8898/149‑1 (MQ=255)
                                        gCCTTGACGCCGATGGCTGATGGGAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGTGGTGCGCGAGCAATTCCACGATGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAgcg                                                                                                            <  7:345054/149‑1 (MQ=255)
                                        gCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGATTTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAgcg                                                                                                            <  2:58868/149‑1 (MQ=255)
                                        gCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAgc                                                                                                             <  1:354799/148‑1 (MQ=255)
                                          cTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAgcgc                                                                                                           >  8:49506/1‑148 (MQ=255)
                                            tCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGaa                                                                                                        >  5:155910/1‑149 (MQ=255)
                                                 ccGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGaa                                                                                                   <  1:3117/149‑1 (MQ=255)
                                                 ccGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGaa                                                                                                   <  1:120847/149‑1 (MQ=255)
                                                           aTGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCt                                                                                          <  3:369294/148‑1 (MQ=255)
                                                                aaGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTgcgc                                                                                      >  6:56534/1‑147 (MQ=255)
                                                                    cgcTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCATTGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCt                                                                                 >  6:205606/1‑148 (MQ=255)
                                                                         gATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACg                                                                            >  8:296019/1‑148 (MQ=255)
                                                                              aGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGaaa                                                                      <  7:296019/149‑1 (MQ=255)
                                                                               gCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCTCACAACCAGCTGCGCCAGCTCAACGAAGAAAc                                                                     >  4:26347/1‑149 (MQ=255)
                                                                                   aaTTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCAt                                                                  <  6:26712/148‑1 (MQ=255)
                                                                                   aaTTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATc                                                                 <  2:227180/149‑1 (MQ=255)
                                                                                     tttGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATccc                                                               >  8:129771/1‑147 (MQ=255)
                                                                                     tttGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCAc                                                               >  1:377155/1‑149 (MQ=255)
                                                                                          ccAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTgg                                                          >  6:25477/1‑149 (MQ=255)
                                                                                                  cGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCACCGAAGAAACCATCACCCTGGCACTGCAg                                                  >  4:38077/1‑149 (MQ=255)
                                                                                                  cGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAg                                                  <  8:140245/149‑1 (MQ=255)
                                                                                                    gaggagGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGcc                                                >  5:13597/1‑149 (MQ=255)
                                                                                                    gaggagGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGcc                                                >  3:123598/1‑149 (MQ=255)
                                                                                                       gaggTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGtt                                             >  2:7120/1‑149 (MQ=255)
                                                                                                             cgcgAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGcc                                                <  1:263779/140‑1 (MQ=255)
                                                                                                             cgcgAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGcc                                                >  2:263779/1‑140 (MQ=255)
                                                                                                             cgcgAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTgg                                        <  5:25477/148‑1 (MQ=255)
                                                                                                              gcgAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGgcg                                      >  3:336046/1‑149 (MQ=255)
                                                                                                              gcgAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGgcg                                      >  1:386009/1‑149 (MQ=255)
                                                                                                               cgAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCg                                               >  7:129823/1‑139 (MQ=255)
                                                                                                               cgAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCg                                               <  8:129823/139‑1 (MQ=255)
                                                                                                                gAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGcgcc                                    >  2:94502/1‑149 (MQ=255)
                                                                                                                      ttCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTc                              >  5:181190/1‑149 (MQ=255)
                                                                                                                        ccACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCg                             <  4:336046/148‑1 (MQ=255)
                                                                                                                          aCGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGa                           >  2:208429/1‑148 (MQ=255)
                                                                                                                            gAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCCGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGaaaa                        <  7:227720/149‑1 (MQ=255)
                                                                                                                              gtaCATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGATTTCAAGCATGACGCAATTGAATCGCGAATCTAACAACCAGATGCGCCAGCTAAACGAAGAAACCATCAACCTTGCACTGCAGACGATGCTTGCGCCTCTGTCGGAAAAg                       >  1:307135/3‑148 (MQ=255)
                                                                                                                               gACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATCGAAGCGCGAATCGAACAACCAGCTGCGCGAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGAAAAGTa                     >  5:274699/1‑149 (MQ=255)
                                                                                                                                 cATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCAAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGAAAAGTACg                   >  6:263007/1‑149 (MQ=255)
                                                                                                                                     ctcttcgttGAGGCGCTCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGAAAAGTACGCCg                >  7:226377/8‑148 (MQ=255)
                                                                                                                                       ccTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGaa                          >  4:50111/1‑136 (MQ=255)
                                                                                                                                       ccTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGaa                          <  3:50111/136‑1 (MQ=255)
                                                                                                                                       ccTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGAAAAGTACGCCGaa              >  2:265795/1‑148 (MQ=255)
                                                                                                                                         tGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGAAAAGTACGCCGAAAAc           >  5:262091/1‑149 (MQ=255)
                                                                                                                                                  agCGCCTCAACCAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGGGAATCGAACAAACAGCTCCGCCAGCTCAACGAAGAACCCATCACACTCGCACTGCAGCCGTAGCTGGCGCCGCTGTCGGAAAAGTACGCCGAAAACGCGGCGgtg  >  8:95400/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGAAAAGGACGTGGCCTTGTACCGCGATGCCAGCAACGTCGCCTTCACGCCGATGGCTGATGGCAAGGCGCTGGATGAAGCTGAATTTGCCCAATTGCCGGAGGAGGTGCGCGAGCAATTCCACGAGGACATTGCCCTGCTTGAGGAGCGCCTCAACGAAGAGCTGGCCAGCCTGCCGCAATGGAAGCGCGAATCGAACAACCAGCTGCGCCAGCTCAACGAAGAAACCATCACCCTGGCACTGCAGCCGTTGCTGGCGCCGCTGTCGGAAAAGTACGCCGAAAACGCGGCGGTG  >  NC_002947/794365‑794659

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