Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,454,292 C→A G122G (GGC→GGA paaA → ring 1,2‑phenylacetyl‑CoA epoxidase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,454,2920CA98.6% 225.1 / NA 73G122G (GGC→GGApaaAring 1,2‑phenylacetyl‑CoA epoxidase subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (39/32);  minor base G (0/1);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGCTGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTTTGAAGCCTGC  >  NC_000913/1454162‑1454424
                                                                                                                                  |                                                                                                                                    
tctATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCtggt                                                                                                                              <  1:440341/139‑1 (MQ=255)
  tataGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCtggtgg                                                                                                                            >  2:795590/1‑139 (MQ=255)
      gcgcCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATgg                                                                                                                        >  1:915331/1‑139 (MQ=255)
      gcgcCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATgg                                                                                                                        >  2:1249095/1‑139 (MQ=255)
         ccgccgAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGgcgc                                                                                                                     >  2:1140443/1‑139 (MQ=255)
             ggAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGGTGGCTGGTGGATCgctgcgcc                                                                                                                 <  2:1195025/138‑7 (MQ=255)
              gAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGa                                                                                                                <  2:555980/139‑1 (MQ=255)
                      gggCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAAt                                                                                                        >  2:576125/1‑139 (MQ=255)
                                 gagaAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGtt                                                                                             >  2:299963/1‑139 (MQ=255)
                                     agaCATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGc                                                                                         <  1:576125/139‑1 (MQ=255)
                                     agaCATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGc                                                                                         <  2:455812/139‑1 (MQ=255)
                                         aTCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCa                                                                                     >  2:959684/1‑139 (MQ=255)
                                          tCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                                    <  1:959684/139‑1 (MQ=255)
                                          tCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                                    <  1:299963/139‑1 (MQ=255)
                                          tCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTCGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                                    <  1:1140443/139‑1 (MQ=255)
                                           cTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGTTGGCGTTATGCCGCAct                                                                                   >  1:527316/1‑139 (MQ=255)
                                            tACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCActt                                                                                  >  1:1288991/1‑139 (MQ=255)
                                            tACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCActt                                                                                  >  2:1022985/1‑139 (MQ=255)
                                              ccAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCActtat                                                                                >  1:839246/1‑137 (MQ=255)
                                                  aaaaTGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATg                                                                            <  2:527316/139‑1 (MQ=255)
                                                  aaaaTGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATg                                                                            <  1:108518/139‑1 (MQ=255)
                                                       gCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCg                                                                       <  1:795590/139‑1 (MQ=255)
                                                        cTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGt                                                                      >  1:10688/1‑139 (MQ=255)
                                                             cggacggaTGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGt                                                                                                   <  1:553470/105‑1 (MQ=255)
                                                             cggacggaTGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGt                                                                                                   >  2:553470/1‑105 (MQ=255)
                                                                     tGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                                    <  2:828914/112‑1 (MQ=255)
                                                                     tGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                                    >  1:828914/1‑112 (MQ=255)
                                                                              cctccATCTTTGATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGTTGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCg                                                                       >  1:934726/1‑116 (MQ=255)
                                                                              cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCg                                                                       <  2:934726/116‑1 (MQ=255)
                                                                              cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGCTGGTGGATGGCGCAGCGa                                                                                                                >  2:861413/1‑74 (MQ=255)
                                                                              cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGCTGGTGGATGGCGCAGCGa                                                                                                                <  1:861413/74‑1 (MQ=255)
                                                                                    tCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCaa                                          >  2:1386021/1‑139 (MQ=255)
                                                                                      tttAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAaag                                        >  1:222532/1‑139 (MQ=255)
                                                                                         aaTTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAaagaag                                     >  1:684981/1‑139 (MQ=255)
                                                                                          aTTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTATTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAagaaga                                    <  2:839246/139‑1 (MQ=255)
                                                                                           ttATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAagag                                   >  2:1414055/1‑139 (MQ=255)
                                                                                              tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGGTGGCGCTGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATg                                                                            <  2:1399192/95‑1 (MQ=255)
                                                                                              tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATg                                                                            >  1:1399192/1‑95 (MQ=255)
                                                                                              tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATg                                                                            <  1:1128335/95‑1 (MQ=255)
                                                                                              tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATg                                                                            >  2:1128335/1‑95 (MQ=255)
                                                                                              tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTg                                                     >  1:1457290/1‑118 (MQ=255)
                                                                                              tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTg                                                     <  2:1457290/118‑1 (MQ=255)
                                                                                              tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCa                                           >  1:706482/1‑128 (MQ=255)
                                                                                              tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCa                                           <  2:706482/128‑1 (MQ=255)
                                                                                              tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGt                                >  1:1468980/1‑139 (MQ=255)
                                                                                               ccGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGtt                               >  2:339209/1‑139 (MQ=255)
                                                                                                cGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGa                                                                                                                           <  2:1066635/46‑1 (MQ=255)
                                                                                                cGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATg                                                                            <  2:673328/93‑1 (MQ=255)
                                                                                                cGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATg                                                                            >  1:673328/1‑93 (MQ=255)
                                                                                                cGACATTGAGGTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGa                                                                                                                           >  1:1066635/1‑46 (MQ=38)
                                                                                                   cATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCa                           >  2:759129/1‑139 (MQ=255)
                                                                                                     ttGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCAcc                         <  1:339209/139‑1 (MQ=255)
                                                                                                     ttGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCAcc                         <  1:1414055/139‑1 (MQ=255)
                                                                                                        aGTTGGGCTGATATCGGTGTGATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCAcc                         >  1:947288/1‑136 (MQ=255)
                                                                                                        aGTTGGGCTGATATCGGTATTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCAcc                         <  2:947288/136‑1 (MQ=255)
                                                                                                          ttGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGt                    >  1:1143576/1‑139 (MQ=255)
                                                                                                            gggCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCa                  <  1:1386021/139‑1 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCaa                                                                                                         <  2:108889/46‑1 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCaa                                                                                                         >  1:108889/1‑46 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTa                                                                                            <  2:843498/59‑1 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTa                                                                                            >  1:843498/1‑59 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                                    <  1:980428/67‑1 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCg                                                                       <  2:1267123/80‑1 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCg                                                                       >  1:1267123/1‑80 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGtt            >  1:977982/1‑139 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGtt            >  2:1333037/1‑139 (MQ=255)
                                                                                                                  atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATAGTCAATCAGGTGGCGTTATGCCGCAc                                                                                    >  2:980428/1‑67 (MQ=255)
                                                                                                                    atCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCg                                                                       >  2:1009481/1‑78 (MQ=255)
                                                                                                                    atCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCg                                                                       <  1:1009481/78‑1 (MQ=255)
                                                                                                                       ggTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTTTGaa        >  2:1113313/1‑138 (MQ=255)
                                                                                                                       ggTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTTTGaa        <  1:1113313/138‑1 (MQ=255)
                                                                                                                        gtgtTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTTTGAAGc      <  1:759129/139‑1 (MQ=255)
                                                                                                                            tATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTTTGAAGCCTGc  >  2:28261/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                    
TCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGCTGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTTTGAAGCCTGC  >  NC_000913/1454162‑1454424

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