Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,037,957 G→T 8.6% P50P (CCC→CCA 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,9570GT 8.6% 471.9 / 24.4 139P50P (CCC→CCArecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base G (61/65);  new base T (6/6);  total (68/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

AACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCA  >  NC_000913/3037809‑3038105
                                                                                                                                                    |                                                                                                                                                    
aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAg                                                                                                                                                      <  4:234180/149‑1 (MQ=255)
 aCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAgg                                                                                                                                                     <  4:359314/149‑1 (MQ=255)
 aCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAgg                                                                                                                                                     >  4:460252/1‑149 (MQ=255)
        gagcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaa                                                                                                                                              <  2:18277/149‑2 (MQ=255)
        gagcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaa                                                                                                                                              <  3:383407/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                                                                                                                                      <  4:644878/149‑1 (MQ=255)
                cGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTaa                                                                                                                                      <  1:651839/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)
                                  cGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCtt                                                                                                                    <  3:529667/149‑1 (MQ=255)
                                   gTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTg                                                                                                                   <  2:122766/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                                                                                                        >  1:229737/1‑149 (MQ=255)
                                               ccGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTAc                                                                                                       >  2:1326621/1‑149 (MQ=255)
                                                  tCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTcc                                                                                                    >  2:411244/1‑149 (MQ=255)
                                                          gAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTgg                                                                                            >  3:762894/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                                                                                                                                                  >  4:196266/1‑83 (MQ=255)
                                                                      cacaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGgca                                                                                                                                                  <  3:196266/83‑1 (MQ=255)
                                                                      cacaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGgc                                                                                >  2:1219411/1‑149 (MQ=255)
                                                                      cacaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGgc                                                                                >  2:1150933/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:1447244/149‑1 (MQ=255)
                                                                                   gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAgg                                                                   <  1:918926/149‑1 (MQ=255)
                                                                                      ccTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCaa                                                                <  2:815019/149‑1 (MQ=255)
                                                                                            cgAAAGGCGTTGTAAAGGATCTCAACGTCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGc                                                          <  4:1357085/149‑1 (MQ=255)
                                                                                            cgAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGATCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGc                                                          <  1:1300310/149‑1 (MQ=255)
                                                                                             gAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCg                                                         <  1:1000528/149‑1 (MQ=255)
                                                                                             gAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCg                                                         >  4:1191393/1‑149 (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)
                                                                                                           aGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAAGGCAGCATACCTTTAAAACTGCGTTCCAGTTCTTGCGCACTTCGTACTCCCCGGCTGGAGTATAAACGGCTCAGCAAGGGAGGAAATTCAGTGGGCAAGTCTGACGt                                           >  3:1445217/1‑149 (MQ=255)
                                                                                                            ggATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGtt                                          <  3:1191393/149‑1 (MQ=255)
                                                                                                             gATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGttt                                         <  4:81336/149‑1 (MQ=255)
                                                                                                              aTCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTc                                        <  2:398252/149‑1 (MQ=255)
                                                                                                               tctcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCa                                       >  3:1357006/1‑149 (MQ=255)
                                                                                                                   aaCGGCCTTTTTGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCAt                                      >  3:265249/1‑146 (MQ=255)
                                                                                                                   aaCGGCCTTTTTGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCAt                                      <  4:265249/146‑1 (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)
                                                                                                                               gACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcgac                       >  2:239369/1‑149 (MQ=255)
                                                                                                                                   ccgcTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAg                   <  4:43725/149‑1 (MQ=255)
                                                                                                                                     gcTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGtt                 >  4:628373/1‑149 (MQ=255)
                                                                                                                                     gcTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGtt                 >  3:193407/1‑149 (MQ=255)
                                                                                                                                      cTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg                >  2:1031842/1‑149 (MQ=255)
                                                                                                                                      cTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg                <  4:762894/149‑1 (MQ=255)
                                                                                                                                      cTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg                <  4:652555/149‑1 (MQ=255)
                                                                                                                                       tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                     <  2:1082195/127‑1 (MQ=255)
                                                                                                                                       tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt               >  3:548263/1‑149 (MQ=255)
                                                                                                                                       tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt               <  1:1164507/149‑1 (MQ=255)
                                                                                                                                       tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt               >  1:519371/1‑149 (MQ=255)
                                                                                                                                       aggtatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt               >  4:701719/6‑149 (MQ=255)
                                                                                                                                       aggtatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt               >  2:1320444/6‑149 (MQ=255)
                                                                                                                                        ggtatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTa              >  3:463455/5‑149 (MQ=255)
                                                                                                                                        ggtatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTa              >  3:122786/5‑149 (MQ=255)
                                                                                                                                         aGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTAt             >  3:592785/1‑149 (MQ=255)
                                                                                                                                           ttGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCt           >  4:220926/1‑149 (MQ=255)
                                                                                                                                           ttGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCt           <  4:193407/149‑1 (MQ=255)
                                                                                                                                           ttGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCt           >  1:274636/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     <  4:34361/149‑1 (MQ=255)
                                                                                                                                                 ccAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTgttgtt     >  1:41151/1‑149 (MQ=255)
                                                                                                                                                 ccAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTgttgtt     <  3:253279/149‑1 (MQ=255)
                                                                                                                                                  cAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGttt    >  3:73349/1‑149 (MQ=255)
                                                                                                                                                    gggCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCa  <  2:896389/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCA  >  NC_000913/3037809‑3038105

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