Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,037,956 A→C 8.3% W51G (TGG→GGG)  recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,037,9560AC 8.3% 477.6 / 33.3 145W51G (TGG→GGG) recJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base A (65/68);  new base C (6/6);  total (71/74)
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.88e-01

Only 100 of 146 total aligned reads displayed.

CAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTT  >  NC_000913/3037808‑3038103
                                                                                                                                                    |                                                                                                                                                   
 aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAg                                                                                                                                                    <  4:234180/149‑1 (MQ=255)
 aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAg                                                                                                                                                    >  4:1307429/1‑149 (MQ=255)
  aCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAgg                                                                                                                                                   <  4:359314/149‑1 (MQ=255)
  aCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAgg                                                                                                                                                   >  4:460252/1‑149 (MQ=255)
         gagcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaa                                                                                                                                            <  3:383407/149‑2 (MQ=255)
         gagcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaa                                                                                                                                            <  2:18277/149‑2 (MQ=255)
          agcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaac                                                                                                                                           <  2:180665/149‑3 (MQ=255)
          agcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATa                                                                                                                                           <  3:1487614/149‑1 (MQ=255)
             gcATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaactga                                                                                                                                        <  1:29321/149‑6 (MQ=255)
               aTCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCttt                                                                                                                                      <  1:430700/149‑1 (MQ=255)
                 cGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTaa                                                                                                                                    <  1:651839/149‑1 (MQ=255)
                 cGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTaa                                                                                                                                    <  4:644878/149‑1 (MQ=255)
                  gCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAac                                                                                                                                   <  2:1115732/149‑1 (MQ=255)
                    cAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAacac                                                                                                                                 >  4:520787/1‑149 (MQ=255)
                      gCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTg                                                                                                                               <  2:1441086/149‑1 (MQ=255)
                            ttAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTcc                                                                                                                         >  3:652555/1‑149 (MQ=255)
                             tAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCa                                                                                                                        >  1:1335321/1‑149 (MQ=255)
                                  ccGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCt                                                                                                                   >  2:1045719/1‑149 (MQ=255)
                                    gTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTg                                                                                                                 <  2:122766/149‑1 (MQ=255)
                                     tGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTgc                                                                                                                <  3:14040/149‑1 (MQ=255)
                                         ggTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCAc                                                                                                            <  2:736548/149‑1 (MQ=255)
                                           tCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTg                                                                                                          >  2:986200/1‑149 (MQ=255)
                                             gcgcCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCg                                                                                                        <  4:843408/149‑1 (MQ=255)
                                              cgcCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGt                                                                                                       <  4:873274/149‑1 (MQ=255)
                                               gcCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTa                                                                                                      <  3:1149116/149‑1 (MQ=255)
                                               gcCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTa                                                                                                      >  1:229737/1‑149 (MQ=255)
                                                ccGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTAc                                                                                                     >  2:1326621/1‑149 (MQ=255)
                                                   tCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTcc                                                                                                  >  2:411244/1‑149 (MQ=255)
                                                     ggCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTcccc                                                                                                >  3:223163/1‑149 (MQ=255)
                                                           gAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTgg                                                                                          >  3:762894/1‑149 (MQ=255)
                                                           gAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTgg                                                                                          >  3:586647/1‑149 (MQ=255)
                                                            aaaTCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGCCGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGc                                                                                         >  1:1291748/1‑149 (MQ=255)
                                                            aaaTCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGc                                                                                         <  1:262705/149‑1 (MQ=255)
                                                            aaaTCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGc                                                                                         >  1:257442/1‑149 (MQ=255)
                                                            aaaTCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCACTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGc                                                                                         <  2:475259/149‑1 (MQ=255)
                                                             aaTCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTaa                                                                                                                                    <  3:1401458/105‑1 (MQ=255)
                                                             aaTCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTaa                                                                                                                                    >  4:1401458/1‑105 (MQ=255)
                                                                    gACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTGAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAAc                                                                                 <  4:123800/149‑1 (MQ=255)
                                                                     aCCACAATAATGCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACg                                                                                >  1:833297/1‑149 (MQ=255)
                                                                      ccACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACgg                                                                               >  4:873498/1‑149 (MQ=255)
                                                                       cacaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGgca                                                                                                                                                <  3:196266/83‑1 (MQ=255)
                                                                       cacaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGgca                                                                                                                                                >  4:196266/1‑83 (MQ=255)
                                                                       cacaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGgc                                                                              >  2:1150933/1‑149 (MQ=255)
                                                                       cacaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGgc                                                                              >  2:1219411/1‑149 (MQ=255)
                                                                          aataatCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCgca                                                                           >  4:253279/1‑149 (MQ=255)
                                                                               tCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAg                                                                      <  2:1491199/149‑1 (MQ=255)
                                                                                 cgcgTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAggg                                                                    <  1:1325786/149‑1 (MQ=255)
                                                                                  gcgTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGa                                                                   <  2:833297/149‑1 (MQ=255)
                                                                                  gcgTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGa                                                                   <  2:1465993/149‑1 (MQ=255)
                                                                                   cgTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAg                                                                  <  4:790649/149‑1 (MQ=255)
                                                                                    gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAgg                                                                 <  1:918926/149‑1 (MQ=255)
                                                                                    gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAgg                                                                 <  1:1447244/149‑1 (MQ=255)
                                                                                       ccTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCaa                                                              <  2:815019/149‑1 (MQ=255)
                                                                                           cgcgAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCa                                                          <  4:1459743/149‑1 (MQ=255)
                                                                                             cgAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGATCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGc                                                        <  1:1300310/149‑1 (MQ=255)
                                                                                              gAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCg                                                       <  1:1000528/149‑1 (MQ=255)
                                                                                                     gTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGt                                                <  2:186853/149‑1 (MQ=255)
                                                                                                        gTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTg                                             >  1:194116/1‑149 (MQ=255)
                                                                                                         tAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGc                                            <  2:668744/149‑1 (MQ=255)
                                                                                                         tAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGc                                            >  3:472056/1‑149 (MQ=255)
                                                                                                          aaaGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGcc                                           <  4:329445/149‑1 (MQ=255)
                                                                                                              gATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGttt                                       <  4:81336/149‑1 (MQ=255)
                                                                                                               aTCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTc                                      <  2:398252/149‑1 (MQ=255)
                                                                                                                 ctcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCAt                                    <  3:52258/149‑1 (MQ=255)
                                                                                                                    aaCGGCCTTTTTGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCAt                                    >  3:265249/1‑146 (MQ=255)
                                                                                                                    aaCGGCCTTTTTGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCAt                                    <  4:265249/146‑1 (MQ=255)
                                                                                                                    aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACTTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGa                                 >  3:847963/1‑149 (MQ=255)
                                                                                                                    aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGc                                            >  2:907726/1‑138 (MQ=255)
                                                                                                                    aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGc                                            <  1:907726/138‑1 (MQ=255)
                                                                                                                    aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                   <  2:1409468/147‑1 (MQ=255)
                                                                                                                    aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGa                                 >  2:83216/1‑149 (MQ=255)
                                                                                                                     aCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGAc                                <  1:1150933/149‑1 (MQ=255)
                                                                                                                       ggCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACtt                              <  1:1045719/149‑1 (MQ=255)
                                                                                                                          cTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCAc                           <  3:1566462/149‑1 (MQ=255)
                                                                                                                                gACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcgac                     >  2:239369/1‑149 (MQ=255)
                                                                                                                                    ccgcTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAg                 <  4:43725/149‑1 (MQ=255)
                                                                                                                                      gcTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGtt               >  3:657258/1‑149 (MQ=255)
                                                                                                                                      gcTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGtt               >  4:628373/1‑149 (MQ=255)
                                                                                                                                      gcTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGtt               <  3:403018/149‑1 (MQ=255)
                                                                                                                                       cTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg              <  4:762894/149‑1 (MQ=255)
                                                                                                                                       cTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg              >  2:1031842/1‑149 (MQ=255)
                                                                                                                                       cTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg              <  4:652555/149‑1 (MQ=255)
                                                                                                                                       cTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg              >  3:1229343/1‑149 (MQ=255)
                                                                                                                                        tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                   <  2:1082195/127‑1 (MQ=255)
                                                                                                                                        tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt             >  1:519371/1‑149 (MQ=255)
                                                                                                                                        tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt             >  3:548263/1‑149 (MQ=255)
                                                                                                                                        tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt             <  1:1164507/149‑1 (MQ=255)
                                                                                                                                        aggtatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt             >  4:701719/6‑149 (MQ=255)
                                                                                                                                        aggtatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt             >  2:1320444/6‑149 (MQ=255)
                                                                                                                                         ggtatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTa            >  3:122786/5‑149 (MQ=255)
                                                                                                                                            ttGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCt         >  1:274636/1‑149 (MQ=255)
                                                                                                                                            ttGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCt         <  4:193407/149‑1 (MQ=255)
                                                                                                                                            ttGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCt         >  4:220926/1‑149 (MQ=255)
                                                                                                                                             tgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtg        >  1:874420/1‑149 (MQ=255)
                                                                                                                                              gctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtgt       <  1:1060806/149‑1 (MQ=255)
                                                                                                                                               ctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtgtt      <  4:223163/149‑1 (MQ=255)
                                                                                                                                                tgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtgttg     <  2:228154/149‑1 (MQ=255)
                                                                                                                                                 gcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtgttgt    >  4:1394304/1‑149 (MQ=255)
                                                                                                                                                  ccAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTgttgtt   >  1:41151/1‑149 (MQ=255)
                                                                                                                                                  ccAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTgttgtt   <  4:34361/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTT  >  NC_000913/3037808‑3038103

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