Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,406,093 A→G 5.1% I39V (ATC→GTC)  smrA → DNA endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,406,0930AG 5.1% 477.6 / 10.6 139I39V (ATC→GTC) smrADNA endonuclease
Reads supporting (aligned to +/- strand):  ref base A (68/63);  new base G (4/3);  total (72/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-02

Only 100 of 146 total aligned reads displayed.

TTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCG  >  NC_000913/1405864‑1406365
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                
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)
                                                         tgttgGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCNACCGATGTCCACTGGCATCCAACNCGTANCCAACGTGCGCCGCAGCGTACCGACCCGCTGCAGCTTGATAATTTCCTcaccc                                                                                                                                                                                                                                                 >  4:540164/1‑205 (MQ=255)
                                                                 cgaccCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                                                                                                                                                    <  2:375816/275‑1 (MQ=255)
                                                                   accCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTa                                                                                                                                                                                                                                                                                  >  3:94469/1‑163 (MQ=255)
                                                                   accCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGTGTGCAACCTGGGGTGCTGGATaa                                                                                                                                                                  >  2:124105/1‑275 (MQ=255)
                                                                   accCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATaa                                                                                                                                                                  >  2:272267/1‑275 (MQ=255)
                                                                     ccGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGc                                                                                                                                                                >  4:213645/1‑275 (MQ=255)
                                                                       ggTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGa                                                                                                                                                                                                                                                                              <  1:94469/163‑1 (MQ=255)
                                                                          cgggttttttgttCGCGTTTTTAATTGGNGNGGGATGTAAAACGAACCGTGACGACAAATCGCTGGTTCTTGCCGCNATGGAAGATGTCCAGCCGCTGAAACGTTCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGNCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcgc                                                                                                                                                           <  3:41667/263‑1 (MQ=255)
                                                                          cTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcg                                                                                                                                                            <  2:518726/274‑1 (MQ=255)
                                                                                      tCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGa                                                                                                                                                                                                                                                                              <  5:742037/148‑1 (MQ=255)
                                                                                        gcgcTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGaca                                                                                                                                                                                                                                                                            <  5:824659/148‑1 (MQ=255)
                                                                                        gcgcTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGaca                                                                                                                                                                                                                                                                            >  5:925313/1‑148 (MQ=255)
                                                                                         cgcTCTTAATTGCGGGCGGAAGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCAATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCAACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGCATCGACACGCTGCAGACTGACAATTTCCTCACCCACTGATTTCTCGACCTCATCCCACTACGCAAGCCGCTGGAGTTTCGGCTGGAAGGGTTGCAACATGGGGTGCGGGATAAGCGCCGCAGACGTaaacaccca                                                                                                                                            >  4:6357/1‑269 (MQ=255)
                                                                                                     cGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGa                                                                                                                                                                                                                                                               >  6:981765/1‑148 (MQ=255)
                                                                                                      gggAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGc                                                                                                                                >  2:432444/1‑274 (MQ=255)
                                                                                                       ggaggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCtt                                                                                                                                                                                                                                                                 >  7:1055243/1‑144 (MQ=255)
                                                                                                            aTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTc                                                                                                                                                                                                                                                        <  5:981765/148‑1 (MQ=255)
                                                                                                               tAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTc                                                                                                                                                                                                                                                     <  6:925313/148‑1 (MQ=255)
                                                                                                                    tGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGc                                                                                                                                                                                                                                                                      >  7:968024/1‑126 (MQ=255)
                                                                                                                    tGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                                                                                                                                                                                                    >  7:406862/1‑144 (MQ=255)
                                                                                                                    tGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                                                                                                                                                                                                    <  8:406862/144‑1 (MQ=255)
                                                                                                                     gAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCCGCGCAGTGGTAAATATCCGCAACAGGCAGGCCTGAATCTTTTGCGcc                                                                                                                >  4:93500/1‑275 (MQ=255)
                                                                                                                      aaCCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                                                                                                                                                                                                    >  5:16231/1‑142 (MQ=255)
                                                                                                                      aaCCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                                                                                                                                                                                                    <  6:16231/142‑1 (MQ=255)
                                                                                                                         ctggacgaaCAATCCCTGTTTCTTGCGGCCCAGGAAGTTTTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGc                                                                                                             <  4:272266/271‑1 (MQ=255)
                                                                                                                         cTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTa                                                                                                                                                                                                                                                                                  <  7:1471/109‑1 (MQ=255)
                                                                                                                         cTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTa                                                                                                                                                                                                                                                                                  >  8:1471/1‑109 (MQ=255)
                                                                                                                          ttGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                                                                                                                                                                                                  <  2:563960/124‑1 (MQ=255)
                                                                                                                           tGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTgcgc                                                                                                                 >  2:485124/1‑268 (MQ=255)
                                                                                                                             acgacAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTtct                                                                                                                                                                                                                                       <  6:111614/148‑1 (MQ=255)
                                                                                                                                   aaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                                                                                                                                                                                                    <  5:828097/129‑1 (MQ=255)
                                                                                                                                   aaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTca                                                                                                                                                                                                                                                    >  6:828097/1‑129 (MQ=255)
                                                                                                                                      cGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatcat                                                                                                                                                                                                                              >  7:1043743/1‑148 (MQ=255)
                                                                                                                                           tttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGaca                                                                                                                                                                                                                                                                            <  5:809652/97‑1 (MQ=255)
                                                                                                                                            ttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                                                                                                                                                                                                   >  5:192866/1‑105 (MQ=255)
                                                                                                                                                tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTtctc                                                                                                                                                                                                                                      <  6:758434/130‑1 (MQ=255)
                                                                                                                                                          gAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                                                                                                                                                                                           <  6:396426/131‑1 (MQ=255)
                                                                                                                                                           aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCt                                                                                                                                                                                                                                                      <  8:203921/103‑1 (MQ=255)
                                                                                                                                                           aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCt                                                                                                                                                                                                                                                      >  7:203921/1‑103 (MQ=255)
                                                                                                                                                           aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGt                                                                                                                                                                                                                    <  5:233830/137‑1 (MQ=255)
                                                                                                                                                           aaGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGt                                                                                                                                                                                                                    >  6:233830/1‑137 (MQ=255)
                                                                                                                                                             gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCt                                                                                                                                                                                                                                                                        >  8:282589/1‑83 (MQ=255)
                                                                                                                                                             gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCt                                                                                                                                                                                                                                                                        <  7:282589/83‑1 (MQ=255)
                                                                                                                                                             gATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                                                                                                                                                                                           <  7:299164/128‑1 (MQ=255)
                                                                                                                                                                gTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCTGGAAGGGTTGCAACATTGGGTTCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCTAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTTTTTAGTTTTTTTcc                                                                     >  4:429715/1‑274 (MQ=255)
                                                                                                                                                                gTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTca                                                                     >  2:445368/1‑275 (MQ=255)
                                                                                                                                                                  ccAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGc                                                                                                                                                                                                                                                                      <  8:1062219/80‑1 (MQ=255)
                                                                                                                                                                  ccAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcgg                                                                                                                                                                                                  >  6:148716/1‑148 (MQ=255)
                                                                                                                                                                     gCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTtc                                                                                                                                                                                                                                        >  8:1032017/1‑107 (MQ=255)
                                                                                                                                                                        gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTa                                                                                                                                                                                                                                                                                  <  7:475093/62‑1 (MQ=255)
                                                                                                                                                                        gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGcc                                                              <  2:213646/274‑1 (MQ=255)
                                                                                                                                                                          tGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                                                                                                                                                       <  5:130956/119‑1 (MQ=255)
                                                                                                                                                                           gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACatcatc                                                                                                                                                                                                                             >  8:652824/1‑112 (MQ=255)
                                                                                                                                                                            aaaCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt                                                                                                                                                                                        <  8:1043743/148‑1 (MQ=255)
                                                                                                                                                                              aCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGa                                                                                                                                                                                                                                                                              <  8:526115/60‑1 (MQ=255)
                                                                                                                                                                                  gCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAttt                                                                                                                                                                                                                                         <  6:1005054/93‑1 (MQ=255)
                                                                                                                                                                                   cTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACa                                                                                                                                                                                 <  8:384259/148‑1 (MQ=255)
                                                                                                                                                                                       cGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                                                                                                                                                     <  7:1089308/92‑1 (MQ=255)
                                                                                                                                                                                       cGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                                                                                                                                                       <  8:1137885/106‑1 (MQ=255)
                                                                                                                                                                                       cGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                                                                                                                                                       >  7:1137885/1‑106 (MQ=255)
                                                                                                                                                                                        caTGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCg                                             <  4:445367/274‑1 (MQ=255)
                                                                                                                                                                                               ccTGGCACGCAACGCAGAACCAACGGGCGCCGACGCTTATCGACACGCTGCAGCTTGATAATTTCCTCAACACCGGANTCCTGGACATCATCCCACGAAGTGAGCCGCTGGAGTGTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTg                                      <  3:117349/274‑1 (MQ=255)
                                                                                                                                                                                               ataagagagagaCGCGTAACCAAAGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGt                                       <  1:3610/263‑1 (MQ=255)
                                                                                                                                                                                               aCTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                                                                                                                                                                                         <  5:527349/96‑1 (MQ=255)
                                                                                                                                                                                               aCTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                                                                                                                                                                                         >  6:527349/1‑96 (MQ=255)
                                                                                                                                                                                                cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAttt                                                                                                                                                                                                                                                         <  7:580590/63‑1 (MQ=255)
                                                                                                                                                                                                cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                                                                                                                                                     <  5:691208/83‑1 (MQ=255)
                                                                                                                                                                                                cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                                                                                                                                                     >  6:691208/1‑83 (MQ=255)
                                                                                                                                                                                                cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGACTTCTCGACATCATCCCACTAAGTCAGCCGCTGGATTTTCTGCGGGAAGGTTTGCAACCTTGGGTGCCGGATAACCTGCGTAGTGGTAACTATCCGCAACAGGCGAGCCTGAACATTTTCCGCCAGCCGGTGGAAGAGTTCCGCCAAATGGTGTTCCGTTTTATTCACAAAGCCATGGCGGCTGGCTCGCGTacagtgg                                     >  4:335335/1‑269 (MQ=255)
                                                                                                                                                                                                      tCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGa                                                                                                                                                                                                                                            >  6:129222/1‑70 (MQ=255)
                                                                                                                                                                                                      tCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGa                                                                                                                                                                                                                                            <  5:129222/70‑1 (MQ=255)
                                                                                                                                                                                                       ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                                                                                                                                                       <  8:201205/58‑1 (MQ=255)
                                                                                                                                                                                                       ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                                                                                                                                                       >  7:201205/1‑58 (MQ=255)
                                                                                                                                                                                                       ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAtt                                                                                                                                                                                                                                          >  8:1107738/1‑71 (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                                                                                                                                                                                                          >  7:1011572/1‑93 (MQ=255)
                                                                                                                                                                                                                 aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATa                                                                                                                                                                   >  5:1011743/1‑132 (MQ=255)
                                                                                                                                                                                                                    aaCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCGCCACCGGATTTCTCg                                                                                                                                                                                                                                     <  6:825968/63‑1 (MQ=255)
                                                                                                                                                                                                                      cGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATCACCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAACACGGTTGTCAGTTTTCTTCAACACACCATGGAGGCTGTTTTGAGTAACGTTGTGGATATTCCTGGTAAAGTGCg               >  4:323791/1‑275 (MQ=255)
                                                                                                                                                                                                                      cGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGCATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATCCACAAAGCCCTGGCGGCTGGTTTGCGTAACGTGCTGATTATTCATGTTTAAg                   >  2:382030/1‑271 (MQ=255)
                                                                                                                                                                                                                       ttGCGCCGCGGCGTATCCACACGGTCCAGCTTGAATATTTCCTCACCACCGGTTTTCGCGACATCATCCCACTAAGTCAGCCGCAGGAGTGTCGGCGGGAAGTGTGGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCg               <  4:382029/273‑1 (MQ=255)
                                                                                                                                                                                                                          cgccgcAGCGTATCTACCGGCTGCAGCTCGGTAAGTTCCTCACCACCGGATTTCTCGACATCATTCCACGAAGTCAGCCGCTGGAGTTGCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGAGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTa                                                                         <  2:327839/213‑1 (MQ=255)
                                                                                                                                                                                                                            ctgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAg                                                                                                                                                                                           >  8:1063576/3‑97 (MQ=39)
                                                                                                                                                                                                                            ctgcAGCGTGTCGATACGCTGCGGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAg                                                                                                                                                                                           <  7:1063576/95‑1 (MQ=39)
                                                                                                                                                                                                                             cgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGaa                                                                                                                                                                                            <  7:273071/95‑1 (MQ=255)
                                                                                                                                                                                                                             cgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGaa                                                                                                                                                                                            >  8:273071/1‑95 (MQ=255)
                                                                                                                                                                                                                                 gCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAAt    <  2:21353/275‑1 (MQ=255)
                                                                                                                                                                                                                                  cGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGa                                                                                                                                  >  6:560684/1‑148 (MQ=255)
                                                                                                                                                                                                                                    tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATa                                                                                                                                                                   >  6:1025601/1‑113 (MQ=255)
                                                                                                                                                                                                                                    tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATa                                                                                                                                                                   <  5:1025601/113‑1 (MQ=255)
                                                                                                                                                                                                                                     aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAatat                                                                                                                                               <  7:944754/132‑1 (MQ=255)
                                                                                                                                                                                                                                     aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAatat                                                                                                                                               >  8:944754/1‑132 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                
TTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCG  >  NC_000913/1405864‑1406365

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