Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 794,955 T→A 5.3% L366Q (CTG→CAG)  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,9550TA 5.3% 241.3 / 2.3 75L366Q (CTG→CAG) PP_0680ATP‑dependent protease‑like protein
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base A (2/2);  total (36/39)
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.87e-01

GGGGCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGC  >  NC_002947/794807‑795098
                                                                                                                                                    |                                                                                                                                               
ggggCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCt                                                                                                                                                 <  7:133156/149‑1 (MQ=255)
 gggCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGTGGCTTCCTGATTCTg                                                                                                                                                <  1:73605/149‑1 (MQ=255)
 gggCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTg                                                                                                                                                >  8:197046/1‑149 (MQ=255)
  ggCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCAGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCAACGTGCCAACGTCGGCTTCCTGATTCTgg                                                                                                                                               <  4:181956/149‑1 (MQ=255)
   gcgcCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTgg                                                                                                                                               >  6:325656/1‑148 (MQ=255)
    cgcCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGaa                                                                                                                                             <  4:141845/149‑1 (MQ=255)
               ttCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAAt                                                                                                                                  >  8:358287/1‑149 (MQ=255)
                   aGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                              <  7:95573/149‑1 (MQ=255)
                   aGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                              <  3:131525/149‑1 (MQ=255)
                    gCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTgg                                                                                                                             >  2:53858/1‑149 (MQ=255)
                     ccGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCACCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTggg                                                                                                                            >  8:251306/1‑149 (MQ=255)
                        cACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGa                                                                                                                         <  5:325656/149‑1 (MQ=255)
                         aCCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACAGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAg                                                                                                                        >  1:22484/1‑149 (MQ=255)
                              aCCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGtt                                                                                                                   >  8:106589/1‑149 (MQ=255)
                                 tATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCg                                                                                                                 >  3:67552/1‑148 (MQ=255)
                                  aTGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTa                                                                                                               >  3:60799/1‑149 (MQ=255)
                                   tGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTAt                                                                                                              >  1:374938/1‑149 (MQ=255)
                                   tGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCACCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGt                                                                                                                <  7:251306/147‑1 (MQ=255)
                                       aaTTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCt                                                                                                                               >  5:307936/1‑128 (MQ=255)
                                       aaTTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCt                                                                                                                               <  6:307936/128‑1 (MQ=255)
                                       aaTTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCAGAATCAGGAAGCCgccgtt                                                                                                                                    <  6:105419/123‑6 (MQ=255)
                                       aaTTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCAGAATCAGGAAGCCgccgtt                                                                                                                                    >  5:105419/1‑118 (MQ=255)
                                                 gCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCaa                                                                                                 <  4:67552/148‑1 (MQ=255)
                                                  ccGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAgc                                                                                               <  7:106589/149‑1 (MQ=255)
                                                        cGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCt                                                                                                                               >  3:75568/1‑111 (MQ=255)
                                                        cGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCt                                                                                                                               <  4:75568/111‑1 (MQ=255)
                                                         gAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGtt                                                                                        <  2:349299/149‑1 (MQ=255)
                                                         gAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGtt                                                                                        <  2:344877/149‑1 (MQ=255)
                                                          aaTACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTg                                                                                       >  7:101865/1‑149 (MQ=255)
                                                            tACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCa                                                                                     <  2:374938/149‑1 (MQ=255)
                                                                  aCCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCaa                                                                                                 <  1:361210/131‑1 (MQ=255)
                                                                  aCCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCaa                                                                                                 >  2:361210/1‑131 (MQ=255)
                                                                  aCCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTcgcg                                                                               >  6:383751/1‑149 (MQ=255)
                                                                      aCCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAg                                                                           <  3:270323/149‑1 (MQ=255)
                                                                        cAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCt                                                                         >  7:316476/1‑149 (MQ=255)
                                                                         aaGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCt                                                                         >  5:118028/1‑148 (MQ=255)
                                                                          aGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGa                                                                                                          >  5:11611/1‑114 (MQ=255)
                                                                          aGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGa                                                                                                          <  6:11611/114‑1 (MQ=255)
                                                                           ggTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGaa                                                                      <  4:60799/149‑1 (MQ=255)
                                                                            gTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTc                                                                                  <  5:114344/136‑1 (MQ=255)
                                                                            gTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTc                                                                                  >  6:114344/1‑136 (MQ=255)
                                                                            gTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAg                                                                     <  8:316476/149‑1 (MQ=255)
                                                                               cgcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATg                                                                  <  8:101865/149‑1 (MQ=255)
                                                                                 cTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGa                                                                <  1:53858/149‑1 (MQ=255)
                                                                                    tataCCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATg                                                                  >  7:390452/1‑144 (MQ=255)
                                                                                    tataCCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATg                                                                  <  8:390452/144‑1 (MQ=255)
                                                                                      taCCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGc                                                           >  6:125179/1‑149 (MQ=255)
                                                                                      taCCTCCTATCGCCAACTGCGCCCAGGCGAGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGc                                                           <  8:200074/149‑1 (MQ=255)
                                                                                         ctcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTc                                                                                                                  >  3:257528/1‑91 (MQ=255)
                                                                                         ctcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTc                                                                                                                  <  4:257528/91‑1 (MQ=255)
                                                                                         ctcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGTAGTCGc                                                           <  7:261964/146‑1 (MQ=255)
                                                                                         ctcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGc                                                           >  8:261964/1‑146 (MQ=255)
                                                                                         ctcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGa                                                        >  5:382137/1‑149 (MQ=255)
                                                                                             tATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCgg                                                    >  7:98234/1‑149 (MQ=255)
                                                                                             tATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCgg                                                    >  7:57681/1‑149 (MQ=255)
                                                                                                 gCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGa                                                 <  5:125179/148‑1 (MQ=255)
                                                                                                  ccAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGaa                                                                                                                                     <  8:282334/63‑1 (MQ=255)
                                                                                                  ccAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGaa                                                                                                                                     >  7:282334/1‑63 (MQ=255)
                                                                                                            ccAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCgtgt                                     >  8:317134/1‑149 (MQ=255)
                                                                                                            ccAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCCTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCgtgt                                     >  1:334762/1‑149 (MQ=255)
                                                                                                             cAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGtt                                    <  7:197046/149‑1 (MQ=255)
                                                                                                                 cgcgcTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCa                                <  5:176556/149‑1 (MQ=255)
                                                                                                                    gcTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGa                                                        <  8:137818/122‑1 (MQ=255)
                                                                                                                    gcTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGa                                                        >  7:137818/1‑122 (MQ=255)
                                                                                                                    gcTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTg                             >  7:252779/1‑149 (MQ=255)
                                                                                                                     cTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGt                            <  8:252779/149‑1 (MQ=255)
                                                                                                                           cGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGAGTATCCACTGTAATCCt                      >  3:219215/1‑149 (MQ=255)
                                                                                                                            gTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTg                     <  6:118028/149‑1 (MQ=255)
                                                                                                                               ccAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCa                   <  6:382137/148‑1 (MQ=255)
                                                                                                                                      cggcTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAgg           >  3:304862/1‑149 (MQ=255)
                                                                                                                                      cggcTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAgg           <  2:334762/149‑1 (MQ=255)
                                                                                                                                      cggcTTCCAGAATCAGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAATCTGAAGAATGAGTAGCCCATAGGTTTGCTGGGCCGTGTTTGCACTg                             >  2:397004/1‑131 (MQ=255)
                                                                                                                                      cggcTTCCAGAATCAGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTg                             <  1:397004/131‑1 (MQ=255)
                                                                                                                                             cTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAg                                                <  5:45519/105‑1 (MQ=255)
                                                                                                                                             cTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAg                                                >  6:45519/1‑105 (MQ=255)
                                                                                                                                                aTTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGc  <  7:317134/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGGGCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGC  >  NC_002947/794807‑795098

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