Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,020,511 C→G 5.4% intergenic (+22/‑29) ssnA → / → ygfM putative chlorohydrolase/aminohydrolase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,5110CG 5.4% 683.0 / 20.9 195intergenic (+22/‑29)ssnA/ygfMputative chlorohydrolase/aminohydrolase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (83/101);  new base G (7/4);  total (90/105)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

CTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCG  >  NC_000913/3020244‑3020773
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                      
cTTTTGGGAGCAGGTTACAAAGCGGATGTAACCATGGGCGTGTACAGCTAGCCGGGGCCGATGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGGTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGaa                                                                                                                                                                                                                                                                  <  4:288903/274‑1 (MQ=255)
           cGGTTACAAAGTTGATTTAACCATTGGGGATTACAACTCGCCGAGGCCGGTGCGGGCAGACAAGATCGCCGGGCATATCGCGTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATa                                                                                                                                                                                                                                                      <  4:296814/275‑1 (MQ=255)
               ttaaaagcgcgtttACCCATTTGGGATTACAACTGCCGGACCCGCCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGGGAGGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGa                                                                                                                                                                                                                                                  <  4:148795/264‑1 (MQ=255)
                                   ttGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTGGATTGCGATTACATTTATGCACAAGCCAGAAAAGCAGCTGCGAGAATGTGGCGTCGGATGTAGGCGCTGGCATAAAGGACAGATGCCATCTTCCCGCAGGGAAGAGGGTTAAATACGGCAAGAAGATGATTGAACAAtt                                                                                                                                                                                                                              >  1:14848/1‑275 (MQ=255)
                                                                                    cATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGCTGCCCTCTTCCCGCAGGGAAGAGGGTTAAATACGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGCCGCGGAAAcgtaaactccccca                                                                                                                                                                              >  4:244922/1‑261 (MQ=255)
                                                                                    cATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAg                                                                                                                                                                                      >  4:282290/1‑266 (MQ=255)
                                                                                                                 cGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGAAGCTTGAACAATTTTTCAGGCCAGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGACACCGTCTGGTTCGCCtggg                                                                                                                                                 >  2:275563/1‑274 (MQ=255)
                                                                                                                         gCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAg                                                                                                                                                                                                                                                                      <  6:1092767/149‑1 (MQ=255)
                                                                                                                         gCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAg                                                                                                                                                                                                                                                                      >  5:964636/1‑149 (MQ=255)
                                                                                                                         gCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCCCTGGCATAAATGACAGATGCCCTCTTCCCGCAg                                                                                                                                                                                                                                                                      <  7:1240407/149‑1 (MQ=255)
                                                                                                                               tGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAaga                                                                                                                                                                                                                                                                <  6:978463/149‑1 (MQ=255)
                                                                                                                               tGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAaga                                                                                                                                                                                                                                                                <  8:655655/149‑1 (MQ=255)
                                                                                                                                ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAagag                                                                                                                                                                                                                                                               <  5:733778/149‑1 (MQ=255)
                                                                                                                                ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAagag                                                                                                                                                                                                                                                               <  7:1345788/149‑1 (MQ=255)
                                                                                                                                ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATGCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAagag                                                                                                                                                                                                                                                               <  5:355645/149‑1 (MQ=255)
                                                                                                                                ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGCTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCGGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCaa                                                                                                                                    <  2:30937/272‑1 (MQ=255)
                                                                                                                                 gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgg                                                                                                                                                                                                                                                              <  8:706799/149‑1 (MQ=255)
                                                                                                                                 gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgg                                                                                                                                                                                                                                                              <  7:557343/149‑1 (MQ=255)
                                                                                                                                 gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgg                                                                                                                                                                                                                                                              <  6:926778/149‑1 (MQ=255)
                                                                                                                                 gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAgg                                                                                                                                                                                                                                                              <  5:923557/149‑1 (MQ=255)
                                                                                                                                 gTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAgg                                                                                                                                                                                                                                                              <  5:1434866/149‑1 (MQ=255)
                                                                                                                                    aaTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGtt                                                                                                                                                                                                                                                           >  7:412952/1‑149 (MQ=255)
                                                                                                                                    aaTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGtt                                                                                                                                                                                                                                                           >  5:1014851/1‑149 (MQ=255)
                                                                                                                                        gtgtGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAt                                                                                                                                                                                                                                                       >  5:1234375/1‑149 (MQ=255)
                                                                                                                                        gtgtGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACt                                                                                                                                       <  2:286951/261‑1 (MQ=255)
                                                                                                                                           tgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAg                                                                                                                                                                                                                                                    <  7:701095/149‑1 (MQ=255)
                                                                                                                                           tgATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAg                                                                                                                                                                                                                                                    <  8:1430183/149‑1 (MQ=255)
                                                                                                                                              tgtgtgaTGAAGACCGGCAGTGTAACTTCGATTGCGCTTCCATTTATCCACAAGCCCGAAAAGGCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGCTGCCTTCTTCCCGCAGGGAAGATGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTa                                                                                                                   <  4:167351/269‑1 (MQ=255)
                                                                                                                                                  cTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaa                                                                                                                                                                                                                                             <  6:669625/149‑1 (MQ=255)
                                                                                                                                                  cTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaa                                                                                                                                                                                                                                             <  5:281626/149‑1 (MQ=255)
                                                                                                                                                  cTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaa                                                                                                                                                                                                                                             <  7:1080579/149‑1 (MQ=255)
                                                                                                                                                   tATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaag                                                                                                                                                                                                                                            >  6:1412073/1‑149 (MQ=255)
                                                                                                                                                   tATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAt                                                                                                                                                                                                                               >  2:226814/1‑162 (MQ=255)
                                                                                                                                                    aTGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAagaaga                                                                                                                                                                                                                                           <  6:711174/149‑1 (MQ=255)
                                                                                                                                                     tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                                                                                                                                          >  7:912878/1‑149 (MQ=255)
                                                                                                                                                     tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                                                                                                                                          >  6:1347344/1‑149 (MQ=255)
                                                                                                                                                      gAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatg                                                                                                                                                                                                                                         >  8:209444/1‑149 (MQ=255)
                                                                                                                                                         gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                                                                                                                                      <  6:901657/149‑1 (MQ=255)
                                                                                                                                                         gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                                                                                                                                      <  7:122197/149‑1 (MQ=255)
                                                                                                                                                         gACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                                                                                                                                      <  5:311383/149‑1 (MQ=255)
                                                                                                                                                           ccGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAgg                                                                                                                                                                                                                                                                     >  2:47606/1‑116 (MQ=255)
                                                                                                                                                             gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgg                                                                                                                                                                                                                                                              >  7:792094/1‑121 (MQ=255)
                                                                                                                                                             gTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgg                                                                                                                                                                                                                                                              <  8:792094/121‑1 (MQ=255)
                                                                                                                                                               cAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAag                                                                                                                                                                                                                                                                 <  4:47606/116‑1 (MQ=255)
                                                                                                                                                               cAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACaa                                                                                                                                                                                                                                <  8:912878/149‑1 (MQ=255)
                                                                                                                                                                  tttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGTATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGCAGAGGGTTAAATAAGGAAAGAAGATGACTGAACAAttt                                                                                                                                                                                                                             >  8:1271064/1‑149 (MQ=255)
                                                                                                                                                                  tttAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgg                                                                                                                                                                                                                                                              >  6:231771/1‑116 (MQ=255)
                                                                                                                                                                      aCTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                                                                                                                                         <  8:1381861/149‑1 (MQ=255)
                                                                                                                                                                            cttGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTCCATTGCGTCGGATGGATGCCCTGGCATAAATGACAGGTGCCCCCTTCCCGCCGGGTAGAGGACTAAATCAGGAAAGTAGATGATTGAACAATTTTTCATGACATACGCTGTCTAACAGGCGCAGGAACTGAAACGCCGCTCCCCAGATGCAGCCGACTGCTGCGCCCGGGGAGGCAGACTCAACGCTACCCCAACCAGTACGGCTAACAAGCTTGACAttcagtcaaac                                                                                      >  4:80065/2‑265 (MQ=255)
                                                                                                                                                                                 gATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGActct                                                                                                                                                                                                              >  8:1231673/1‑149 (MQ=255)
                                                                                                                                                                                       aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                                                                                                                                        <  6:431610/149‑1 (MQ=255)
                                                                                                                                                                                       aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATACCCTCTTCCCTGCGGGAAGAggg                                                                                                                                                                                                                                                             <  5:1240612/96‑1 (MQ=255)
                                                                                                                                                                                       aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATACCCTCTTCCCTGCGGGAAGAggg                                                                                                                                                                                                                                                             >  6:1240612/1‑96 (MQ=255)
                                                                                                                                                                                        tttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCa                                                                                                                                                                                                                         <  5:651145/131‑1 (MQ=255)
                                                                                                                                                                                         ttATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACa                                                                                                                                                                                                      <  7:209444/149‑1 (MQ=255)
                                                                                                                                                                                              cacaAGTCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCGATCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGgcgc                                                                                                                                                                                                 <  8:1311063/149‑1 (MQ=255)
                                                                                                                                                                                               acaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCt                                                                                                                                                                                                >  7:645973/1‑149 (MQ=255)
                                                                                                                                                                                                caAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTg                                                                                                                                                                                               <  6:510382/149‑1 (MQ=255)
                                                                                                                                                                                                caAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTg                                                                                                                                                                                               >  5:1191352/1‑149 (MQ=255)
                                                                                                                                                                                                   gCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGATGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGaa                                                                                                                                                                                            >  8:1050279/1‑149 (MQ=255)
                                                                                                                                                                                                   gCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGAc                                                              <  2:113112/275‑1 (MQ=255)
                                                                                                                                                                                                     cAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACt                                                                                                                                                                                          <  7:878414/149‑1 (MQ=255)
                                                                                                                                                                                                     cAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACt                                                                                                                                                                                          >  5:1326741/1‑149 (MQ=255)
                                                                                                                                                                                                      aGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGACGCCGTCTGGTTCGCCGGGCGCAGCAAACTCAATGCTACACCAACCCGTACCGATACAAAGATTGCCATTTCCTTACAGGATCTGGACCTGGACTGGGTTgcccg                                                            >  4:254082/1‑270 (MQ=255)
                                                                                                                                                                                                        aaaaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGaa                                                                                                                                                                                       <  7:1440293/149‑1 (MQ=255)
                                                                                                                                                                                                                  gCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTAc                                                                                                                                                                             >  6:7277/1‑149 (MQ=255)
                                                                                                                                                                                                                     aGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTc                                                                                                                                                                                                                          <  8:527002/101‑1 (MQ=255)
                                                                                                                                                                                                                     aGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTc                                                                                                                                                                                                                          >  7:527002/1‑101 (MQ=255)
                                                                                                                                                                                                                      gTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAgg                                                                                                                                                                         >  8:1102885/1‑149 (MQ=255)
                                                                                                                                                                                                                       tATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGa                                                                                                                                                                        <  7:391406/149‑1 (MQ=255)
                                                                                                                                                                                                                        aTGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGAt                                                                                                                                                                       <  5:565979/149‑1 (MQ=255)
                                                                                                                                                                                                                               gTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAAAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAAGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTgca                                   >  2:91894/1‑272 (MQ=255)
                                                                                                                                                                                                                                           cgcTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCg                                                                                                                                                    <  5:1357573/149‑1 (MQ=255)
                                                                                                                                                                                                                                              tGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCgggg                                                                                                                                                 >  6:1456570/1‑149 (MQ=255)
                                                                                                                                                                                                                                                gCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgcc                                                                                                                                                                                  <  5:801502/114‑1 (MQ=255)
                                                                                                                                                                                                                                                gCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgcc                                                                                                                                                                                  >  6:801502/1‑114 (MQ=255)
                                                                                                                                                                                                                                                   tAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGAtt                                                                                                                                                                                                                                      >  5:429608/1‑59 (MQ=255)
                                                                                                                                                                                                                                                   tAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcagc                                                                                                                                            >  5:1095533/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    aaaTGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACa                                                                                                                                                                                                      <  7:1416152/90‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaaTGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACa                                                                                                                                                                                                      >  8:1416152/1‑90 (MQ=255)
                                                                                                                                                                                                                                                    aaaTGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgcagca                                                                                                                                           >  8:1163759/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      aTGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCaaa                                                                                                                                         >  5:464731/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      aTGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTACCGATAAAAAGATTGCCATTTCCt                                                                                           >  1:35135/1‑195 (MQ=255)
                                                                                                                                                                                                                                                          cAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGCTTGGCGCAATGTCTCGCTTTCAGCCattcggt        >  4:184789/1‑268 (MQ=255)
                                                                                                                                                                                                                                                           aGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTGTGGTTCGCCGGGGGCAGCAAACTCaa                                                                                                                                    <  8:74746/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              tGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGc                                                                                                                                 <  4:107664/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               acccTCTTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCt                                                                                                                                >  7:653399/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                 ccTCTTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac                                                                                                                              >  6:607951/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 ccTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac                                                                                                                              >  8:583076/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 ccTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac                                                                                                                              >  7:758698/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  ctctTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca                                                                                                                             >  7:295597/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  ctctTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca                                                                                                                             >  7:1296616/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   tctTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTacac                                                                                                                            >  6:366418/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   tctTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTacac                                                                                                                            <  8:412952/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ttCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCa                                                                                                                          >  5:1063385/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       cccGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAc                                                                                                                        <  5:10345/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cccGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAc                                                                                                                        <  7:1028373/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ccGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAcc                                                                                                                       >  5:689143/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        ccGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAcc                                                                                                                       >  6:957508/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           gCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCg                                                                                                                     >  6:1064976/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                      
CTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCG  >  NC_000913/3020244‑3020773

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