Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 794,952 T→A 6.6% L364Q (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,9520TA 6.6% 189.8 / 4.2 61L364Q (CTG→CAG) PP_0680ATP‑dependent protease‑like protein
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base A (2/2);  total (27/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 = 5.96e-01

TGACGGTGGGGCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGC  >  NC_002947/794803‑795101
                                                                                                                                                     |                                                                                                                                                     
tGACGGTGGGGCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCt                                                                                                                                                       >  1:1092512/1‑150 (MQ=255)
tGACGGTGGGGCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCt                                                                                                                                                       >  2:965212/1‑150 (MQ=255)
       ggggCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTg                                                                                                                                                <  1:586338/150‑1 (MQ=255)
       ggggCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTg                                                                                                                                                <  1:310989/150‑1 (MQ=255)
          gcgcCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCCCCGTGCCAACGGCGGCTTCCTGATTCTGGaa                                                                                                                                             >  2:9269/1‑150 (MQ=255)
          gcgcCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGaa                                                                                                                                             <  2:1104720/150‑1 (MQ=255)
          gcgcCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGGATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGaa                                                                                                                                             <  1:910024/150‑1 (MQ=255)
               ggtggtGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCga                                                                                                                                        <  1:126261/150‑1 (MQ=255)
                    tgtTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCAAACGGCGGCTTCCTGATTCTGGAAGCCGAGaaaa                                                                                                                                   >  1:311711/1‑150 (MQ=255)
                      ttCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGcc                                                                                                                                          <  2:973048/141‑1 (MQ=255)
                      ttCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGcc                                                                                                                                          >  1:973048/1‑141 (MQ=255)
                         gAGCTGCACTCTACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                              <  2:311711/150‑1 (MQ=255)
                         gAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                              <  1:893052/150‑1 (MQ=255)
                         gAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                              <  2:1164557/150‑1 (MQ=255)
                              gCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGa                                                                                                                         <  2:1075938/150‑1 (MQ=255)
                                aCCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCAGAATCAGGAAGcc                                                                                                                                          <  2:837736/131‑1 (MQ=255)
                                aCCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCAGAATCAGGAAGcc                                                                                                                                          >  1:837736/1‑131 (MQ=255)
                                   cGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGt                                                                                                                    >  2:1086559/1‑150 (MQ=255)
                                         aTGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTAt                                                                                                              >  2:1057696/1‑150 (MQ=255)
                                          tGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTc                                                                                                                  >  1:360315/1‑145 (MQ=255)
                                          tGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTc                                                                                                                  <  2:360315/145‑1 (MQ=255)
                                          tGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATg                                                                                                             <  2:695141/150‑1 (MQ=255)
                                                tttgtttgGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGAcgc                                                                                                       >  2:511932/1‑150 (MQ=255)
                                                  tgtttgGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACgcgc                                                                                                     <  2:815045/150‑1 (MQ=255)
                                                        gCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAgc                                                                                               >  2:749719/1‑150 (MQ=255)
                                                                 aaTACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGCAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGc                                                                                      <  1:689272/150‑1 (MQ=255)
                                                                 aaTACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAgcgt                                                                                             <  1:787303/143‑1 (MQ=255)
                                                                 aaTACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAgcgt                                                                                             >  2:787303/1‑143 (MQ=255)
                                                                  aTACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCa                                                                                     <  2:92370/150‑1 (MQ=255)
                                                                        cACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCga                                                                                                                                        >  1:873177/1‑93 (MQ=255)
                                                                        cACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCga                                                                                                                                        <  2:873177/93‑1 (MQ=255)
                                                                         aCNGACCAAGGTGCGCTGTATACCTCCNATCGCCAACTGCNCCCAGGNGNGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCgcgc                                                                              >  1:692439/1‑150 (MQ=255)
                                                                                 aGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGaa                                                                      >  1:1019998/1‑150 (MQ=255)
                                                                                 aGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGAACTGTGCTGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGaa                                                                      <  1:318095/150‑1 (MQ=255)
                                                                                  ggTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAg                                                                     <  2:1019998/150‑1 (MQ=255)
                                                                                   cggtgcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGa                                                                    <  2:520269/146‑1 (MQ=255)
                                                                                    tGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGAt                                                                   <  2:988121/150‑1 (MQ=255)
                                                                                    tGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCAGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGAt                                                                   <  1:511932/150‑1 (MQ=255)
                                                                                     gcgcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTAAACGCCGAGAAAATGCTGTCGGTGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATg                                                                  <  1:749719/150‑1 (MQ=255)
                                                                                     gcgcTGCATACCTCCTATAGACCGATGAGCCTGAGTGTGCCGGACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATg                                                                  <  2:918303/150‑1 (MQ=255)
                                                                                       gcTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGa                                                                >  1:251910/1‑150 (MQ=255)
                                                                                                ctcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGAt                                                       >  1:1189917/1‑150 (MQ=255)
                                                                                                ctcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGAt                                                       >  2:507035/1‑150 (MQ=255)
                                                                                                 tcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATc                                                      >  1:294197/1‑150 (MQ=255)
                                                                                                 tcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATc                                                      >  1:1041107/1‑150 (MQ=255)
                                                                                                    tATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGt                                                   <  2:905434/150‑1 (MQ=255)
                                                                                                         ccAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCt                                              >  2:790456/1‑150 (MQ=255)
                                                                                                             cTGCGCCCAGGCGCGCTGGACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTgggn                                          <  2:692439/150‑2 (MQ=255)
                                                                                                                 gcCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCgtg                                      <  1:507035/150‑1 (MQ=255)
                                                                                                                   ccAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGtt                                    >  2:502221/1‑150 (MQ=255)
                                                                                                                         gcgcTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGa                                                                       <  2:461852/109‑1 (MQ=255)
                                                                                                                         gcgcTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGa                                                                       >  1:461852/1‑109 (MQ=255)
                                                                                                                            cTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTc                           >  2:339461/1‑150 (MQ=255)
                                                                                                                                 ccGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCt                      <  2:1189917/150‑1 (MQ=255)
                                                                                                                                     gCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCgtg                                      >  2:1195068/1‑130 (MQ=255)
                                                                                                                                     gCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCgtg                                      <  1:1195068/130‑1 (MQ=255)
                                                                                                                                       ttttctcggcTTCCAGAATCAGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGcc                >  1:728917/7‑150 (MQ=255)
                                                                                                                                       cAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGcc                >  2:242067/1‑150 (MQ=255)
                                                                                                                                         aCGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTgg                                            >  1:128955/1‑120 (MQ=255)
                                                                                                                                         aCGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTgg                                            <  2:128955/120‑1 (MQ=255)
                                                                                                                                           ggcggcTTCCTGATTCGGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAg            <  2:381810/150‑1 (MQ=255)
                                                                                                                                             cggcTTCCAGAATCAGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGt          <  2:728917/150‑1 (MQ=255)
                                                                                                                                                    cTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTGCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCg   >  2:419091/1‑150 (MQ=255)
                                                                                                                                                     tGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGc  <  1:790456/150‑1 (MQ=255)
                                                                                                                                                     tGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGc  <  1:1086559/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACGGTGGGGCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGC  >  NC_002947/794803‑795101

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