Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,824,518 A→C 18.7% E44A (GAG→GCG) ‡ PP_2478 → isoquinoline 1‑oxidoreductase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,824,5180AC18.7% 127.6 / 25.3 64E44A (GAG→GCG) ‡PP_2478isoquinoline 1‑oxidoreductase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base C (5/7);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

AGGTTGAGGAGGCGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTATTGCCATGA  >  NC_002947/2824370‑2824658
                                                                                                                                                    |                                                                                                                                            
aGGTTGAGGAGGCGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCACGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGa                                                                                                                                              <  1:119048/149‑1 (MQ=255)
  gTTGAGGAGGCGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAg                                                                                                                                             <  6:207923/148‑1 (MQ=255)
  gTTGAGGAGGCGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAg                                                                                                                                             >  5:207923/1‑148 (MQ=255)
     gaggagGCGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAggcgg                                                                                                                                         >  2:133023/1‑149 (MQ=255)
     gaggagGCGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAggcgg                                                                                                                                         >  2:138130/1‑149 (MQ=255)
      aggaggCGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAggcgg                                                                                                                                         <  1:133023/148‑1 (MQ=255)
      aggaggCGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAggcgg                                                                                                                                         <  1:138130/148‑1 (MQ=255)
            cGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCTCGGCGGNTGCg                                                                                                                                    <  5:123349/147‑1 (MQ=255)
            cGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCTCGGCGGCTGCg                                                                                                                                    >  6:123349/1‑147 (MQ=255)
                agaTGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCTCGGCGGCTGCgcttttg                                                                                                                              <  5:157707/149‑7 (MQ=255)
                    gAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCTCGGCGGCTGCgc                                                                                                                                   >  3:332704/1‑139 (MQ=255)
                    gAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCTCGGCGGCTGCgc                                                                                                                                   <  4:332704/140‑2 (MQ=255)
                    gAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGc                                                                                                                          >  8:102248/1‑149 (MQ=255)
                          gCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGACGCCGCCGCGGCAGGGCGGCTTggg                                                                                                                     >  2:372101/1‑147 (MQ=255)
                            cTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGAt                                                                                                                  <  5:311855/149‑1 (MQ=255)
                            cTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGAt                                                                                                                  <  1:270436/149‑1 (MQ=255)
                                 gAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTc                                                                                                             >  4:301492/1‑149 (MQ=255)
                                        cTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAg                                                                                                      <  5:316527/149‑1 (MQ=255)
                                            tGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACTTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCTCGGCGGCTGCgctttt                                                                                                                               <  1:285548/120‑6 (MQ=255)
                                            tGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCTCGGCGGCTGCgctttt                                                                                                                               >  2:285548/1‑115 (MQ=255)
                                            tGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGt                                                                                                   >  2:170721/1‑148 (MQ=255)
                                                atGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCTCGGCGGCTGc                                                                                                                                     <  6:292661/110‑1 (MQ=255)
                                                atGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACgg                                                                                              >  1:14858/1‑149 (MQ=255)
                                                atGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACgg                                                                                              >  2:67797/1‑149 (MQ=255)
                                                atGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACg                                                                                               <  7:159557/148‑1 (MQ=255)
                                                     ccGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCt                                                                                         >  1:15249/1‑149 (MQ=255)
                                                       gTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCttt                                                                                       <  7:150132/149‑1 (MQ=255)
                                                         cgccgcCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGt                                                                                     <  3:167006/149‑1 (MQ=255)
                                                          gccgccTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTc                                                                                    <  1:170721/149‑1 (MQ=255)
                                                                      aGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGcccc                                                                        >  5:118131/1‑149 (MQ=255)
                                                                               gTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTg                                                               <  5:85224/149‑1 (MQ=255)
                                                                                tATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTgg                                                              >  1:176044/1‑149 (MQ=255)
                                                                                          aGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGa                                                     >  4:288266/1‑148 (MQ=255)
                                                                                               cTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGt                                                             <  1:122877/135‑1 (MQ=255)
                                                                                               cTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGt                                                             >  2:122877/1‑135 (MQ=255)
                                                                                                     gtTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGaaaa                                         >  5:247947/1‑149 (MQ=255)
                                                                                                        gAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAAtt                                       >  1:342311/1‑148 (MQ=255)
                                                                                                        gAGTACCTGGCTACCGCCGCTCGTGGCCAAAAGCGCCCCCGCCCAGGGCGGAGCCCGACCTCCGGGTGGTGCTGACTCCGCCGAGGGTTACGGTGCCTATCTACCCCTAGGCCCCGTTGCCGGGGTCACGGTGATCTCggcccaaa                                         >  8:121825/1‑139 (MQ=255)
                                                                                                         aGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGa                                     <  5:47235/149‑1 (MQ=255)
                                                                                                           tACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGaaa                                   <  2:342311/149‑1 (MQ=255)
                                                                                                             ccTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAAt                                  <  2:176044/148‑1 (MQ=255)
                                                                                                              cTGGCTACCGCCTCTGGTTTCCAAAAGCTCAGCCGCCGAGGCGGCTGCGTCAGGGAGTATTTGTGATAACTACGCGGAGTGTTACGGTTACTTTGTCCGCATTGTCCCCGATGGCGTGGTCACGGTGATCTCGACGAAAATTGAAATgt                                <  2:121957/149‑2 (MQ=255)
                                                                                                              cTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATgg                                <  3:288266/149‑1 (MQ=255)
                                                                                                                  cTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCa                            <  2:14858/149‑1 (MQ=255)
                                                                                                                  cTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCa                            <  3:301492/149‑1 (MQ=255)
                                                                                                                    aCCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAGATGGGCCAgg                          >  7:271799/1‑149 (MQ=255)
                                                                                                                      cgccgcTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGt                        >  7:32347/1‑149 (MQ=255)
                                                                                                                        ccgcTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTAc                                                                                                >  6:366695/1‑75 (MQ=255)
                                                                                                                        ccgcTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTAc                                                                                                <  5:366695/75‑1 (MQ=255)
                                                                                                                        ccgcTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTgc                      >  3:21548/1‑149 (MQ=255)
                                                                                                                         cgcTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTgcg                     <  6:247947/149‑1 (MQ=255)
                                                                                                                                ggCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGCCGAGGTCACGGTGATCACGCCGAAAATTGAAATGGGCAAAGGTGCGCAAACCg              >  6:281000/1‑149 (MQ=255)
                                                                                                                                  ccAAAAGGGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGt            <  1:372101/149‑1 (MQ=255)
                                                                                                                                  ccAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCgg             >  3:314547/1‑148 (MQ=255)
                                                                                                                                   cAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTa           >  5:170735/1‑149 (MQ=255)
                                                                                                                                    aaaaGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGCTCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGAGCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAACCCgggag          >  4:286168/1‑146 (MQ=255)
                                                                                                                                     aaaGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTAtt         <  1:67797/149‑1 (MQ=255)
                                                                                                                                      tgccgcAGCCGCCTCGGCGGCTGCAGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTAtt         >  7:153314/4‑148 (MQ=255)
                                                                                                                                      tgccgcAGCCGCCTCGGCGGCTGCAGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTAtt         <  8:153314/145‑1 (MQ=255)
                                                                                                                                         cgcAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTATTGCCa     >  3:220000/1‑149 (MQ=255)
                                                                                                                                          gcAGCCGCCTCGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAAt                                  <  2:195180/119‑1 (MQ=255)
                                                                                                                                          gcAGCCGCCTCGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAAt                                  >  1:195180/1‑119 (MQ=255)
                                                                                                                                          gcAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTATTGCCat    <  7:102248/149‑1 (MQ=255)
                                                                                                                                            aGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCTGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTATTGCCatga  >  1:288045/1‑149 (MQ=255)
                                                                                                                                            aGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTATTGCCatga  >  5:247976/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AGGTTGAGGAGGCGTGAGATGAAAAAGCCTAACGAAGTGACTGTTGATATGTCCCGTCGCCGCCTCTTGCAGGGCAGCGGTATCGCCCTCAGTGGCTTGGTGTTGAGTACCTGGCTACCGCCGCTGGTGGCCAAAAGCGCAGCCGCCGAGGCGGCTGCGGCAGGTCGGCTTGGTGATCACTCCGCGGAGGGTTACGGTGCCTTTGTCCGCATTGGCCCCGATGGCGTGGTCACGGTGATCTCGCCGAAAATTGAAATGGGCCAGGGTGCGCAAACCGGTATTGCCATGA  >  NC_002947/2824370‑2824658

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