Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,658,771 G→T 5.4% G259C (GGC→TGC)  sahR → methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,658,7710GT 5.4% 361.6 / 5.4 111G259C (GGC→TGC) sahRmethionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (53/52);  new base T (4/2);  total (57/54)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 113 total aligned reads displayed.

AAGGGCTCGACAACGTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGCCTGTACGT  >  NC_002947/5658622‑5658916
                                                                                                                                                     |                                                                                                                                                 
 aGGGCTCGACAACGTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGTCAGCTGGCCAAACGGGTGAAGGCAggcgg                                                                                                                                                  >  1:598173/1‑150 (MQ=255)
     cTCGACAACGTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGc                                                                                                                                              <  1:2226725/150‑1 (MQ=255)
      tCGACAACGTGGACCTGCAGTTGGCCGATGCACTGGGGGCAACGGATGTGGCCGCCGACTGCGGCGTGATGAACATAGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGAGTGAAGGCAGGCGGCAGcc                                                                                                                                             >  2:775625/1‑150 (MQ=255)
      tCGACAACGTGAACCTGCAGTTGGCCGATGCACTGGGAGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGAAGCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGcc                                                                                                                                             >  1:1955365/1‑150 (MQ=255)
         aCAACGTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCctgc                                                                                                                                          <  2:315490/150‑1 (MQ=255)
          cAACGTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCctgct                                                                                                                                         <  2:623733/150‑1 (MQ=255)
           aaCGTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCctgctg                                                                                                                                        <  1:514004/150‑1 (MQ=255)
           aaCGTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCctgctg                                                                                                                                        <  2:819117/150‑1 (MQ=255)
              gTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCa                                                                                                                                                >  1:1091756/1‑139 (MQ=255)
              gTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCa                                                                                                                                                <  2:1091756/139‑1 (MQ=255)
              gTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTc                                                                                                                                     >  2:796345/1‑150 (MQ=255)
              gTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTc                                                                                                                                     >  2:1185063/1‑150 (MQ=255)
                  aCCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCg                                                                                                                                 >  1:2187691/1‑150 (MQ=255)
                   ccTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGAACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAACCTGCTGGTCACCGa                                                                                                                                >  1:6903/1‑150 (MQ=255)
                          ttGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGc                                                                                                                         <  1:1967765/150‑1 (MQ=255)
                           tGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCa                                                                                                                        <  2:240216/150‑1 (MQ=255)
                             gCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGc                                                                                                                      >  1:1720646/1‑150 (MQ=255)
                               cGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCa                                                                                                                    >  2:2296642/1‑150 (MQ=255)
                                     aCTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCctgctg                                                                                                                                        <  2:814802/124‑1 (MQ=255)
                                     aCTGGGTGCAACGGATGTGGCCGACGAGTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGGCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGAGGCAGCGTGCTGGTCAGCGAACTGTGCAGCCATGACCa                                                                                                              >  1:194193/1‑150 (MQ=255)
                                       tGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAgg                                                                                                            >  1:53973/1‑150 (MQ=255)
                                        gggTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAggg                                                                                                           >  1:2114646/1‑150 (MQ=255)
                                         ggTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAgggg                                                                                                          >  1:1805867/1‑150 (MQ=255)
                                         ggTGCAACGGATGTGGCCGCCGACTCAGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAgggg                                                                                                          <  1:1811457/150‑1 (MQ=255)
                                           tGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTg                                                                                                        >  2:1503397/1‑150 (MQ=255)
                                            gCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTgg                                                                                                       <  1:188124/150‑1 (MQ=255)
                                               aCGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCa                                                                                                                        <  2:405997/130‑1 (MQ=255)
                                               aCGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCa                                                                                                                        >  1:405997/1‑130 (MQ=255)
                                                 ggATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGgcgcg                                                                                                  >  2:1805249/1‑150 (MQ=255)
                                                     gtgGCCGCCGACTGCGTTGTGCTGAACGTGGTGATTGACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGaa                                                                                              <  2:194193/150‑1 (MQ=255)
                                                     gtgGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGaa                                                                                              <  1:1816500/150‑1 (MQ=255)
                                                      tgGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAg                                                                                             <  1:1640356/150‑1 (MQ=255)
                                                       ggCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGc                                                                                            >  1:1019702/1‑150 (MQ=255)
                                                          cgccgACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTg                                                                                         <  1:1185063/150‑1 (MQ=255)
                                                           gccgACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTgc                                                                                        >  1:482422/1‑150 (MQ=255)
                                                              gACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTgcggc                                                                                     <  1:113767/150‑1 (MQ=255)
                                                              gACTGCGTTGTGCTGAACATGGGGCTGGACCATGTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACTTAGCTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCATgcggt                                                                                     <  2:407216/150‑2 (MQ=255)
                                                                ttGCGTTGTGCTGAACGTGATGCCGCACCCTGTCATCTACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGa                                                                                   <  2:1720646/149‑1 (MQ=255)
                                                                 tGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGAc                                                                                  >  1:1026643/1‑150 (MQ=255)
                                                                 tGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGAc                                                                                  >  1:116337/1‑150 (MQ=255)
                                                                 tGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCCGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGAc                                                                                  <  1:2380405/150‑1 (MQ=255)
                                                                    gTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACctc                                                                               >  2:20146/1‑150 (MQ=255)
                                                                    gTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACctc                                                                               >  1:868512/1‑150 (MQ=255)
                                                                     ttGTGCTGAACATGGTGCTGCACCATTTCAGCGAGACGGCCCTGGCCTTGCGCCAGCGGGCCAAACGGGTGAAGGCAGGGGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACctct                                                                              >  2:1480903/1‑150 (MQ=255)
                                                                     ttGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACctct                                                                              >  1:1760882/1‑150 (MQ=255)
                                                                       gtgCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTctgg                                                                            <  2:1275671/150‑1 (MQ=255)
                                                                          cTGAACATGGTGCTGCACCATTTCAGTGACCCGGCCCTGGCCTTGCTCCAGCTGGCCAAACGTGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACAAGGGGTGGGCGCGCGAAGCCTGCGGCGCCCTctggctg                                                                         <  2:1955365/150‑1 (MQ=255)
                                                                          cTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTctggctg                                                                         >  1:984558/1‑150 (MQ=255)
                                                                          cTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTctggctg                                                                         >  1:1323060/1‑150 (MQ=255)
                                                                          cTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTctggctg                                                                         <  1:20146/150‑1 (MQ=255)
                                                                           tGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTctggctgg                                                                        <  1:399090/150‑1 (MQ=255)
                                                                           tGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTctggctgg                                                                        <  1:2296642/150‑1 (MQ=255)
                                                                           tGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTctggctgg                                                                        <  1:897595/150‑1 (MQ=255)
                                                                             aaCATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGc                                                                                                                                              >  1:707302/1‑78 (MQ=255)
                                                                             aaCATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGc                                                                                                                                              <  2:707302/78‑1 (MQ=255)
                                                                             aaCATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGc                                                                      <  2:984558/150‑1 (MQ=255)
                                                                             aaCATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAACCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGc                                                                      >  2:1324051/1‑150 (MQ=255)
                                                                                     gctgcACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACa                                                              <  1:1600367/150‑1 (MQ=255)
                                                                                      ctgcACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAg                                                             <  2:721939/150‑1 (MQ=255)
                                                                                       tgcACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAgg                                                            >  2:1301172/1‑150 (MQ=255)
                                                                                        gcACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGga                                                           >  2:936883/1‑150 (MQ=255)
                                                                                          aCCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCct                                                                                                                                            <  2:447949/67‑1 (MQ=255)
                                                                                          aCCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCct                                                                                                                                            >  1:447949/1‑67 (MQ=255)
                                                                                          aCCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTTACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGgacg                                                         >  2:860565/1‑150 (MQ=255)
                                                                                           ccATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGgacga                                                        >  1:1816364/1‑150 (MQ=255)
                                                                                           ccATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGgacga                                                        >  1:372254/1‑150 (MQ=255)
                                                                                            cATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGgacgac                                                       <  2:1816364/150‑1 (MQ=255)
                                                                                            cATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGgacgac                                                       >  2:2237364/1‑150 (MQ=255)
                                                                                             aTTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGAcc                                                      >  2:2401291/1‑150 (MQ=255)
                                                                                               ttCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTg                                                    <  2:1323060/150‑1 (MQ=255)
                                                                                               ttCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTg                                                    >  1:2271159/1‑150 (MQ=255)
                                                                                                  aGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGgcg                                                 <  1:2086333/150‑1 (MQ=255)
                                                                                                   gcgcCCCGGCCCTTGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGgcgc                                                >  2:1901561/1‑150 (MQ=255)
                                                                                                              cTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAAc                                     <  2:2187691/150‑1 (MQ=255)
                                                                                                              cTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAAc                                     <  1:796345/150‑1 (MQ=255)
                                                                                                                   cTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGga                                                           <  1:89952/123‑1 (MQ=255)
                                                                                                                   cTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGga                                                           >  2:89952/1‑123 (MQ=255)
                                                                                                                   cTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAAc                                     <  1:500244/145‑1 (MQ=255)
                                                                                                                   cTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAAc                                     >  2:500244/1‑145 (MQ=255)
                                                                                                                   cTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGc                                >  2:1685683/1‑150 (MQ=255)
                                                                                                                        gcCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGATCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGAcc                                                      >  2:414741/1‑123 (MQ=255)
                                                                                                                        gcCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGAcc                                                      <  1:414741/123‑1 (MQ=255)
                                                                                                                            gCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGAGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGc                       <  2:598173/150‑1 (MQ=255)
                                                                                                                              tGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAg                                                             >  1:987403/1‑110 (MQ=255)
                                                                                                                                  cAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGGCGGGCTGGCCCATGc                 >  2:490739/1‑150 (MQ=255)
                                                                                                                                  cAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGc                 >  1:2189847/1‑150 (MQ=255)
                                                                                                                                    aaCGGGTGAAGGCAGGCGGCAGCCTGCTAGTCACCGAACAGGTTACCCTGTGCCGGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCgg               >  1:2108613/1‑150 (MQ=255)
                                                                                                                                     aCGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGa              <  2:461003/150‑1 (MQ=255)
                                                                                                                                          tGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGcc         <  2:2398402/150‑1 (MQ=255)
                                                                                                                                           gAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGCCt        <  2:1019702/150‑1 (MQ=255)
                                                                                                                                             cagcaggcTGCCGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGga                                                           >  2:71716/3‑97 (MQ=255)
                                                                                                                                             cagcaggcTGCCGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGga                                                           <  1:71716/95‑1 (MQ=255)
                                                                                                                                             cagcaggcTGCCGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGCCTGt      >  2:376266/3‑150 (MQ=255)
                                                                                                                                             cagcaggcTGCCGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCCGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGcc         >  2:628534/3‑147 (MQ=255)
                                                                                                                                             cagcaggcTGCCGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCCGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGcc         <  1:628534/145‑1 (MQ=255)
                                                                                                                                             cagcaggcTGCCGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGGGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGCCTGt      >  1:1535424/3‑150 (MQ=255)
                                                                                                                                              ggcaggcGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGCCTGTa     >  2:1372709/1‑150 (MQ=255)
                                                                                                                                                 aggcGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGCCTGTACGt  <  2:868512/150‑1 (MQ=255)
                                                                                                                                                 aggcGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGCCTGTACGt  <  1:814133/150‑1 (MQ=255)
                                                                                                                                                 aggcGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCCGCCGGGCTGGCCCATGCGGACAGCCTGTACGt  <  1:2401291/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AAGGGCTCGACAACGTGAACCTGCAGTTGGCCGATGCACTGGGTGCAACGGATGTGGCCGCCGACTGCGTTGTGCTGAACATGGTGCTGCACCATTTCAGCGACCCGGCCCTGGCCTTGCGCCAGCTGGCCAAACGGGTGAAGGCAGGCGGCAGCCTGCTGGTCACCGAACTGTGCAGCCATGACCAGGGGTGGGCGCGCGAAGCCTGCGGCGACCTCTGGCTGGGCTTCGAACAGGACGACCTGGCGCGCTGGGCCAACGCTGCCGGGCTGGCCCATGCGGACAGCCTGTACGT  >  NC_002947/5658622‑5658916

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