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

GGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATCCGGCCCGGCGATACCACT  >  NC_000913/2823557‑2824040
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
ggAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑Caa                                                                                                                                                                                                                                           <  2:357068/251‑1 (MQ=255)
 gAAGCTTCCGGTCCGTAGATTGCGACGATACGGCCCAGCGGCAGACCACCTGCCCCAAGCGCGATATCCATTGAAAGCGAACCGGTAGAGATTGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCGTTCTGTTTGGTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAat                                                                                                                                                                                                                                          <  2:232588/251‑1 (MQ=255)
     gttCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGt                                                                                                                                                                                                                                      <  2:267126/250‑1 (MQ=255)
           gTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGCTATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCCCATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGtt                                                                                                                                                                                                                                >  1:477251/1‑251 (MQ=255)
             ccGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGTAGCCTTTACCAAATTGTTTCTCAATCTTTCCCAGTTCTGCCGCCAACGCTTGCTGGTTGTTTTCGTCGATAGCCAGGTGTACTCCTGGCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTg                                                                                                                                                                                                                               <  1:431389/250‑1 (MQ=255)
               gTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAg                                                                                                                                                                                                                            <  2:57474/251‑1 (MQ=255)
               gTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAg                                                                                                                                                                                                                            <  2:401906/251‑1 (MQ=255)
                tAGATTTCGACTATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTGCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAAGGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCGGTCATGCCGGGTAATACCGGATGCT‑CAATATGTTCTGTTGAAg                                                                                                                                                                                                                            <  2:437647/250‑1 (MQ=255)
                           gaTACGGCACATCGGCCGACCACCTGCCCCAAGCGCGATATCCAGTGAACGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTACTACCGGCTGGT‑CAATATGTTCTGTTGAAGCAATTATACttt                                                                                                                                                                                                                >  2:388182/1‑249 (MQ=255)
                            aTACGGCCCATCGGCAGACCGCCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTGTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGTTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTa                                                                                                                                                                                                               <  1:189444/251‑1 (MQ=255)
                            aTACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGCTATCCAGTGACAGCGAACCGGTAGAGATGGTTTCCCCATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGt                                                                                                                                                                                                                >  1:29313/1‑250 (MQ=255)
                                     aTCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGCCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTACTACCGGATGGT‑CAATCTGTTCTGTTGAAGCAATTATACTGTATGCTCATAc                                                                                                                                                                                                      >  2:2923/1‑251 (MQ=255)
                                     aTCGGAAGAACACCTGCACCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGCTGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATAc                                                                                                                                                                                                      >  1:357395/1‑251 (MQ=255)
                                      ccGGCAGACCACCTGCGCCAAGCGCGATATCCCGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTGCACCCAGGCGCATGATGGAGCCTTTACGAAATTCTTTCTCAATCTTGCCCACTGCCTCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCTATTTTTACTCCTGTCAGGCCGGGCGTTACCGGAT‑GTGCAGTATGGGCTGGTGAAGCGATGATACTTTATGCTCATACa                                                                                                                                                                                                     <  2:224156/250‑1 (MQ=255)
                                               ccccTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAAGCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGTTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                                                                                                                            <  2:219877/249‑1 (MQ=255)
                                               caccTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                                                                                                                            <  1:221193/251‑1 (MQ=255)
                                               caccTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                                                                                                                            <  2:246544/251‑1 (MQ=255)
                                                    gCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCATGCGCATGATGGAGCCTTTACCAATTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACTCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGTGTAATTCCGGATGGT‑CAAGATGTTCTGTGGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGtt                                                                                                                                                                                         <  1:287752/249‑1 (MQ=255)
                                                    gCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGtttt                                                                                                                                                                                       <  2:55385/251‑1 (MQ=255)
                                                    gCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAAGCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATTCTGTATGCTTATACAGTATCAAGTGtt                                                                                                                                                                                         <  2:13620/249‑1 (MQ=255)
                                                         aaGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGt                                                                                                                                                                                     >  2:44575/1‑248 (MQ=255)
                                                             gcgATATCCAGTGAAAGCGAACCGTTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTGTCTGTTTGTGTTCGTCGATAGCCATTTTTACTCCTGTCATGCCTGGTAATACCGGATGTT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAtt                                                                                                                                                                              <  1:227833/251‑1 (MQ=255)
                                                               gATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgt                                                                                                                                                                            >  1:166271/1‑251 (MQ=255)
                                                                    ccAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTAGACTGGATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGcc                                                                                                                                                                        <  1:466575/250‑1 (MQ=255)
                                                                      aGTGAAAGCGAACCGGTAGAGATGGTTTCACCATCCATGGACCGGTCTTCAACAACGCGCACGATGGGGCCTTTACCAAATTGTTTGTCAATCTGGCACAGTGCTGACGCCAACGCTTCATGTTTGTTTTCGTCGATAGCAATTTTACCTCCTGTCATGCCTGGTAATACCGGATGGT‑CATTATGTTCTGTTGAAGCAATTATACTGTATGCTCATCCAGTATAAAGTGTTTTGGAGAAATTGTTGCcac                                                                                                                                                                      >  2:436816/1‑250 (MQ=255)
                                                                          aaaGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGGTTTGTAGAAATTGTTGCCACAAGGt                                                                                                                                                                 <  2:277771/251‑1 (MQ=255)
                                                                                        gagaTGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttg                                                                                                                                                     >  2:351512‑M1/1‑246 (MQ=255)
                                                                                             ggTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactgg                                                                                                                                              >  2:82460‑M1/1‑241 (MQ=255)
                                                                                                    acaTCCATGGAACGGTCGTCACCCAGGCGCATGATGGAGACGTTCCCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACTCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGG‑CAATATGTTCGGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTATGTAGAAATTGGTGCCTCAAGGGCTGCAATGCttgggactggcctttct                                                                                                                                       <  2:72204‑M1/251‑18 (MQ=255)
                                                                                                             gAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgc                                                                                                                              <  1:351512‑M1/251‑27 (MQ=255)
                                                                                                                              gcgcATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaac                                                                                                             >  2:222147‑M1/1‑208 (MQ=255)
                                                                                                                                 cATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagc                                                                                                          <  1:47537‑M1/251‑47 (MQ=255)
                                                                                                                                       ggAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGC‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCACGTGTTTTGTAGAACTTGTTGCCACACGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcg                                                                                                    >  1:447831‑M1/1‑199 (MQ=255)
                                                                                                                                                           tctcAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgttt                                                                                 <  2:477251‑M1/250‑72 (MQ=255)
                                                                                                                                                                 tCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTCTACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggat                                                                          >  1:194956‑M1/1‑173 (MQ=255)
                                                                                                                                                                    ggCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttc                                                                       <  1:319125‑M1/251‑82 (MQ=255)
                                                                                                                                                                            gctgcCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCTATTTTTACTCCTGTCATGCCGGGTAATACCGGATCGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAACCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaa                                                                 >  1:4254/1‑249 (MQ=255)
                                                                                                                                                                               gccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CACTATGTTCTGTTGAAGCAATTATACTGTATGCTCATCCAGTATCAAGTGTTTTGTAGAAATTGTTGCCACaa                                                                                                                                                                    >  1:64669/1‑147 (MQ=255)
                                                                                                                                                                               gccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggccttgctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcggttgcaggatttcgcgcaaacg                                                              <  1:188785‑M1/249‑91 (MQ=255)
                                                                                                                                                                               gccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTGGTAGAAATTGTTGCCACaa                                                                                                                                                                    <  2:64669/147‑1 (MQ=255)
                                                                                                                                                                               gccgccAACGCTTTCTGTTTGTTTTCGTCGATAGCCATGTTTACTCCTGTCATGCCGGGTAATACCGGATGGG‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGGCTGCAATGCttgagactggcctttctgacagatgctttcttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggc                                                            <  1:208878‑M1/251‑93 (MQ=255)
                                                                                                                                                                                 agccaACGCTTTCTGTATGTTTTCGTCTATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGCTGGT‑CCATATGTTCTGTTGAAGCAATTATACTGTAGGCTCATACAGTATCAAGTGTTTTGTAGAAATCGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgc                                                          >  2:295041‑M1/2‑157 (MQ=255)
                                                                                                                                                                                     aaCGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggc                                                      >  1:80210‑M1/1‑153 (MQ=255)
                                                                                                                                                                                             cTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacc                                              >  2:73980‑M1/1‑145 (MQ=255)
                                                                                                                                                                                                   gtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatc                                        >  1:157884‑M1/1‑139 (MQ=255)
                                                                                                                                                                                                     tttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacgtactcgcggaacagcgcttcgatatgcagccccggcgtttgcaggagttcgcgcaaacggcgcaggtcttcaaccggaatct                                       <  2:405341‑M1/250‑114 (MQ=255)
                                                                                                                                                                                                         gtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtgtgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcga                                  <  1:54851‑M1/251‑119 (MQ=255)
                                                                                                                                                                                                            gATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaccc                               >  2:365010‑M1/1‑130 (MQ=255)
                                                                                                                                                                                                                            ccTGTCATGCCGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggc               >  1:14940‑M1/1‑114 (MQ=255)
                                                                                                                                                                                                                                     ccGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcgtaacagcgcttcgatattcagctcctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgtt                           <  1:77418‑M1/230‑126 (MQ=255)
                                                                                                                                                                                                                                     ccGGGCAATACCGGCTGGT‑CAATATGTTCTGTTGAAGCAATTCTACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgtt                           >  2:77418‑M1/1‑105 (MQ=255)
                                                                                                                                                                                                                                      cGGGTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggccaatttcacg     >  1:151368‑M1/1‑104 (MQ=255)
                                                                                                                                                                                                                                         gTAATACCGGATGGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgaaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggccaatttcacgacc  >  2:448785‑M1/1‑101 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
GGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGT‑CAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATCCGGCCCGGCGATACCACT  >  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: ‑