Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,047,208 C→T 10.0% R81R (CGC→CGT polA → 5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,047,2080CT10.0% 199.3 / 10.0 69R81R (CGC→CGTpolA5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base T (5/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GCGAGCCGACCGGTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTC  >  NC_000913/4047063‑4047356
                                                                                                                                                 |                                                                                                                                                    
gCGAGCCGACCGGTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATcgcccg                                                                                                                                                   >  2:473860/1‑149 (MQ=255)
    gCCGACCGGTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCaa                                                                                                                                               >  1:500416/1‑149 (MQ=255)
     ccGACCGGTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAAt                                                                                                                                              >  1:445081/1‑149 (MQ=255)
        aCCGGTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGcc                                                                                                                                           <  4:634871/149‑1 (MQ=255)
               cGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGTCCGGCATTGGCGGGCGAt                                                                                                                                    <  1:581217/149‑1 (MQ=255)
                   gTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGc                                                                                                                                <  3:512075/149‑1 (MQ=255)
                     aTGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGGACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGt                                                                                                                              <  2:524902/149‑1 (MQ=255)
                     aTGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGt                                                                                                                              >  4:558000/1‑149 (MQ=255)
                     aTGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGt                                                                                                                              <  2:592621/149‑1 (MQ=255)
                        gtgtCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGca                                                                                                                           >  4:687443/1‑149 (MQ=255)
                         tgtCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGcac                                                                                                                          <  1:219481/149‑1 (MQ=255)
                             cTCAACATGTTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                      >  2:417472/1‑149 (MQ=255)
                             cTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                      <  1:488553/149‑1 (MQ=255)
                             cTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                      <  2:445081/149‑1 (MQ=255)
                             cTCAACATGCTGCGCAGTCTGATAATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAAt                                                                                                                      <  2:274679/149‑1 (MQ=255)
                                aaCATGCTGCGCAGACTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGa                                                                                                                   >  3:48620/1‑149 (MQ=255)
                                  cATGCTGCGCAGTCTGAAAATGCAATATAAACCGACGCATGCAGCGGTGGTGTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACTAATCACATCGCCGTCCAATGCCGGACGATCTGCGTGCACAAATCGAAc                                                                                                                 <  4:265101/149‑1 (MQ=255)
                                     gctgcGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGGCGATCTGCGTGCACAAATCg                                                                                                                    >  2:7306/1‑143 (MQ=255)
                                     gctgcGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCg                                                                                                                    <  1:7306/143‑1 (MQ=255)
                                                 gATCATGCAATATAATCCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATgg                                                                                                  <  3:160572/149‑1 (MQ=255)
                                                  aTCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGt                                                                                                 >  2:640787/1‑149 (MQ=255)
                                                  aTCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGt                                                                                                 <  1:417472/149‑1 (MQ=255)
                                                   tCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGtt                                                                                                <  3:212187/149‑1 (MQ=255)
                                                      tGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTaaa                                                                                             <  2:583831/149‑1 (MQ=255)
                                                      tGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGAAGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTaaa                                                                                             <  1:640787/149‑1 (MQ=255)
                                                         aaTATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCg                                                                                          >  4:623880/1‑149 (MQ=255)
                                                           tataAACCGACGCATGCAGCGGCGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGCTAAAGCGAt                                                                                        <  1:594439/149‑1 (MQ=255)
                                                             taAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATgg                                                                                      >  2:362962/1‑149 (MQ=255)
                                                                aCCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGAc                                                                                   >  2:281283/1‑149 (MQ=255)
                                                                        aTGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCgctg                                                                           >  4:672127/1‑149 (MQ=255)
                                                                         tGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCgctgc                                                                          <  3:594766/149‑1 (MQ=255)
                                                                         tGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCgctgc                                                                          <  3:520048/149‑1 (MQ=255)
                                                                               ggtggtCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCgg                                                                    <  1:473860/149‑1 (MQ=255)
                                                                                   gtCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTc                                                                >  4:665042/1‑149 (MQ=255)
                                                                                    tCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCt                                                               >  4:216752/1‑149 (MQ=255)
                                                                                        tGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCg                                                           <  3:248875/149‑1 (MQ=255)
                                                                                                    aaaaaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTaa                                                                                              <  4:463100/102‑1 (MQ=255)
                                                                                                    aaaaaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTaa                                                                                              >  3:463100/1‑102 (MQ=255)
                                                                                                      aaaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGgacgac                                             <  4:463766/149‑1 (MQ=255)
                                                                                                       aaCCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGacgacg                                            >  2:620375/1‑149 (MQ=255)
                                                                                                         ccTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGtt                                          <  3:542728/149‑1 (MQ=255)
                                                                                                         ccTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGtt                                          <  3:687443/149‑1 (MQ=255)
                                                                                                           tttCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTAt                                        <  3:216752/149‑1 (MQ=255)
                                                                                                           tttCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTAt                                        >  3:354508/1‑149 (MQ=255)
                                                                                                             tCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCg                                      <  3:672127/149‑1 (MQ=255)
                                                                                                                tgatgaACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTa                                   >  4:151525/1‑149 (MQ=255)
                                                                                                                   tgaACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTActc                                >  1:354794/1‑149 (MQ=255)
                                                                                                                    gaACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTActct                               >  4:261763/1‑149 (MQ=255)
                                                                                                                       cTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGgc                            >  1:348089/1‑149 (MQ=255)
                                                                                                                        tGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGgcg                           <  3:452290/149‑1 (MQ=255)
                                                                                                                          tttGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGgcgcg                         >  1:689174/1‑149 (MQ=255)
                                                                                                                           ttGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGgcgcgc                        >  2:8993/1‑149 (MQ=255)
                                                                                                                              aaCATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGaa                     <  3:665042/149‑1 (MQ=255)
                                                                                                                              aaCATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGaa                     >  1:169986/1‑149 (MQ=255)
                                                                                                                                 aTTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGcc                  >  3:270314/1‑149 (MQ=255)
                                                                                                                                  ttACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCg                 <  2:500416/149‑1 (MQ=255)
                                                                                                                                         tCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGCTGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGc          >  2:441825/1‑149 (MQ=255)
                                                                                                                                          cacaTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGcc         <  1:281283/149‑1 (MQ=255)
                                                                                                                                           agaTCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCg        >  3:574986/3‑149 (MQ=255)
                                                                                                                                           acaTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCg        <  4:270314/149‑1 (MQ=255)
                                                                                                                                            gaTCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCgg       >  2:351720/2‑149 (MQ=255)
                                                                                                                                            caTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGTCTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCgg       <  1:478552/149‑1 (MQ=255)
                                                                                                                                            caTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCgg       <  2:169986/149‑1 (MQ=255)
                                                                                                                                            caTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCgg       <  4:48620/149‑1 (MQ=255)
                                                                                                                                            caTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCACGAAGCCGAAAAAGCCgg       <  3:679352/149‑1 (MQ=255)
                                                                                                                                             aTCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCgg       <  1:539696/148‑1 (MQ=255)
                                                                                                                                             aTCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCgg       >  2:539696/1‑148 (MQ=255)
                                                                                                                                              tCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCAGGCGATGGTTAAATCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCTTAGAAGCGGACGACGTTATCGGTACTCTTGCGCGCGAAGCCGAAAAAGCCGGGc     >  4:658818/1‑149 (MQ=255)
                                                                                                                                              tCGTCCGGCATTGGCGGGCGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGc     >  2:225857/1‑149 (MQ=255)
                                                                                                                                              tCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGc     >  2:362696/1‑149 (MQ=255)
                                                                                                                                                tccccccAATGCGGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGAGGGTGAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGGAGCGGACGACGGTAACGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGt   >  2:308132/2‑149 (MQ=255)
                                                                                                                                                 cccgccAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTc  >  3:658801/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCGAGCCGACCGGTGCGATGTATGGTGTCCTCAACATGCTGCGCAGTCTGATCATGCAATATAAACCGACGCATGCAGCGGTGGTCTTTGACGCCAAGGGAAAAACCTTTCGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTC  >  NC_000913/4047063‑4047356

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: