Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,803 C→T 6.7% intergenic (‑33/+53) yciF ← / ← yciG putative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,315,8030CT 6.7% 571.7 / 22.4 178intergenic (‑33/+53)yciF/yciGputative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand):  ref base C (77/89);  new base T (3/9);  total (80/98)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

Only 100 of 179 total aligned reads displayed.

ATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATGAG  >  NC_000913/1315538‑1316076
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                 
            ttCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGCGGTGCGCATGACAAGCCTGACTTAATCTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCCGTGCCCGGGTTAATTGTTTTTCTGCGCCGTAGGTATCTGAAAGCAGGTGAATAAATACATTTTAATTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGACATACGGCCCGCAGATGCACTTATca                                                                                                                                                                                                                                                              >  2:160599/1‑274 (MQ=255)
            ttCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACTTCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCCCTGCGAGCCATACGGCCCGCAACTGAAATTATCt                                                                                                                                                                                                                                                              >  4:275329/1‑275 (MQ=255)
              cATGGGTTTAATTTTCCGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCGGTCCATGGGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                                            <  4:249482/275‑1 (MQ=255)
                      tAAGTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGGGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                                                                                                                                                                                                    <  2:85508/275‑1 (MQ=255)
                                                         tCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGa                                                                                                                                                                                                                                                                         >  2:113928/1‑219 (MQ=255)
                                                             tACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCa                                                                                                                                                                                                                                                                     <  4:113928/219‑1 (MQ=255)
                                                                       cTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                                                                                                                                                                                                  >  4:143866/1‑244 (MQ=255)
                                                                           ccATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAAt                                                                                                                                                                                                                              <  2:143866/244‑1 (MQ=255)
                                                                                                                     ttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATAc                                                                                                                                                                                                                                                                                   <  8:541495/149‑1 (MQ=255)
                                                                                                                      tttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGTGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACg                                                                                                                                                                                                                                                                                  >  5:327271/1‑149 (MQ=255)
                                                                                                                           tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                                                                                                                                                             >  7:49944/1‑149 (MQ=255)
                                                                                                                           tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTATTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                                                                                                                                                             >  8:952366/1‑149 (MQ=255)
                                                                                                                           tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCACGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGGCATCATATTCATATTTTTCTCCAGTGAAATCACAGCGAGCCATACGGCCCg                                                                                                                                                                                                                                                                             >  8:912898/1‑149 (MQ=255)
                                                                                                                                  ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgat                                                                                                                                                                                                                                                                      <  6:967256/149‑5 (MQ=255)
                                                                                                                                  ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgat                                                                                                                                                                                                                                                                      <  6:1453502/149‑5 (MQ=255)
                                                                                                                                  ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGc                                                                                                                                                                                                                                                                      >  6:900121/1‑149 (MQ=255)
                                                                                                                                  ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGc                                                                                                                                                                                                                                                                      <  5:1231505/149‑1 (MQ=255)
                                                                                                                                  ttGCTCTTGCGAGTTTTGCCAGTGCACGGGTTAATTGTTTTTCTGCGCTGTAGGCATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgat                                                                                                                                                                                                                                                                      <  6:445842/149‑5 (MQ=255)
                                                                                                                                     ctctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                                                                                                                                                                                                   >  5:779300/1‑149 (MQ=255)
                                                                                                                                        ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCaa                                                                                                                                                                                                                                                                    <  5:506301/145‑1 (MQ=255)
                                                                                                                                        ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTAt                                                                                                                                                                                                                                                                <  6:1322305/149‑1 (MQ=255)
                                                                                                                                        ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTAt                                                                                                                                                                                                                                                                >  6:704970/1‑149 (MQ=255)
                                                                                                                                            gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                                            >  8:146038/1‑149 (MQ=255)
                                                                                                                                            gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                                            >  6:231175/1‑149 (MQ=255)
                                                                                                                                            gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                                            >  7:1340266/1‑149 (MQ=255)
                                                                                                                                            gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                                            >  8:377833/1‑149 (MQ=255)
                                                                                                                                              gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTAt                                                                                                                                                                                                                                                          <  5:723431/149‑1 (MQ=255)
                                                                                                                                               ttttGCCAGTGACCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATc                                                                                                                                                                                                                                                         >  7:828348/1‑149 (MQ=255)
                                                                                                                                                tttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCa                                                                                                                                                                                                                                                        >  5:948388/1‑149 (MQ=255)
                                                                                                                                                 ttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAg                                                                                                                                                                                                                                                       <  5:385088/149‑1 (MQ=255)
                                                                                                                                                   gCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCt                                                                                                                                                                                                                                                     <  6:408574/149‑1 (MQ=255)
                                                                                                                                                    ccAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                                                                                                                                                                                                    >  7:721927/1‑149 (MQ=255)
                                                                                                                                                    ccAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                                                                                                                                                                                                    <  6:182799/149‑1 (MQ=255)
                                                                                                                                                       gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgt                                                                                                                                                                                                                                                                         >  5:893456/1‑124 (MQ=255)
                                                                                                                                                       gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGata                                                                                                                                                                                                                                                 <  7:929805/149‑1 (MQ=255)
                                                                                                                                                         gCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATAtt                                                                                                                                                                                                                                               >  7:1442270/1‑149 (MQ=255)
                                                                                                                                                              ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaataa                                                                                                                                                                                                                                          <  7:1191450/149‑1 (MQ=255)
                                                                                                                                                                   aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAc                                                                                                                                                                                                                                     <  5:221247/149‑1 (MQ=255)
                                                                                                                                                                   aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAc                                                                                                                                                                                                                                     <  7:108983/149‑1 (MQ=255)
                                                                                                                                                                   aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAc                                                                                                                                                                                                                                     <  8:1340266/149‑1 (MQ=255)
                                                                                                                                                                     tGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACgg                                                                                                                                                                                                                                   <  7:1402773/149‑1 (MQ=255)
                                                                                                                                                                      gTTTTTCTGCGCTGTAGGTATCTGAGAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                                                                                                                                                                                                  <  7:955403/149‑1 (MQ=255)
                                                                                                                                                                      gTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                                                                                                                                                                                                  <  7:227598/149‑1 (MQ=255)
                                                                                                                                                                       tttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTg                                                                                                                                                                                                                                 <  5:1363574/149‑1 (MQ=255)
                                                                                                                                                                                gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                                                                                                                                                                                                        <  7:382571/149‑1 (MQ=255)
                                                                                                                                                                                    tAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGc                                                                                                                                                                                                                    >  5:1112364/1‑149 (MQ=255)
                                                                                                                                                                                        tATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGa                                                                                                                                                                                                                >  5:968238/1‑149 (MQ=255)
                                                                                                                                                                                         aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                                                                                                                                                               <  7:717138/149‑1 (MQ=255)
                                                                                                                                                                                         aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                                                                                                                                                               <  5:356697/149‑1 (MQ=255)
                                                                                                                                                                                          tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                                                                                                                                                              >  5:518005/1‑149 (MQ=255)
                                                                                                                                                                                          tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                                                                                                                                                              <  7:657403/149‑1 (MQ=255)
                                                                                                                                                                                           cTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAttt                                                                                                                                                                                                             >  7:1253888/1‑149 (MQ=255)
                                                                                                                                                                                             gAAAGCAGGTGAATAAATACATCTTCGATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttat                                                                                                                                                                                                           >  7:810902/1‑149 (MQ=255)
                                                                                                                                                                                              aaaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTa                                                                                                                                                                                                                                                           >  1:31065/1‑100 (MQ=255)
                                                                                                                                                                                              aaaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                                                                                                                                                                                          <  6:848872/149‑1 (MQ=255)
                                                                                                                                                                                              aaaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatt                                                                                                                                                                                                          <  6:337153/149‑1 (MQ=255)
                                                                                                                                                                                               aaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                                                                                                                                                         <  7:139936/149‑1 (MQ=255)
                                                                                                                                                                                                aGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTac                                                                                                                                                                                                        <  8:347523/149‑1 (MQ=255)
                                                                                                                                                                                                        aataaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGACTGAATTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGTTTCAGATGCGCGTTTCGGATCATTGTTAAAATTACCGTAGCTATGCTGACCGCCTTTACGGCCTGCGTCTGTTGCATTCTAACGGTCTTCGGCGAAATTTcat                                                                  >  2:294529/1‑273 (MQ=255)
                                                                                                                                                                                                         ataaataCATCTTTAATGGTTGTCATATGCATATTTTTCTCCAGTGAAAGCACTGCAGGCCATCCGGCCCGCAAATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTg                                                                 <  2:47242/275‑1 (MQ=255)
                                                                                                                                                                                                         ataaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCttt                                                                                                                                                                                               <  5:1053421/149‑1 (MQ=255)
                                                                                                                                                                                                           aaataCATCTTAAAGTGTCTTCATATTCATATTTTTATCCAGTAAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa                                                               <  2:236103/275‑1 (MQ=255)
                                                                                                                                                                                                            aataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTCGCCTGATTTATTACCACCGCTttgtt                                                                                                                                                                                            <  8:49944/149‑1 (MQ=255)
                                                                                                                                                                                                             ataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                                            >  6:1245915/1‑84 (MQ=255)
                                                                                                                                                                                                              taCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGa                                                                                                                                                                                          <  7:947703/149‑1 (MQ=255)
                                                                                                                                                                                                                 aTCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCg                                                                                                                                                                                       >  6:519626/1‑149 (MQ=255)
                                                                                                                                                                                                                  tCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                                                                                                                                                                                                        >  3:55334/1‑115 (MQ=255)
                                                                                                                                                                                                                  tCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTCTTCACGGTCTTCGGCGAAATTTCCTTAtcca                                                            >  4:342101/1‑271 (MQ=255)
                                                                                                                                                                                                                    ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCTTGATTTATTACCACCGCTTGGTTGACCGCCt                                                                                                                                                                                    >  5:301401/1‑149 (MQ=255)
                                                                                                                                                                                                                    ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCt                                                                                                                                                                                    >  5:1307902/1‑149 (MQ=255)
                                                                                                                                                                                                                      cAATGGTCTTCATATTCATATTTTTCTCCATTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGc                                                                                                                                                                                                                    <  1:55334/115‑1 (MQ=255)
                                                                                                                                                                                                                           gtattcatatgaataTTTTTCTCCATTGAAATCCCTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAATACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCgg                                               <  2:333249/272‑1 (MQ=255)
                                                                                                                                                                                                                           gTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                                                                                                                                                               <  8:578276/115‑1 (MQ=255)
                                                                                                                                                                                                                            ttttcatnttcataTTTTTTGCCATGTAAATCCCTGCCGGCCTTCCGGCCCGCAGATGCAAGTATCTTGATCAGCCGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTGCCACCGTTTTGTTGACTGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGc                                              <  4:320004/273‑1 (MQ=255)
                                                                                                                                                                                                                                    ttcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAg                                                                                                                                                                    >  7:1057848/1‑149 (MQ=255)
                                                                                                                                                                                                                                     tcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGa                                                                                                                                                                   >  8:169415/1‑149 (MQ=255)
                                                                                                                                                                                                                                       atattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATg                                                                                                                                                                 >  7:1358223/1‑149 (MQ=255)
                                                                                                                                                                                                                                        tattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                                                                                                                                                                >  5:684183/1‑149 (MQ=255)
                                                                                                                                                                                                                                          tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                                                                                                                                                              <  6:500487/149‑1 (MQ=255)
                                                                                                                                                                                                                                               ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGc                                                                                                                                                         <  8:721927/149‑1 (MQ=255)
                                                                                                                                                                                                                                               ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGc                                                                                                                                                         >  8:271506/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   aGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                                                                                                                                                     >  7:16124/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   aGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTa                                                                                                               <  4:289502/187‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                                                                                                                                                      >  8:11204/1‑147 (MQ=255)
                                                                                                                                                                                                                                                       aaaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtt                                                                                                                                                 <  6:301035/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt                                                                                                                                                <  7:599762/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt                                                                                                                                                >  5:184044/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tCACTGTGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttttt                                                                                                                                              >  8:443337/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt                                                                                                                                                                                 >  7:317275/1‑114 (MQ=255)
                                                                                                                                                                                                                                                          tCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt                                                                                                                                                                                 <  8:317275/114‑1 (MQ=255)
                                                                                                                                                                                                                                                          tCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttttt                                                                                                                                              >  5:955855/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           cACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                                                                                                                             <  6:48135/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                                                                                                                                                            <  8:1358860/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                cgggCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                                                                                                                                                     >  7:1291486/4‑136 (MQ=255)
                                                                                                                                                                                                                                                                cgggCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGAt                                                                                                                                                     <  8:1291486/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                cGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt                                                                                                                                        >  6:1162405/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       tACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                                                                                                                                                            <  6:1240257/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                                                                                                                                                            >  5:1240257/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                        aCGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACcgccgc                                                                                                                                <  5:231175/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aCGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATGAg  >  2:278275/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                 
ATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATGAG  >  NC_000913/1315538‑1316076

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