Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,026,752 A→T 12.5% D21V (GAC→GTC)  ocd → ornithine cyclodeaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,026,7520AT12.5% 146.3 / 11.7 56D21V (GAC→GTC) ocdornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base A (28/21);  new base T (2/5);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

TGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTAT  >  NC_002947/5026610‑5026899
                                                                                                                                              |                                                                                                                                                   
tGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTAt                                                                                                                                               <  8:231848/149‑1 (MQ=255)
    gACATGGCTCATCACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGAt                                                                                                                                           <  4:172156/149‑1 (MQ=255)
    gACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                            <  5:100979/148‑1 (MQ=255)
    gACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                            <  2:161553/148‑1 (MQ=255)
     aCATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGagc                                                                                                                                          <  8:46853/149‑3 (MQ=255)
        tGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAg                                                                                                                                       >  6:62027/1‑149 (MQ=255)
            tCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                            <  3:166111/140‑1 (MQ=255)
            tCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                            >  4:166111/1‑140 (MQ=255)
                 acCAACGCATCATCCTGGAAGCTGTGCGCGCCCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGCTGTAGTGGAAAACg                                                                                                                              >  5:131660/1‑149 (MQ=255)
                 acCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCACCGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACg                                                                                                                              >  7:84196/1‑149 (MQ=255)
                   cAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATacagac                                                                                                                                                    >  5:55876/1‑125 (MQ=255)
                   cAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATacagac                                                                                                                                                    <  6:55876/125‑1 (MQ=255)
                     aCGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATacagaca                                                                                                                                                   >  5:152829/1‑124 (MQ=255)
                     aCGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATacagaca                                                                                                                                                   <  6:152829/124‑1 (MQ=255)
                          tcatcCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCAt                                                                                                                     <  1:327817/149‑1 (MQ=255)
                            atcCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCAttt                                                                                                                         <  3:138509/143‑1 (MQ=255)
                            atcCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCAttt                                                                                                                         >  4:138509/1‑143 (MQ=255)
                                tGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTgg                                                                                                               >  8:326784/1‑149 (MQ=255)
                                   aaGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCa                                                                                                            <  7:326784/149‑1 (MQ=255)
                                         agcgcGACCGCCGCCCAGCACGTGCTTAAAGCGGTGCCATGTGAGCGAAGTCCTTCTGGTTAAGCAACGCGAAATCCGTGAGTGCGTAACGCGCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCa                                                                                                      <  6:131660/148‑1 (MQ=255)
                                               aCCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTgg                                                                                                <  6:221081/149‑1 (MQ=255)
                                               aCCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTgg                                                                                                >  5:221081/1‑149 (MQ=255)
                                                   ccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCAt                                                                                            >  3:70262/1‑149 (MQ=255)
                                                   ccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCAt                                                                                            >  3:157381/1‑149 (MQ=255)
                                                       ccAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGATCATGCCg                                                                                        >  8:41197/1‑149 (MQ=255)
                                                        cAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCACCGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGt                                                                                       <  8:84196/149‑1 (MQ=255)
                                                                     aaGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGt                                                                           >  5:368044/1‑148 (MQ=255)
                                                                         ggTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGATCATGCCGTCAATCCTCCGGTTAGAc                                                                      <  7:41197/149‑1 (MQ=255)
                                                                             ccATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACAtct                                                                  >  7:61630/1‑149 (MQ=255)
                                                                                tgtgAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTcc                                                                >  4:323919/1‑148 (MQ=255)
                                                                                tgtgAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTcc                                                                >  8:120267/1‑148 (MQ=255)
                                                                                  tgAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAg                                                             <  8:61630/149‑1 (MQ=255)
                                                                                         aGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATg                                                      <  4:157381/149‑1 (MQ=255)
                                                                                           tCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGc                                                    >  1:364865/1‑149 (MQ=255)
                                                                                              ttCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACAtct                                                                  >  1:195756/1‑132 (MQ=255)
                                                                                              ttCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACAtct                                                                  <  2:195756/132‑1 (MQ=255)
                                                                                                  gCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTg                                              >  5:332034/1‑148 (MQ=255)
                                                                                                    ttAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGAc                                           >  7:120737/1‑149 (MQ=255)
                                                                                                      aaCCAACGCGAAATCCGTGAGTGCTTAAAGCTCGATACAGACTGTATCGATGTAGTGGAAAACGCATTTGCATCGCTTGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGAc                                           <  7:253134/147‑1 (MQ=255)
                                                                                                        ccAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTcc                                                                <  4:135734/124‑1 (MQ=255)
                                                                                                        ccAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTcc                                                                >  3:135734/1‑124 (MQ=255)
                                                                                                        ccAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTc                                        >  8:81162/1‑148 (MQ=255)
                                                                                                             gcgAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAAcc                                  >  7:67564/1‑149 (MQ=255)
                                                                                                               gAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGc                                <  5:62027/149‑1 (MQ=255)
                                                                                                               gAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGACGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGc                                >  4:305830/1‑149 (MQ=255)
                                                                                                                 aaTCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCt                               >  3:19299/1‑148 (MQ=255)
                                                                                                                     cGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCt                          >  8:208887/1‑149 (MQ=255)
                                                                                                                       tGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGc                        <  7:208887/149‑1 (MQ=255)
                                                                                                                         aGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCa                      <  8:120737/149‑1 (MQ=255)
                                                                                                                              cactacatCGATACTGTCTGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCt                 >  3:29855/8‑149 (MQ=255)
                                                                                                                                ctacatCGATACTGTCTGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTg               <  4:29855/144‑1 (MQ=255)
                                                                                                                                  cGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCtt                >  2:298278/1‑146 (MQ=255)
                                                                                                                                  cGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTgag             <  4:19299/149‑1 (MQ=255)
                                                                                                                                   gCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTgag             >  7:95747/1‑148 (MQ=255)
                                                                                                                                     tCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAAt          >  3:52316/1‑149 (MQ=255)
                                                                                                                                       gATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACTATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAAttt        >  8:105516/1‑149 (MQ=255)
                                                                                                                                       gATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCt                               <  8:25830/126‑1 (MQ=255)
                                                                                                                                       gATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCt                               >  7:25830/1‑126 (MQ=255)
                                                                                                                                        aTACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTc       <  8:95747/149‑1 (MQ=255)
                                                                                                                                         tACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCg      <  3:323919/149‑1 (MQ=255)
                                                                                                                                             gacagTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTAt  >  4:37304/1‑149 (MQ=255)
                                                                                                                                              acagTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTAt  >  3:274457/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTAT  >  NC_002947/5026610‑5026899

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: ‑