Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,406,098 C→T 5.8% D40D (GAC→GAT smrA → DNA endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,406,0980CT 5.8% 447.7 / 13.9 139D40D (GAC→GATsmrADNA endonuclease
Reads supporting (aligned to +/- strand):  ref base C (69/62);  new base T (4/4);  total (73/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01

Only 100 of 143 total aligned reads displayed.

TTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCG  >  NC_000913/1405864‑1406365
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                           
tttGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                                                                                                                                                     >  2:393742/1‑275 (MQ=255)
tttGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCAAACGCGTAACCAACGTGCGCCGCAGCGTAGAGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                                                                                                                                                     >  1:34750/1‑275 (MQ=255)
  tGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCAACGGATTtccccac                                                                                                                                                                                                                                   >  4:375815/1‑270 (MQ=255)
              cGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                                                                                                                                                       >  1:117349/1‑275 (MQ=255)
              cGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCAATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCCACACGCTGAAGCTTTACAATTTCCTCACCCCCGGATTTCTCGACATCACCCCACTa                                                                                                                                                                                                                       >  4:518725/1‑275 (MQ=255)
                                                 tattatgttgttGGACCGCCCCGGTTGGNTNTTCATCTCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGNCGACAAATCGCTTTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGNCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCa                                                                                                                                                                                    <  3:34750/269‑1 (MQ=255)
                                                         tgttgGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCNACCGATGTCCACTGGCATCCAACNCGTANCCAACGTGCGCCGCAGCGTACCGACCCGCTGCAGCTTGATAATTTCCTcaccc                                                                                                                                                                                                                                                 >  4:540164/1‑205 (MQ=255)
                                                                 cgaccCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                                                                                                                                                    <  2:375816/275‑1 (MQ=255)
                                                                   accCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGTGTGCAACCTGGGGTGCTGGATaa                                                                                                                                                                  >  2:124105/1‑275 (MQ=255)
                                                                   accCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATaa                                                                                                                                                                  >  2:215690/1‑275 (MQ=255)
                                                                     ccGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGc                                                                                                                                                                >  4:213645/1‑275 (MQ=255)
                                                                       ggTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCggg                                                                                                                                                                                              <  2:504780/243‑1 (MQ=255)
                                                                          cgggttttttgttCGCGTTTTTAATTGGNGNGGGATGTAAAACGAACCGTGACGACAAATCGCTGGTTCTTGCCGCNATGGAAGATGTCCAGCCGCTGAAACGTTCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGNCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcgc                                                                                                                                                           <  3:41667/263‑1 (MQ=255)
                                                                                        gcgcTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGaca                                                                                                                                                                                                                                                                            >  5:925313/1‑148 (MQ=255)
                                                                                        gcgcTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGaca                                                                                                                                                                                                                                                                            <  5:824659/148‑1 (MQ=255)
                                                                                         cgcTCTTAATTGCGGGCGGAAGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCAATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCAACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGCATCGACACGCTGCAGACTGACAATTTCCTCACCCACTGATTTCTCGACCTCATCCCACTACGCAAGCCGCTGGAGTTTCGGCTGGAAGGGTTGCAACATGGGGTGCGGGATAAGCGCCGCAGACGTaaacaccca                                                                                                                                            >  4:6357/1‑269 (MQ=255)
                                                                                                     cGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGa                                                                                                                                                                                                                                                               >  6:981765/1‑148 (MQ=255)
                                                                                                      gggAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGc                                                                                                                                >  2:432444/1‑274 (MQ=255)
                                                                                                       ggaggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCtt                                                                                                                                                                                                                                                                 <  8:1055243/144‑1 (MQ=255)
                                                                                                               tAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTc                                                                                                                                                                                                                                                     <  6:925313/148‑1 (MQ=255)
                                                                                                                    tGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGc                                                                                                                                                                                                                                                                      >  7:968024/1‑126 (MQ=255)
                                                                                                                    tGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                                                                                                                                                                                                    >  7:406862/1‑144 (MQ=255)
                                                                                                                     gAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCCGCGCAGTGGTAAATATCCGCAACAGGCAGGCCTGAATCTTTTGCGcc                                                                                                                >  4:93500/1‑275 (MQ=255)
                                                                                                                      aaCCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                                                                                                                                                                                                    <  6:16231/142‑1 (MQ=255)
                                                                                                                      aaCCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGc                                                                                                                                                                >  2:362583/1‑226 (MQ=255)
                                                                                                                       aCCTTACCGGCAAACCGCGTTTTCTTTACGCCATGGAAGATGTCCGGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCa                                                                                                               <  4:393741/274‑1 (MQ=255)
                                                                                                                         ctggacgaaCAATCCCTGTTTCTTGCGGCCCAGGAAGTTTTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGc                                                                                                             <  4:272266/271‑1 (MQ=255)
                                                                                                                          ttGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                                                                                                                                                                                                  <  2:563960/124‑1 (MQ=255)
                                                                                                                           tGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTgcgc                                                                                                                 >  2:485124/1‑268 (MQ=255)
                                                                                                                               gacAAATCGCTTTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGc                                                                                                             <  4:485123/268‑1 (MQ=255)
                                                                                                                                      cGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatcat                                                                                                                                                                                                                              >  7:1043743/1‑148 (MQ=255)
                                                                                                                                           tttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGaca                                                                                                                                                                                                                                                                            >  6:809652/1‑97 (MQ=255)
                                                                                                                                            ttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                                                                                                                                                                                                   <  6:192866/105‑1 (MQ=255)
                                                                                                                                            ttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                                                                                                                                                                                                   >  5:192866/1‑105 (MQ=255)
                                                                                                                                                tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTtctc                                                                                                                                                                                                                                      >  5:758434/1‑130 (MQ=255)
                                                                                                                                                          gAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                                                                                                                                                                                           <  6:396426/131‑1 (MQ=255)
                                                                                                                                                           aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCt                                                                                                                                                                                                                                                      >  7:203921/1‑103 (MQ=255)
                                                                                                                                                           aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCt                                                                                                                                                                                                                                                      <  8:203921/103‑1 (MQ=255)
                                                                                                                                                           aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGt                                                                                                                                                                                                                    <  5:233830/137‑1 (MQ=255)
                                                                                                                                                             gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCt                                                                                                                                                                                                                                                                        <  7:282589/83‑1 (MQ=255)
                                                                                                                                                             gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                                                                                                                                                                                           >  8:299164/1‑128 (MQ=255)
                                                                                                                                                              aTGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGtt                                                                                                                                                                                                      >  5:652683/1‑148 (MQ=255)
                                                                                                                                                                gTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCTGGAAGGGTTGCAACATTGGGTTCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCTAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTTTTTAGTTTTTTTcc                                                                     >  4:429715/1‑274 (MQ=255)
                                                                                                                                                                  ccAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                                                                                                                                                                                                      <  8:1062219/80‑1 (MQ=255)
                                                                                                                                                                  ccAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                                                                                                                                                                                                      >  7:1062219/1‑80 (MQ=255)
                                                                                                                                                                  ccAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTtc                                                                                                                                                                                                                                        <  8:237245/110‑1 (MQ=255)
                                                                                                                                                                     gCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTtc                                                                                                                                                                                                                                        <  7:1032017/107‑1 (MQ=255)
                                                                                                                                                                        gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGcc                                                              <  2:213646/274‑1 (MQ=255)
                                                                                                                                                                          tGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                                                                                                                                                       <  5:130956/119‑1 (MQ=255)
                                                                                                                                                                          tGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                                                                                                                                                       >  6:130956/1‑119 (MQ=255)
                                                                                                                                                                          tGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg                                                                                                                                                                                          <  5:148716/148‑1 (MQ=255)
                                                                                                                                                                           gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACatcatc                                                                                                                                                                                                                             <  7:652824/112‑1 (MQ=255)
                                                                                                                                                                           gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg                                                                                                                                                                                         >  7:384259/1‑148 (MQ=255)
                                                                                                                                                                            aaaCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt                                                                                                                                                                                        <  8:189745/148‑1 (MQ=255)
                                                                                                                                                                            aaaCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt                                                                                                                                                                                        <  8:1043743/148‑1 (MQ=255)
                                                                                                                                                                             aaCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt                                                                                                                                                                                       <  6:652683/148‑1 (MQ=255)
                                                                                                                                                                              ccGTGCTACCTATTTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCgg                                                       <  2:6357/274‑1 (MQ=255)
                                                                                                                                                                                  gCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAttt                                                                                                                                                                                                                                         <  6:1005054/93‑1 (MQ=255)
                                                                                                                                                                                   cTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACa                                                                                                                                                                                 <  8:384259/148‑1 (MQ=255)
                                                                                                                                                                                       cGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                                                                                                                                                     >  8:1089308/1‑92 (MQ=255)
                                                                                                                                                                                       cGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                                                                                                                                                       >  7:1137885/1‑106 (MQ=255)
                                                                                                                                                                                       cGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                                                                                                                                                       <  8:1137885/106‑1 (MQ=255)
                                                                                                                                                                                           gTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTaa                                          >  2:203810/1‑275 (MQ=255)
                                                                                                                                                                                               ataagagagagaCGCGTAACCAAAGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGt                                       <  1:3610/263‑1 (MQ=255)
                                                                                                                                                                                               aCTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                                                                                                                                                                                         >  6:527349/1‑96 (MQ=255)
                                                                                                                                                                                               aCTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                                                                                                                                                                                         <  5:527349/96‑1 (MQ=255)
                                                                                                                                                                                                cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAttt                                                                                                                                                                                                                                                         >  8:580590/1‑63 (MQ=255)
                                                                                                                                                                                                cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                                                                                                                                                     >  6:691208/1‑83 (MQ=255)
                                                                                                                                                                                                cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                                                                                                                                                     <  5:691208/83‑1 (MQ=255)
                                                                                                                                                                                                    cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATg                                                                                                                                                                               <  6:932316/133‑1 (MQ=255)
                                                                                                                                                                                                    cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATg                                                                                                                                                                               >  5:932316/1‑133 (MQ=255)
                                                                                                                                                                                                      tCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGa                                                                                                                                                                                                                                            >  6:129222/1‑70 (MQ=255)
                                                                                                                                                                                                      tCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGa                                                                                                                                                                                                                                            <  5:129222/70‑1 (MQ=255)
                                                                                                                                                                                                       ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                                                                                                                                                       >  7:201205/1‑58 (MQ=255)
                                                                                                                                                                                                       ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAtt                                                                                                                                                                                                                                          <  7:1107738/71‑1 (MQ=255)
                                                                                                                                                                                                        cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTg                                                                                                                                                                                                           <  6:504991/101‑1 (MQ=255)
                                                                                                                                                                                                        cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTg                                                                                                                                                                                                           >  5:504991/1‑101 (MQ=255)
                                                                                                                                                                                                        cAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGattatt                             >  4:21353/1‑275 (MQ=255)
                                                                                                                                                                                                         aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAAc                                                                                                                                                                                  <  8:552522/125‑1 (MQ=255)
                                                                                                                                                                                                         aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAAc                                                                                                                                                                                  >  7:552522/1‑125 (MQ=255)
                                                                                                                                                                                                          aCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCa                           >  2:74449/1‑275 (MQ=255)
                                                                                                                                                                                                                 aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTgg                                                                                                                                                                                                          <  8:1011572/93‑1 (MQ=255)
                                                                                                                                                                                                                 aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATa                                                                                                                                                                   >  5:1011743/1‑132 (MQ=255)
                                                                                                                                                                                                                 aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATAATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATa                                                                                                                                                                   <  6:1011743/132‑1 (MQ=255)
                                                                                                                                                                                                                    aaCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCGCCACCGGATTTCTCg                                                                                                                                                                                                                                     <  6:825968/63‑1 (MQ=255)
                                                                                                                                                                                                                    aaCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                                                                                                                                                     >  5:825968/1‑63 (MQ=255)
                                                                                                                                                                                                                      cGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATCACCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAACACGGTTGTCAGTTTTCTTCAACACACCATGGAGGCTGTTTTGAGTAACGTTGTGGATATTCCTGGTAAAGTGCg               >  4:323791/1‑275 (MQ=255)
                                                                                                                                                                                                                          cgccgcAGCGTATCTACCGGCTGCAGCTCGGTAAGTTCCTCACCACCGGATTTCTCGACATCATTCCACGAAGTCAGCCGCTGGAGTTGCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGAGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTa                                                                         <  2:327839/213‑1 (MQ=255)
                                                                                                                                                                                                                            ctgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAg                                                                                                                                                                                           <  7:1063576/95‑1 (MQ=39)
                                                                                                                                                                                                                             cgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGaa                                                                                                                                                                                            <  7:273071/95‑1 (MQ=255)
                                                                                                                                                                                                                             cgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGaa                                                                                                                                                                                            >  8:273071/1‑95 (MQ=255)
                                                                                                                                                                                                                                    tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATa                                                                                                                                                                   >  6:1025601/1‑113 (MQ=255)
                                                                                                                                                                                                                                    tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATa                                                                                                                                                                   <  5:1025601/113‑1 (MQ=255)
                                                                                                                                                                                                                                    tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTTGAACAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTTGCGGATGGTTTGCGTAACGTCCTGATTATTCTTGGTAAAGGGCGGGATGATAAACCg  >  4:328045/1‑274 (MQ=255)
                                                                                                                                                                                                                                     aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAatat                                                                                                                                               >  8:944754/1‑132 (MQ=255)
                                                                                                                                                                                                                                       cGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTg                                                                                                                                                              <  7:915342/115‑1 (MQ=255)
                                                                                                                                                                                                                                       cGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTg                                                                                                                                                              >  8:915342/1‑115 (MQ=255)
                                                                                                                                                                                                                                       cGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATgg                                                                                     >  4:366675/1‑188 (MQ=255)
                                                                                                                                                                                                                                        gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCt                                                                                                                              >  7:1095992/1‑146 (MQ=255)
                                                                                                                                                                                                                                        gACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCt                                                                                                                              <  8:1095992/146‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                           
TTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCG  >  NC_000913/1405864‑1406365

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