Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,823,802 A→G intergenic (‑33/+47) recA ← / ← pncC DNA recombination/repair protein RecA/NMN aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,823,8020AG98.1% 185.3 / NA 53intergenic (‑33/+47)recA/pncCDNA recombination/repair protein RecA/NMN aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/28);  minor base C (1/0);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATCCGGCCCGGCGATACCACT  >  NC_000913/2823557‑2824040
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
ggAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCaa                                                                                                                                                                                                                                           <  2:357068/251‑1 (MQ=255)
 gAAGCTTCCGGTCCGTAGATTGCGACGATACGGCCCAGCGGCAGACCACCTGCCCCAAGCGCGATATCCATTGAAAGCGAACCGGTAGAGATTGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCGTTCTGTTTGGTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAat                                                                                                                                                                                                                                          <  2:232588/251‑1 (MQ=255)
     gttCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGt                                                                                                                                                                                                                                      <  2:267126/250‑1 (MQ=255)
           gTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGCTATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCCCATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGtt                                                                                                                                                                                                                                >  1:477251/1‑251 (MQ=255)
             ccGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGTAGCCTTTACCAAATTGTTTCTCAATCTTTCCCAGTTCTGCCGCCAACGCTTGCTGGTTGTTTTCGTCGATAGCCAGGTGTACTCCTGGCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTg                                                                                                                                                                                                                               <  1:431389/250‑1 (MQ=255)
               gTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAg                                                                                                                                                                                                                            <  2:57474/251‑1 (MQ=255)
               gTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAg                                                                                                                                                                                                                            <  2:401906/251‑1 (MQ=255)
                tAGATTTCGACTATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTGCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAAGGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCGGTCATGCCGGGTAATACCGGATGCTCAATATGTTCTGTTGAAg                                                                                                                                                                                                                            <  2:437647/250‑1 (MQ=255)
                           gaTACGGCACATCGGCCGACCACCTGCCCCAAGCGCGATATCCAGTGAACGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTACTACCGGCTGGTCAATATGTTCTGTTGAAGCAATTATACttt                                                                                                                                                                                                                >  2:388182/1‑249 (MQ=255)
                            aTACGGCCCATCGGCAGACCGCCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTGTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGTTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTa                                                                                                                                                                                                               <  1:189444/251‑1 (MQ=255)
                            aTACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGCTATCCAGTGACAGCGAACCGGTAGAGATGGTTTCCCCATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGt                                                                                                                                                                                                                >  1:29313/1‑250 (MQ=255)
                                     aTCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGCCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTACTACCGGATGGTCAATCTGTTCTGTTGAAGCAATTATACTGTATGCTCATAc                                                                                                                                                                                                      >  2:2923/1‑251 (MQ=255)
                                     aTCGGAAGAACACCTGCACCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGCTGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATAc                                                                                                                                                                                                      >  1:357395/1‑251 (MQ=255)
                                      ccGGCAGACCACCTGCGCCAAGCGCGATATCCCGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTGCACCCAGGCGCATGATGGAGCCTTTACGAAATTCTTTCTCAATCTTGCCCACTGCCTCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCTATTTTTACTCCTGTCAGGCCGGGCGTTACCGGATGTGCAGTATGGGCTGGTGAAGCGATGATACTTTATGCTCATACa                                                                                                                                                                                                     <  2:224156/250‑1 (MQ=255)
                                               ccccTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAAGCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGTTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                                                                                                                            <  2:219877/249‑1 (MQ=255)
                                               caccTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                                                                                                                            <  1:221193/251‑1 (MQ=255)
                                               caccTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                                                                                                                            <  2:246544/251‑1 (MQ=255)
                                                    gCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCATGCGCATGATGGAGCCTTTACCAATTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACTCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGTGTAATTCCGGATGGTCAAGATGTTCTGTGGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGtt                                                                                                                                                                                         <  1:287752/249‑1 (MQ=255)
                                                    gCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGtttt                                                                                                                                                                                       <  2:55385/251‑1 (MQ=255)
                                                    gCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAAGCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATTCTGTATGCTTATACAGTATCAAGTGtt                                                                                                                                                                                         <  2:13620/249‑1 (MQ=255)
                                                         aaGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGt                                                                                                                                                                                     >  2:44575/1‑248 (MQ=255)
                                                             gcgATATCCAGTGAAAGCGAACCGTTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTGTCTGTTTGTGTTCGTCGATAGCCATTTTTACTCCTGTCATGCCTGGTAATACCGGATGTTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAtt                                                                                                                                                                              <  1:227833/251‑1 (MQ=255)
                                                               gATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgt                                                                                                                                                                            >  1:166271/1‑251 (MQ=255)
                                                                    ccAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTAGACTGGATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGcc                                                                                                                                                                        <  1:466575/250‑1 (MQ=255)
                                                                      aGTGAAAGCGAACCGGTAGAGATGGTTTCACCATCCATGGACCGGTCTTCAACAACGCGCACGATGGGGCCTTTACCAAATTGTTTGTCAATCTGGCACAGTGCTGACGCCAACGCTTCATGTTTGTTTTCGTCGATAGCAATTTTACCTCCTGTCATGCCTGGTAATACCGGATGGTCATTATGTTCTGTTGAAGCAATTATACTGTATGCTCATCCAGTATAAAGTGTTTTGGAGAAATTGTTGCcac                                                                                                                                                                      >  2:436816/1‑250 (MQ=255)
                                                                          aaaGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGGTTTGTAGAAATTGTTGCCACAAGGt                                                                                                                                                                 <  2:277771/251‑1 (MQ=255)
                                                                                        gagaTGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttg                                                                                                                                                     >  2:351512‑M1/1‑246 (MQ=255)
                                                                                             ggTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactgg                                                                                                                                              >  2:82460‑M1/1‑241 (MQ=255)
                                                                                                    acaTCCATGGAACGGTCGTCACCCAGGCGCATGATGGAGACGTTCCCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACTCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGGCAATATGTTCGGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTATGTAGAAATTGGTGCCTCAAGGGCTGCAATGCttgggactggcctttct                                                                                                                                       <  2:72204‑M1/251‑18 (MQ=255)
                                                                                                             gAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgc                                                                                                                              <  1:351512‑M1/251‑27 (MQ=255)
                                                                                                                              gcgcATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaac                                                                                                             >  2:222147‑M1/1‑208 (MQ=255)
                                                                                                                                 cATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagc                                                                                                          <  1:47537‑M1/251‑47 (MQ=255)
                                                                                                                                       ggAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGCCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCACGTGTTTTGTAGAACTTGTTGCCACACGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcg                                                                                                    >  1:447831‑M1/1‑199 (MQ=255)
                                                                                                                                                           tctcAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgttt                                                                                 <  2:477251‑M1/250‑72 (MQ=255)
                                                                                                                                                                 tCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTCTACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggat                                                                          >  1:194956‑M1/1‑173 (MQ=255)
                                                                                                                                                                    ggCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttc                                                                       <  1:319125‑M1/251‑82 (MQ=255)
                                                                                                                                                                            gctgcCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCTATTTTTACTCCTGTCATGCCGGGTAATACCGGATCGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAACCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaa                                                                 >  1:4254/1‑249 (MQ=255)
                                                                                                                                                                               gccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCACTATGTTCTGTTGAAGCAATTATACTGTATGCTCATCCAGTATCAAGTGTTTTGTAGAAATTGTTGCCACaa                                                                                                                                                                    >  1:64669/1‑147 (MQ=255)
                                                                                                                                                                               gccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggccttgctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcggttgcaggatttcgcgcaaacg                                                              <  1:188785‑M1/249‑91 (MQ=255)
                                                                                                                                                                               gccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTGGTAGAAATTGTTGCCACaa                                                                                                                                                                    <  2:64669/147‑1 (MQ=255)
                                                                                                                                                                               gccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATGTTTACTCCTGTCATGCCGGGTAATACCGGATGGGCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGGCTGCAATGCttgagactggcctttctgacagatgctttcttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggc                                                            <  1:208878‑M1/251‑93 (MQ=255)
                                                                                                                                                                                 agccaACGCTTTCTGTATGTTTTCGTCTATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGCTGGTCCATATGTTCTGTTGAAGCAATTATACTGTAGGCTCATACAGTATCAAGTGTTTTGTAGAAATCGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgc                                                          >  2:295041‑M1/2‑157 (MQ=255)
                                                                                                                                                                                     aaCGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggc                                                      >  1:80210‑M1/1‑153 (MQ=255)
                                                                                                                                                                                             cTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacc                                              >  2:73980‑M1/1‑145 (MQ=255)
                                                                                                                                                                                                   gtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatc                                        >  1:157884‑M1/1‑139 (MQ=255)
                                                                                                                                                                                                     tttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacgtactcgcggaacagcgcttcgatatgcagccccggcgtttgcaggagttcgcgcaaacggcgcaggtcttcaaccggaatct                                       <  2:405341‑M1/250‑114 (MQ=255)
                                                                                                                                                                                                         gtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtgtgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcga                                  <  1:54851‑M1/251‑119 (MQ=255)
                                                                                                                                                                                                            gATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaccc                               >  2:365010‑M1/1‑130 (MQ=255)
                                                                                                                                                                                                                            ccTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggc               >  1:14940‑M1/1‑114 (MQ=255)
                                                                                                                                                                                                                                     ccGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcgtaacagcgcttcgatattcagctcctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgtt                           <  1:77418‑M1/230‑126 (MQ=255)
                                                                                                                                                                                                                                     ccGGGCAATACCGGCTGGTCAATATGTTCTGTTGAAGCAATTCTACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgtt                           >  2:77418‑M1/1‑105 (MQ=255)
                                                                                                                                                                                                                                      cGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggccaatttcacg     >  1:151368‑M1/1‑104 (MQ=255)
                                                                                                                                                                                                                                         gTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgaaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggccaatttcacgacc  >  2:448785‑M1/1‑101 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
GGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATCCGGCCCGGCGATACCACT  >  NC_000913/2823557‑2824040

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: