Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,796 A→G 6.2% intergenic (‑26/+60) 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,7960AG 6.2% 601.7 / 17.1 181intergenic (‑26/+60)yciF/yciGputative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand):  ref base A (79/90);  new base G (2/10);  total (81/100)
Fisher's exact test for biased strand distribution p-value = 6.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 185 total aligned reads displayed.

ATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCG  >  NC_000913/1315538‑1316053
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                 
            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)
                                   cGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAAtt                                                                                                                                                                                                                 <  2:275329/274‑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)
                                                                                                              aCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGa                                                                                                                                                                                                                                                                   <  6:204582/149‑1 (MQ=255)
                                                                                                              aCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGa                                                                                                                                                                                                                                                                   <  7:1330870/149‑1 (MQ=255)
                                                                                                                ttAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGc                                                                                                                                                                                                                                                                 <  6:219567/149‑1 (MQ=255)
                                                                                                                ttAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGc                                                                                                                                                                                                                                                                 >  5:219567/1‑149 (MQ=255)
                                                                                                                 tAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc                                                                                                                                                                                                                                                                <  6:422366/149‑1 (MQ=255)
                                                                                                                  aaTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa                                                                                                                                                                                                                                                               <  6:122322/149‑1 (MQ=255)
                                                                                                                  aaTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa                                                                                                                                                                                                                                                               <  8:883429/149‑1 (MQ=255)
                                                                                                                    tttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                                                                                                                                                                             <  6:182058/149‑1 (MQ=255)
                                                                                                                      tttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGTGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACg                                                                                                                                                                                                                                                           >  5:327271/1‑149 (MQ=255)
                                                                                                                           tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                                                                                                                                      >  7:49944/1‑149 (MQ=255)
                                                                                                                           tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTACGAGCCATACGGCCCg                                                                                                                                                                                                                                                      <  5:705024/149‑1 (MQ=255)
                                                                                                                           tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTATTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                                                                                                                                      >  8:952366/1‑149 (MQ=255)
                                                                                                                              gATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAg                                                                                                                                                                                                                                                   <  7:208561/149‑1 (MQ=255)
                                                                                                                               aTGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgt                                                                                                                                                                                                                                                  <  7:1340543/149‑2 (MQ=255)
                                                                                                                                 gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtga                                                                                                                                                                                                                                                <  7:1269856/149‑4 (MQ=255)
                                                                                                                                 gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                                                                                                                                <  8:1319145/149‑1 (MQ=255)
                                                                                                                                  ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgat                                                                                                                                                                                                                                               <  6:1453502/149‑5 (MQ=255)
                                                                                                                                  ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGc                                                                                                                                                                                                                                               >  6:900121/1‑149 (MQ=255)
                                                                                                                                  ttGCTCTTGCGAGTTTTGCCAGTGCACGGGTTAATTGTTTTTCTGCGCTGTAGGCATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgat                                                                                                                                                                                                                                               <  6:445842/149‑5 (MQ=255)
                                                                                                                                     ctctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAAt                                                                                                                                                                                                                                            >  5:779300/1‑149 (MQ=255)
                                                                                                                                        ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCaa                                                                                                                                                                                                                                             <  5:506301/145‑1 (MQ=255)
                                                                                                                                        ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCaa                                                                                                                                                                                                                                             >  6:506301/1‑145 (MQ=255)
                                                                                                                                        ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTAt                                                                                                                                                                                                                                         >  6:704970/1‑149 (MQ=255)
                                                                                                                                        ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTAt                                                                                                                                                                                                                                         <  6:1322305/149‑1 (MQ=255)
                                                                                                                                            gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                     >  6:231175/1‑149 (MQ=255)
                                                                                                                                            gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                     >  8:146038/1‑149 (MQ=255)
                                                                                                                                            gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                     <  7:256711/149‑1 (MQ=255)
                                                                                                                                            gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                     >  8:377833/1‑149 (MQ=255)
                                                                                                                                               ttttGCCAGTGACCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATc                                                                                                                                                                                                                                  >  7:828348/1‑149 (MQ=255)
                                                                                                                                                 ttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAg                                                                                                                                                                                                                                <  5:385088/149‑1 (MQ=255)
                                                                                                                                                    ccAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                                                                                                                                                                             >  7:721927/1‑149 (MQ=255)
                                                                                                                                                       gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGata                                                                                                                                                                                                                          <  7:688107/149‑1 (MQ=255)
                                                                                                                                                       gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGata                                                                                                                                                                                                                          <  7:929805/149‑1 (MQ=255)
                                                                                                                                                            cGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaat                                                                                                                                                                                                                     >  5:1195339/1‑149 (MQ=255)
                                                                                                                                                             gggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaata                                                                                                                                                                                                                    >  8:304712/1‑149 (MQ=255)
                                                                                                                                                              ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaataa                                                                                                                                                                                                                   <  7:1191450/149‑1 (MQ=255)
                                                                                                                                                              ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaataa                                                                                                                                                                                                                   <  5:851752/149‑1 (MQ=255)
                                                                                                                                                                ttAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAAtt                                                                                                                                                                                                                 <  7:615306/149‑1 (MQ=255)
                                                                                                                                                                      gTTTTTCTGCGCTGTAGGTATCTGAGAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGt                                                                                                                                                                                                           <  7:955403/149‑1 (MQ=255)
                                                                                                                                                                       tttttCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTg                                                                                                                                                                                                          <  5:1363574/149‑1 (MQ=255)
                                                                                                                                                                                gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                                                                                                                                                                                 <  7:382571/149‑1 (MQ=255)
                                                                                                                                                                                      ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCt                                                                                                                                                                                           >  7:1424536/1‑149 (MQ=255)
                                                                                                                                                                                        tATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGa                                                                                                                                                                                         >  5:968238/1‑149 (MQ=255)
                                                                                                                                                                                         aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                                                                                                                                        <  5:356697/149‑1 (MQ=255)
                                                                                                                                                                                          tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                                                                                                                                       >  5:518005/1‑149 (MQ=255)
                                                                                                                                                                                          tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                                                                                                                                       <  7:1077421/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)
                                                                                                                                                                                               aaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                                                                                                                                  >  5:1236849/1‑149 (MQ=255)
                                                                                                                                                                                               aaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                                                                                                                                  <  7:139936/149‑1 (MQ=255)
                                                                                                                                                                                                aGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTac                                                                                                                                                                                 <  8:347523/149‑1 (MQ=255)
                                                                                                                                                                                                 gCAGGTGAATAAATACATCTTCAATGGCCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacc                                                                                                                                                                                >  6:374085/1‑149 (MQ=255)
                                                                                                                                                                                                  cAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                                                                                                                                                                               <  6:968238/149‑1 (MQ=255)
                                                                                                                                                                                                        aataaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGACTGAATTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGTTTCAGATGCGCGTTTCGGATCATTGTTAAAATTACCGTAGCTATGCTGACCGCCTTTACGGCCTGCGTCTGTTGCATTCTAACGGTCTTCGGCGAAATTTcat                                           >  2:294529/1‑273 (MQ=255)
                                                                                                                                                                                                         ataaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCttt                                                                                                                                                                        <  5:1053421/149‑1 (MQ=255)
                                                                                                                                                                                                          taaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg                                                                                                                                                                       <  8:914589/149‑1 (MQ=255)
                                                                                                                                                                                                            aataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTCGCCTGATTTATTACCACCGCTttgtt                                                                                                                                                                     <  8:49944/149‑1 (MQ=255)
                                                                                                                                                                                                             ataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                     <  5:1245915/84‑1 (MQ=255)
                                                                                                                                                                                                             ataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                     >  6:1245915/1‑84 (MQ=255)
                                                                                                                                                                                                             ataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgttg                                                                                                                                                                    >  8:502927/1‑149 (MQ=255)
                                                                                                                                                                                                              taCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGa                                                                                                                                                                   <  7:947703/149‑1 (MQ=255)
                                                                                                                                                                                                               aCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGAc                                                                                                                                                                  <  5:900706/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)
                                                                                                                                                                                                                      cAATGGTCTTCATATTCATATTTTTCTCCATTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGc                                                                                                                                                                                             <  1:55334/115‑1 (MQ=255)
                                                                                                                                                                                                                       aaTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt                                                                                                                                                          <  5:742409/149‑1 (MQ=255)
                                                                                                                                                                                                                           gTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                                                                                                                                        >  7:578276/1‑115 (MQ=255)
                                                                                                                                                                                                                           gTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                                                                                                                                        <  8:578276/115‑1 (MQ=255)
                                                                                                                                                                                                                            ttttcatnttcataTTTTTTGCCATGTAAATCCCTGCCGGCCTTCCGGCCCGCAGATGCAAGTATCTTGATCAGCCGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTGCCACCGTTTTGTTGACTGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGc                       <  4:320004/273‑1 (MQ=255)
                                                                                                                                                                                                                              ttcatattcataTTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTCCCGCCGCTATGCTGCCCGCCTTTACGGCCTGCGTCGGATTCCTTCTCCCGGTCTTCGGCGGAATTTCCGGAACCACCACGATGtt                           >  2:320004/1‑269 (MQ=255)
                                                                                                                                                                                                                                catattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCt                                                                                                                                                 <  8:171071/149‑1 (MQ=255)
                                                                                                                                                                                                                                 atattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                                                                                                                                >  7:1442521/1‑149 (MQ=255)
                                                                                                                                                                                                                                          tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAtt                                                                                                                                                                                                >  2:374216/1‑92 (MQ=255)
                                                                                                                                                                                                                                          tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                                                                                                                                       <  5:224011/149‑1 (MQ=255)
                                                                                                                                                                                                                                               ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGc                                                                                                                                  <  8:721927/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTa                                                                                        <  4:289502/187‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                                                                                                                               <  7:11204/147‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                                                                                                                               >  8:11204/1‑147 (MQ=255)
                                                                                                                                                                                                                                                       aaaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtt                                                                                                                          <  8:16124/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        aaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt                                                                                                                         >  5:184044/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tCACTGTGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttttt                                                                                                                       >  8:443337/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          tCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttttt                                                                                                                       >  5:955855/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           cACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                                                                                                      <  6:48135/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                                                                                                                                     <  8:1358860/100‑1 (MQ=255)
                                                                                                                                                                                                                                                            tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                                                                                                                                     >  7:1358860/2‑101 (MQ=255)
                                                                                                                                                                                                                                                                cGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt                                                                                                                 >  6:1162405/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                                                                               <  6:334953/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                 
ATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCG  >  NC_000913/1315538‑1316053

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