Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 794,955 T→A 6.3% I365N (ATT→AAT)  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 6.3% 192.9 / 4.2 63I365N (ATT→AAT) PP_0680ATP‑dependent protease‑like protein
Reads supporting (aligned to +/- strand):  ref base T (25/34);  new base A (2/2);  total (27/36)
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.37e-01

GGGGCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGC  >  NC_002947/794810‑795101
                                                                                                                                                 |                                                                                                                                                  
ggggCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTg                                                                                                                                                <  1:310989/150‑1 (MQ=255)
ggggCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTg                                                                                                                                                <  1:586338/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                                                                                                                              <  2:1164557/150‑1 (MQ=255)
                  gAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTg                                                                                                                              <  1:893052/150‑1 (MQ=255)
                       gCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGa                                                                                                                         <  2:1075938/150‑1 (MQ=255)
                         aCCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCAGAATCAGGAAGcc                                                                                                                                          >  1:837736/1‑131 (MQ=255)
                         aCCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCAGAATCAGGAAGcc                                                                                                                                          <  2:837736/131‑1 (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:1041107/1‑150 (MQ=255)
                                                                                          tcctATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATc                                                      >  1:294197/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                                                                       >  1:461852/1‑109 (MQ=255)
                                                                                                                  gcgcTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGa                                                                       <  2:461852/109‑1 (MQ=255)
                                                                                                                     cTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTc                           >  2:339461/1‑150 (MQ=255)
                                                                                                                          ccGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCt                      <  2:1189917/150‑1 (MQ=255)
                                                                                                                              gCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCgtg                                      <  1:1195068/130‑1 (MQ=255)
                                                                                                                              gCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCgtg                                      >  2:1195068/1‑130 (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)
                                                                                                                                                 |                                                                                                                                                  
GGGGCGCCGGTGGTGTTCGAGCCGCACCCGACCTATGACAATTTGTTTGGCCGTATCGAATACAGCACCGACCAAGGTGCGCTGTATACCTCCTATCGCCAACTGCGCCCAGGCGCGCTGCACCGTGCCAACGGCGGCTTCCTGATTCTGGAAGCCGAGAAAATGCTGGGTGAGCCGTTCGTATGGGACGCGCTCAAGCGTGCGTTGCAGTCGCGCAAGCTGAAGATGGAGTCGCCGATCGGTGAGCTGGGCCGTGTTGCCACTGTCAGCCTGCAGCCGCAGGTGATTCCGC  >  NC_002947/794810‑795101

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