Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 339,376 G→A G351E (GGG→GAG)  yahG → DUF1116 domain‑containing protein YahG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,3760GA97.9% 157.7 / NA 48G351E (GGG→GAG) yahGDUF1116 domain‑containing protein YahG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/25);  minor base C (1/0);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGGGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGC‑GTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATT  >  NC_000913/339132‑339611
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
cAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACCTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGCAGATTCGGAGCTgg                                                                                                                                                                                                                                         >  2:192337/1‑250 (MQ=255)
 aGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGCTATAAGCCGGAAGATTCGGAGCTGGa                                                                                                                                                                                                                                        >  2:352506/1‑250 (MQ=255)
    gcgAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATc                                                                                                                                                                                                                                    <  1:275410/251‑1 (MQ=255)
       aaGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGc                                                                                                                                                                                                                                 <  1:97555/251‑1 (MQ=255)
                     gTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCa                                                                                                                                                                                                                    <  2:355268/250‑1 (MQ=255)
                             cagcGACTACTTCTCCGGCCCGACGTGGATGGCGCTGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCCGGTGATCGGCCCGATGATTGCCGGTTATAAGCCGGGAGATTCGGAGCGGGGTAGCGGCGACAGCGCCCTCACCGAAAccc                                                                                                                                                                                                           >  2:246296/1‑250 (MQ=255)
                                cGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACg                                                                                                                                                                                                        <  1:334886/251‑1 (MQ=255)
                                   ctactTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGGAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGt                                                                                                                                                                                                      >  2:135152/1‑250 (MQ=255)
                                       ttCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTAc                                                                                                                                                                                                         <  2:74413/243‑1 (MQ=255)
                                       ttCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTAc                                                                                                                                                                                                         >  1:74413/1‑243 (MQ=255)
                                       ttCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTgg                                                                                                                                                                                                 >  1:192421/1‑251 (MQ=255)
                                        tctcCGGCCCGACGTGGATGGCGATGTGTAAAGCCGCGATGGATGCGGCGCCCGGCATCGAATACAGCCCCGCGGTCACCCCCATGACGCGTAACGACGACGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGAGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGACGGAAGATTCGGCGCTGGATATCGGCGAAAGCGCCATCACCGAAACCTACGGTATTGGc                                                                                                                                                                                                >  1:38948/1‑251 (MQ=255)
                                          tcCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCgg                                                                                                                                                                                              >  1:268844/1‑251 (MQ=255)
                                                             cGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACAACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGACCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTTCTATGGCGaacg                                                                                                                                                                             >  1:372485/1‑246 (MQ=255)
                                                                     aaaGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTAt                                                                                                                                                                   >  1:396398/1‑251 (MQ=255)
                                                                          ggcgATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCGCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCTCCGAAACCTACGGTATTGGCGGATTTGCGATGGCGACCGCTCCCGCTAtcgtcg                                                                                                                                                              <  1:333006/251‑1 (MQ=255)
                                                                           gcgATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTAGCGTcgc                                                                                                                                                             <  1:283401/251‑1 (MQ=255)
                                                                             gatggatgCGGCGCATGGCATCGAATACAGCACCGTGGGCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGGTTGCCGGTTATAAGACGGAAGGTTAGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGGATTGGCGGGTTTGCTATGGCGAcagaccccccgatcagggccc                                                                                                                                                           >  2:385436/1‑231 (MQ=255)
                                                                                ggatgCGGTGCCTGGCCTCGAATACTGCACCGTGGCAACCTCCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCTAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCTCCCGCAATCGTCGCGCtg                                                                                                                                                         <  2:195986/250‑1 (MQ=255)
                                                                                   ggcggcgCATGGCACCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCTCCATCTCCGAAACCTACGTTAGTGGCGGATTTGCTATGGCGACCGCGGCCGCTATCGGCGCGCtggtgg                                                                                                                                                     <  2:376335/250‑1 (MQ=255)
                                                                                    gcggcgCATGCCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCtggtgg                                                                                                                                                     <  1:15215/250‑1 (MQ=255)
                                                                                         gcATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGAACGCGCCCGCCATCGTCGCCCTGGTGggcggcc                                                                                                                                               >  2:210617/1‑250 (MQ=255)
                                                                                               cATCGAATCCAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACggggg                                                                                                                                         >  1:361961/1‑248 (MQ=255)
                                                                                                aTCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGGGGGCGGCACggggg                                                                                                                                         >  1:81638/1‑247 (MQ=255)
                                                                                                       aCAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCCATCGTCGCGCTGGGGGGCGGCACCGTGGgagaaggt                                                                                                                                 >  2:150371/1‑249 (MQ=255)
                                                                                                                tGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTtc                                                                                                                        >  1:69030/1‑251 (MQ=255)
                                                                                                                    caccacCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTccc                                                                                                                    <  1:312036/251‑1 (MQ=255)
                                                                                                                        cccaTGGCGCGTACCGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCa                                                                                                                <  2:308045/250‑1 (MQ=255)
                                                                                                                         ccaTGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGGGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTAATGGCGGATTTGGTATGGGGACCGCCCCCCCCATCGTCGCGCTGGTGGGCGGGCCGGGGGAAGAAACTATTGATTTCCCCCGTCa                                                                                                                >  2:81422/1‑250 (MQ=255)
                                                                                                                                    aaCGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGAGTTCTCCCGTCAGATGCGCGAAAt                                                                                                    <  2:402728/251‑1 (MQ=255)
                                                                                                                                      cGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCa                                                                                                  <  2:109903/251‑1 (MQ=255)
                                                                                                                                            cGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTc                                                                                             <  2:396398/250‑1 (MQ=255)
                                                                                                                                            cGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGAACGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTTTGGCGACCGCGCCCGCCATCGGCGCGCCGGGGGGCGGCCCGGTGGAAGAAGCTATTGATTTCTCCCGTCAGGTGCGCGAAATCccccccc                                                                                            >  2:10452/1‑248 (MQ=255)
                                                                                                                                               gTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTg                                                                                         <  2:45725/251‑1 (MQ=255)
                                                                                                                                                            gTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGCGGGAGAAGCTTTTGATTTCTCCCGTCAGGTGCGCGCAATCACCCTCGGTGGAAACCCCAACgg                                                                            >  1:250663/1‑250 (MQ=255)
                                                                                                                                                                                  ggTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGttt                                                       <  1:181783/250‑1 (MQ=255)
                                                                                                                                                                                            cccggcGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTTATTtc                                                                                                                        >  1:178852/1‑175 (MQ=255)
                                                                                                                                                                                            cccggcGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTtc                                                                                                                        <  2:178852/175‑1 (MQ=255)
                                                                                                                                                                                              cggcGCAGCAGGTGATCGGCCCGATGGTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGC‑GTGCCg                                         <  1:331111/251‑1 (MQ=255)
                                                                                                                                                                                                       aggtgagcgggacGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATTGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCTCCATTCCGCTGCGCGGTTTTATGGGC‑GTGCCGGCGGCAATc                                <  2:162774/239‑1 (MQ=255)
                                                                                                                                                                                                              cGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTTCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGGGGAAGAAGCTATTGGTTTCTCCCGTCAGATGCGCGGAATCACCCCCGGTGGAAACCCCCACGTCACCATTCCGCTGCGCGGTTTTTTGGGC‑GTGCCCTCGGCAATCGACATc                          >  2:119695/1‑250 (MQ=255)
                                                                                                                                                                                                                  ccGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTAGCGCTGGTGGGCGGCCCGGTGGAAGAAGCTATTGATTTCTCCCGTCCGATGCGCGCAATCACCCTCGGTGAAAACCCCAACGACACCATTCCGCTGCTCGGTTTTATGGGC‑GTGCCGTCGGCAATCGACATCACCcg                     >  1:399092/1‑251 (MQ=255)
                                                                                                                                                                                                                    gATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAAACGACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGGGGGCGGCACGGGGGAAGAAACGAATGAATTCTGCCGTCAGATGGGCGAAAACACCCGCGGGGAAAAACCCAAAGTCACCAGTCCGCTGAGCGGGTTTATGGGCAGT‑CCGTCGGGAAATGACATCACCcgcg                   >  2:224014/1‑251 (MQ=255)
                                                                                                                                                                                                                       gTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGC‑GTGCCGTCGGCAATCGACATCACCCGCGTgg                <  2:268844/251‑1 (MQ=255)
                                                                                                                                                                                                                       gTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGC‑GTGCCGTCGGCAATCGACATCACCCGCGTgg                <  1:150590/251‑1 (MQ=255)
                                                                                                                                                                                                                          tGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTGATGGGC‑GGGCCGTCGGCAATCGACAGCACCCGCGTGGGt              <  2:237967/250‑1 (MQ=255)
                                                                                                                                                                                                                            ccGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGC‑GTGCCGTCGGCAATCGACATCACCCGCGTGGGTagc           <  1:5645/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGC‑GTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCAtt  <  1:210617/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
CAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGGGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGC‑GTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATT  >  NC_000913/339132‑339611

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: ‑