Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 328,598 C→T V46I (GTC→ATC)  betB ← betaine aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,5980CT98.0% 165.4 / ‑4.2 51V46I (GTC→ATC) betBbetaine aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/26);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTT  >  NC_000913/328359‑328806
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                 
gCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTcgcg                                                                                                                                                                                                          <  2:50327/249‑1 (MQ=255)
             gggATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCAc                                                                                                                                                                                           >  1:213306/1‑251 (MQ=255)
                aTCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGCCGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGGCCGACCGATATCCCCTCGCCCGGGGGCCTGCCCggg                                                                                                                                                                                        >  2:155509/1‑250 (MQ=255)
                aTCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACggg                                                                                                                                                                                        >  2:218413/1‑250 (MQ=255)
                                aCTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCg                                                                                                                                                                       <  2:110671/251‑1 (MQ=255)
                                   ccAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGt                                                                                                                                                                      <  2:301785/249‑1 (MQ=255)
                                   ccAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCGTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGtt                                                                                                                                                                     <  2:259190/250‑1 (MQ=255)
                                     aGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGc                                                                                                                                                                  <  2:228201/251‑1 (MQ=255)
                                     aGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGATGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGc                                                                                                                                                                  <  1:376727/251‑1 (MQ=255)
                                        aCGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCgg                                                                                                                                                               <  2:99508/251‑1 (MQ=255)
                                            ccGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTa                                                                                                                                                           <  2:87352/251‑1 (MQ=255)
                                               cACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATg                                                                                                                                                        >  1:121138/1‑251 (MQ=255)
                                                 ccGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCCGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGt                                                                                                                                                      >  2:186089/1‑251 (MQ=255)
                                                         gATATCGAAGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAgg                                                                                                                                              >  1:103078/1‑251 (MQ=255)
                                                            atCGACGGTTGCGGGTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCCGCCCGGTCGCCCGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAggcgc                                                                                                                                           >  2:178664/1‑248 (MQ=255)
                                                            atCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATGTCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGCCCGGGTTAATGGTCTCGAAGGTg                                                                                                                                            <  2:379817/250‑1 (MQ=255)
                                                            agcGACGGTTGAGGTTTCCGAATATGCGTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTg                                                                                                                                            <  2:30004/248‑1 (MQ=255)
                                                                 cGGTGGAGTTTTCCGAATATGCTGTTCCGGTGTCGACGGTTTCCAGTTTTGCGAGTTCGGCATTGCGGTCACGCAGAAGATCAACGGCCCGACGCAGAATTCGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCTATATCCGCGCGCCCGGCGGCCTGCACGGGCTCCAGCACGTTACCGTTGGCCGGTTTAAGGGTCGCGAAGTTGCGACCg                                                                                                                                      <  2:259508/251‑1 (MQ=255)
                                                                     tGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCtgg                                                                                                                                  >  1:185603/1‑251 (MQ=255)
                                                                      gAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACTCAGAATATCAACTGCCCGACGCAGAATACGCGAGCGGTGCATGGCGGTCATCGACCCACAGATTTTTTGCCCCGGCTGGGCGCGTTTCACGGCGCGAACGATATCCTCCCGCCCGGCGGCCTGCAAGGTCGCCAGTCAGTAACCGTTCGCCGGGCTAAAAGTCTCGGAGGTTCGACCGCtggg                                                                                                                                 >  2:95958/1‑250 (MQ=255)
                                                                      gAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCtggg                                                                                                                                 >  1:29241/1‑250 (MQ=255)
                                                                        ggtgtCCGAATATGCTTTTCCGGTGTCTAGGGTTTCCAGGTTTGCGAGTTCGTCATTGCGTTCACGAAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGTCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCtggtgg                                                                                                                               <  2:229982/247‑1 (MQ=255)
                                                                          tttCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCg                                                                                                                             >  2:348196/1‑251 (MQ=255)
                                                                          tttCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGGGGCg                                                                                                                             >  1:119705/1‑251 (MQ=255)
                                                                                     gCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAacc                                                                                                                  <  1:338182/251‑1 (MQ=255)
                                                                                           ccGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCGGGTGGCGGAGGTATAACCACCAt                                                                                                             <  1:162182/250‑1 (MQ=255)
                                                                                             ggTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGGGGCGGAGGTATAACCACCATGt                                                                                                           >  1:8910/1‑250 (MQ=255)
                                                                                                   gAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCAGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGCTATCCTCGCGCCCGGCGGCCTGCACGTTCGCCACCACGTTACCGCTGGCCGGGGTAATGGTCTCGAAGTGGCTCCCCCTGGGGGCGGAGGTATAACCACCATGTATCTaa                                                                                                     >  2:15336/1‑250 (MQ=255)
                                                                                                    aGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGGGCGACCGCTGGGGGCGGAGGTATAACCACCATGGATATaa                                                                                                     >  2:376469/1‑249 (MQ=255)
                                                                                                        tttCCCGTTTTGCGAGTTCGTCCTTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGCTTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGCTCGATATTCTCGCGCCCGGCGGCCTGCCCGGTCGCCCGCACGTTACCGTTTGCCGGGTTAATGGTCTCGAAGGTGCGGCCGCTGGTGGCGGAGGGATAACCACCATGTTTATAAAgatgc                                                                                               >  2:165553/1‑251 (MQ=255)
                                                                                                         ttCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATTAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAgctgct                                                                                              >  1:48437/1‑251 (MQ=255)
                                                                                                                tttGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGc                                                                                       >  1:186359/1‑251 (MQ=255)
                                                                                                                           gTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGAc                                                                             >  1:47872/1‑250 (MQ=255)
                                                                                                                             cATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCTGGACATc                                                                          <  2:186359/251‑1 (MQ=255)
                                                                                                                                         gCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCgg                                                              >  1:206132/1‑251 (MQ=255)
                                                                                                                                                  cAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGtata                                                                                                        >  1:351730/1‑200 (MQ=255)
                                                                                                                                                  cAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGtata                                                                                                        <  2:351730/200‑1 (MQ=255)
                                                                                                                                                   aaCGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTg                                                    <  1:44629/251‑1 (MQ=255)
                                                                                                                                                      ggCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGGGGCGGAGGTATACCCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTACTCGGGGGGt                                                 >  2:121649/1‑251 (MQ=255)
                                                                                                                                                       gCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGt                                                 >  1:329337/1‑250 (MQ=255)
                                                                                                                                                          cGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGAt                                             <  1:35641/251‑1 (MQ=255)
                                                                                                                                                              gCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGtct                                                                                                                                                    >  1:277919/1‑144 (MQ=255)
                                                                                                                                                              gCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGtct                                                                                                                                                    <  2:277919/144‑1 (MQ=255)
                                                                                                                                                              gCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTAACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTTGGTAGATGCTg                                         >  1:198756/1‑251 (MQ=255)
                                                                                                                                                                 aaatgcgcgAGCGCTCCATGGCGTGCAGCGACGCCCCGATGTGTTGCCCCTGCTGGTCGCGTGTCACGGCGCGATCGATATCCTCGCGCCCGTCGGCCTGCACTGTCGCCGGCACGTTACCGTTGGCCGGGGTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTAGATAAAGGTGCTGGTCTGCCATTCGGGACATCACGGCTCCTCGGGTAATCGGTGGGTATATGCGGAGt                                      <  2:33608/246‑1 (MQ=255)
                                                                                                                                                                         gAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTTGGTAGATGCTGAGTGATAAAGt                              <  2:262510/251‑1 (MQ=255)
                                                                                                                                                                            cgcTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGc                           <  1:30946/251‑1 (MQ=255)
                                                                                                                                                                                    ggCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGTTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAg                    <  1:165553/250‑1 (MQ=255)
                                                                                                                                                                                    ggCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAgg                   <  1:178664/251‑1 (MQ=255)
                                                                                                                                                                                             cGACGCCAAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCTATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTTGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCTGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGTAATTAgc           <  2:135739/250‑1 (MQ=255)
                                                                                                                                                                                                     ggATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGATATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCGCGGTTAATCGGTGGGTAGATGCTGAGTGAGAAAGTGGCGGGTCAGTGAATTAGCGCGGGtttt  <  2:244193/250‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                 
GCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTT  >  NC_000913/328359‑328806

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