Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 100,117 T→G 11.4% N98H (AAC→CAC) ‡ PP_0094 ← 5'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947100,1170TG11.4% 163.2 / 11.4 60N98H (AAC→CAC) ‡PP_00945'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base G (4/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ATCAGCTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTGCAGCAGGG  >  NC_002947/99969‑100263
                                                                                                                                                    |                                                                                                                                                  
aTCAGCTCGACCTTGTTGGTGCGCGTGCCGTCCTGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGtt                                                                                                                                                    <  8:267145/149‑1 (MQ=255)
   aGCTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAg                                                                                                                                                 >  4:111640/1‑149 (MQ=255)
    gCTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAgg                                                                                                                                                <  3:251633/149‑1 (MQ=255)
        gACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGcc                                                                                                                                             <  4:43398/148‑1 (MQ=255)
        gACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGcc                                                                                                                                             <  6:224461/148‑1 (MQ=255)
           cTTGTTGGTGCGCGTGGCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTc                                                                                                                                           >  6:327392/1‑147 (MQ=255)
               ttgGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGCAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGAACGGGCAGGGCCGCACGct                                                                                                                                      <  5:320884/148‑2 (MQ=255)
                              tCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTc                                                                                                                      <  5:149926/149‑1 (MQ=255)
                                  gTTCACTGCCGTAGATCACCTTGAAGTGATTATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAac                                                                                                                  <  3:216509/149‑1 (MQ=255)
                                  gTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAac                                                                                                                  <  3:92063/149‑1 (MQ=255)
                                       cTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAg                                                                                                             <  5:184583/149‑1 (MQ=255)
                                         gCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGt                                                                                                           >  7:197433/1‑149 (MQ=255)
                                           cGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGttt                                                                                                         >  8:181592/1‑149 (MQ=255)
                                                 tCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTa                                                                                                    >  8:206503/1‑148 (MQ=255)
                                                     cTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGAc                                                                                               >  6:291971/1‑149 (MQ=255)
                                                                  tCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCgg                                                                                    >  6:201380/1‑147 (MQ=255)
                                                                    gAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACt                                                                                                                                   <  3:148/98‑1 (MQ=255)
                                                                    gAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACt                                                                                                                                   >  4:148/1‑98 (MQ=255)
                                                                          aaaGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGGAAAGACCGGTCCCGCGGAAGAGCTcgc                                                                          >  5:195879/1‑149 (MQ=255)
                                                                                ccGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGa                                                                    <  7:206503/149‑1 (MQ=255)
                                                                                  gAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAg                                                                  >  3:22822/1‑149 (MQ=255)
                                                                                    gCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGt                                                                 >  4:3074/1‑148 (MQ=255)
                                                                                        tttCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGaa                                                                                  >  1:15167/1‑127 (MQ=255)
                                                                                        tttCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGaa                                                                                  <  2:15167/127‑1 (MQ=255)
                                                                                           cgcgGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCa                                                         >  2:275596/1‑149 (MQ=255)
                                                                                                gCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCa                                                    <  3:3074/149‑1 (MQ=255)
                                                                                                gCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCa                                                    <  5:327392/149‑1 (MQ=255)
                                                                                                  gAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAgg                                                  >  3:69687/1‑149 (MQ=255)
                                                                                                      tCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCtt                                              >  3:151030/1‑149 (MQ=255)
                                                                                                       cccATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTg                                             >  2:66484/1‑149 (MQ=255)
                                                                                                        ccATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCtt                                              <  8:197433/147‑1 (MQ=255)
                                                                                                        ccATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTgg                                            >  8:107275/1‑149 (MQ=255)
                                                                                                        ccATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGATTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTgg                                            >  8:288432/1‑149 (MQ=255)
                                                                                                          aTGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCt                                          >  5:293813/1‑149 (MQ=255)
                                                                                                          aTGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCt                                          >  3:117058/1‑149 (MQ=255)
                                                                                                              tttGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTtcgt                                      <  4:21615/149‑1 (MQ=255)
                                                                                                               ttGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTtcgtc                                     >  3:228603/1‑149 (MQ=255)
                                                                                                                  gAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGa                                   <  6:293813/148‑1 (MQ=255)
                                                                                                                     gTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCt                               >  6:20597/1‑149 (MQ=255)
                                                                                                                     gTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGACGAAGc                                <  4:69687/148‑1 (MQ=255)
                                                                                                                      tGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCt                                          >  5:333537/1‑137 (MQ=255)
                                                                                                                      tGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCt                                          <  6:333537/137‑1 (MQ=255)
                                                                                                                      tGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTg                              >  7:135950/1‑149 (MQ=255)
                                                                                                                      tGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTg                              <  3:111640/149‑1 (MQ=255)
                                                                                                                             gTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACt                       >  4:67035/1‑149 (MQ=255)
                                                                                                                                cAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCa                    <  5:291971/149‑1 (MQ=255)
                                                                                                                                 aGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCAt                   <  5:201380/149‑1 (MQ=255)
                                                                                                                                    ttgaggCCCTGAACGGGCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGGTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCCTCGTCGAAGCTGTAGAACTGCATGa                 >  7:175576/5‑148 (MQ=255)
                                                                                                                                    ttgaggCCCTGAACGGGCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCAcc                                                            <  6:286809/101‑1 (MQ=255)
                                                                                                                                    ttgaggCCCTGAACGGGCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCAcc                                                            >  5:286809/5‑105 (MQ=255)
                                                                                                                                    ttgaggCCCTGAACGGGCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGa                 >  8:148245/5‑148 (MQ=255)
                                                                                                                                      gaggCCCTGAACGGGCAGGGCCGCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTgg                                            <  4:83654/117‑1 (MQ=255)
                                                                                                                                      gaggCCCTGAACGGGCAGGGCCGCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTgg                                            >  3:83654/3‑119 (MQ=255)
                                                                                                                                       cGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGc             <  7:181592/149‑1 (MQ=255)
                                                                                                                                         gCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGt                             >  7:111779/1‑131 (MQ=255)
                                                                                                                                         gCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGt                             <  8:111779/131‑1 (MQ=255)
                                                                                                                                         gCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCt           >  2:779/1‑149 (MQ=255)
                                                                                                                                           ccTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTgc         <  4:117058/149‑1 (MQ=255)
                                                                                                                                              gCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTgcag       >  1:334372/1‑148 (MQ=255)
                                                                                                                                                 cGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTGCAGCAgg   <  2:334372/149‑1 (MQ=255)
                                                                                                                                                  gTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTGCAGCAggg  <  4:22822/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATCAGCTCGACCTTGTTGGTGCGCGTGCCGTCCAGTTCACTGCCGTAGATCACCTTGAAGTGATGATCGAACGCAAAGTGCCGAGCGATTTCGCGGGCGAACTCCCATGGTTTGGAGGTGGCAATGTACAGCGTGCGGCCCTGCCCGTTCAGGGCCTCAAGCAACTCCGGCACACCTTCGAACACGAGGTTTTCGTAAAGACCGGTCACGCGGAAGCGCTCGCGGTAGAAGTTCACCGCATCCCAGGCCTTGGCTTCGTCGAAGCTGTAGAACTGCATGAAAGCCTGCAGCAGGG  >  NC_002947/99969‑100263

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