Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,002,280 G→C 9.2% A26G (GCC→GGC)  amaA ← L‑pipecolate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,002,2800GC 9.2% 225.9 / 13.0 76A26G (GCC→GGC) amaAL‑pipecolate oxidase
Reads supporting (aligned to +/- strand):  ref base G (39/30);  new base C (3/4);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

ACGTTGCGCCCGGAGCCACCGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCT  >  NC_002947/6002133‑6002429
                                                                                                                                                   |                                                                                                                                                     
aCGTTGCGCCCGGAGCCACCGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGcc                                                                                                                                                     <  1:288357/150‑1 (MQ=255)
          cGGAGCCACCGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGGCCGGTGa                                                                                                                                           <  2:1254089/150‑1 (MQ=255)
                 aCCGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTAGCGCTCACCGGCCAGGGCCTGGGCGGCGAcggtcgg                                                                                                                                    >  2:824234/1‑150 (MQ=255)
                   cGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGtt                                                                                                                                  >  2:1067982/1‑150 (MQ=255)
                   cGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGtt                                                                                                                                  <  2:277243/150‑1 (MQ=255)
                         cGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGa                                                                                                                            <  1:1067982/150‑1 (MQ=255)
                          gATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTTGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGAc                                                                                                                           <  1:819395/150‑1 (MQ=255)
                          gATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGAc                                                                                                                           >  1:919011/1‑150 (MQ=255)
                          gATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGAc                                                                                                                           <  2:1291415/150‑1 (MQ=255)
                               tCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCCGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGtt                                                                                                                      <  1:824234/150‑1 (MQ=255)
                                  aGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTccc                                                                                                                   >  2:168387/1‑150 (MQ=255)
                                  aGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTccc                                                                                                                   >  2:828719/1‑150 (MQ=255)
                                       tCCAGCACAATCACCGCCCTGCACTGTTCGAGCAGGTGGATGGCCGCCGACCGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTAACCGGGCAGGGCCTGGGCGGCGACGGCCGATGTGCTGTGGGGTTACCCCAAg                                                                                                              <  2:558011/150‑1 (MQ=255)
                                         cAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCa                                                                                                            >  2:686538/1‑150 (MQ=255)
                                         cAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCa                                                                                                            >  1:770357/1‑150 (MQ=255)
                                         cAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCa                                                                                                            >  1:682234/1‑150 (MQ=255)
                                           gCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCAtt                                                                                                          <  1:1214067/150‑1 (MQ=255)
                                           gCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCAtt                                                                                                          >  2:812487/1‑150 (MQ=255)
                                           gCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCAtt                                                                                                          >  1:612774/1‑150 (MQ=255)
                                             aCGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCAttgt                                                                                                        <  2:612774/150‑1 (MQ=255)
                                               gATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGgccggtg                                                                                                                                            >  2:280371/1‑110 (MQ=255)
                                               gATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGgccggtg                                                                                                                                            <  1:280371/112‑3 (MQ=255)
                                                  cACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACg                                                                                                   <  1:1197390/150‑1 (MQ=255)
                                                         ttGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCg                                                                                            >  1:922704/1‑150 (MQ=255)
                                                           gCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGa                                                                                          <  1:813705/150‑1 (MQ=255)
                                                            cccTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGAc                                                                                         >  1:1172483/1‑150 (MQ=255)
                                                             ccTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACa                                                                                        >  2:1367824/1‑150 (MQ=255)
                                                             ccTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACa                                                                                        >  2:700338/1‑150 (MQ=255)
                                                             ccTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGgccggt                                                                                                                                             <  1:517241/97‑2 (MQ=255)
                                                             ccTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGgccggt                                                                                                                                             >  2:517241/1‑96 (MQ=255)
                                                             ccTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGGCCGGTGa                                                                                                                                           >  2:87264/1‑99 (MQ=255)
                                                             ccTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGCCCAGGCCCTGGCCGGTGa                                                                                                                                           <  1:87264/99‑1 (MQ=255)
                                                             ccTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGATCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCGCACAAGCATTGTTGACGCAGTTCGGACa                                                                                        >  1:454641/1‑150 (MQ=255)
                                                                  tCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGcc                                                                                   <  2:180998/150‑1 (MQ=255)
                                                                     aGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACAAACTTCGGACATGGCCGGc                                                                                <  1:825405/150‑1 (MQ=255)
                                                                        aGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTcc                                                                             >  2:284267/1‑150 (MQ=255)
                                                                             gATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCaa                                                                        <  1:812487/150‑1 (MQ=255)
                                                                                tgccgccGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAAttt                                                                     >  2:94063/2‑150 (MQ=255)
                                                                                   cgccgACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGtt                                                                  >  2:597057/1‑150 (MQ=255)
                                                                                            gCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTggg                                                         >  1:893060/1‑150 (MQ=255)
                                                                                             ccGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACgtgt                                                                                                                       >  2:597166/1‑87 (MQ=255)
                                                                                             ccGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACgtgt                                                                                                                       <  1:597166/87‑1 (MQ=255)
                                                                                             ccGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTgggg                                                        >  1:350076/1‑150 (MQ=255)
                                                                                              cGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGc                                                                                    <  1:178181/121‑1 (MQ=255)
                                                                                              cGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGc                                                                                    >  2:178181/1‑121 (MQ=255)
                                                                                              cGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGc                                                       >  1:18468/1‑150 (MQ=255)
                                                                                                 tAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATAGCCGGCGCCAGCAATTTGTTCTTGTTGGGGCCGc                                                    >  2:824081/1‑150 (MQ=255)
                                                                                                    tACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCAttgt                                                                                                        <  1:1104543/95‑1 (MQ=255)
                                                                                                    tACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCAttgt                                                                                                        >  2:1104543/1‑95 (MQ=255)
                                                                                                    tACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTCTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAgg                                                 <  2:922704/150‑1 (MQ=255)
                                                                                                              gATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCttt                                       <  2:855246/150‑1 (MQ=255)
                                                                                                              gATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTCCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCttt                                       >  2:966727/1‑150 (MQ=255)
                                                                                                               aTCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTc                                      <  2:18468/150‑1 (MQ=255)
                                                                                                                   cacaGACATCGGCTTTGTGCTCACCGGCCAGGGCCTAGGCGGAGACGGTCGAACTGGCGGCACACCCTGACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCgg                                  <  1:168387/150‑1 (MQ=255)
                                                                                                                   cacaGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCgg                                  >  2:131850/1‑150 (MQ=255)
                                                                                                                   cacaGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCgg                                  <  2:753912/150‑1 (MQ=255)
                                                                                                                   cacaGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCgg                                  <  1:131850/150‑1 (MQ=255)
                                                                                                                   cacaGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCgg                                  <  2:1330319/150‑1 (MQ=255)
                                                                                                                    acaGACATCGGCCTTGTGCTCACCGTCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCggg                                 >  1:1321805/1‑150 (MQ=255)
                                                                                                                    acaGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGAc                                                                                                    >  2:243804/1‑83 (MQ=255)
                                                                                                                    acaGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGAc                                                                                                    <  1:243804/83‑1 (MQ=255)
                                                                                                                    acaGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCggg                                 >  1:390555/1‑150 (MQ=255)
                                                                                                                      aGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTa                               >  1:453286/1‑150 (MQ=255)
                                                                                                                          aTCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAAcc                           >  1:862285/1‑150 (MQ=255)
                                                                                                                             ggCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGc                        <  1:1199521/150‑1 (MQ=255)
                                                                                                                             ggCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGc                        <  2:1172483/150‑1 (MQ=255)
                                                                                                                               ccTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTc                      >  1:16918/1‑150 (MQ=255)
                                                                                                                               ccTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTc                      >  2:718398/1‑150 (MQ=255)
                                                                                                                               ccTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTc                      >  2:1120529/1‑150 (MQ=255)
                                                                                                                                      cTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAgg               <  2:362149/150‑1 (MQ=255)
                                                                                                                                      cTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAgg               <  1:1037952/150‑1 (MQ=255)
                                                                                                                                       tCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAggg              >  1:913043/1‑150 (MQ=255)
                                                                                                                                         aCCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACCGGGGt            >  2:1249853/1‑150 (MQ=255)
                                                                                                                                            ggCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCgg         <  1:1120529/150‑1 (MQ=255)
                                                                                                                                             gCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGgt        >  2:488591/1‑150 (MQ=255)
                                                                                                                                               cAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGgtgt      >  1:400098/1‑150 (MQ=255)
                                                                                                                                                   gCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCt  >  1:408969/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACGTTGCGCCCGGAGCCACCGTGGCCGATCTTCCAGGCCTCCAGCACGATCACCGACTTGCCCTGTTCGAGCAGGTGGATGGCCGCCGACAGGCCGGTAATACCGCCGCCGATCACACAGACATCGGCCTTGTGCTCACCGGCCAGGGCCTGGGCGGCGACGGTCGGTTTGCTGACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCT  >  NC_002947/6002133‑6002429

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: