Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,007,239 T→C 6.5% Q41Q (CAA→CAG PP_1790 ← acylneuraminate cytidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,007,2390TC 6.5% 177.9 / 2.7 61Q41Q (CAA→CAGPP_1790acylneuraminate cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base C (2/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

ACTGGTCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCTTT  >  NC_002947/2007091‑2007383
                                                                                                                                                      |                                                                                                                                                
aCTGGTCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGt                                                                                                                                                  <  7:7970/149‑1 (MQ=255)
 cTGGTCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGtt                                                                                                                                                 >  3:45285/1‑149 (MQ=255)
    gTCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCTTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGc                                                                                                                                               <  1:212012/148‑1 (MQ=255)
    gTCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTg                                                                                                                                                <  7:41845/147‑1 (MQ=255)
     tCAGTGCTTGTCTGGCTTGCAGATAGTTCTTCTCGTCTCTTATGTAGTTTGATTTTATTTTTTTGAAGCGATGAGGCGTATTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCtt                                                                                                                                             <  4:199972/149‑1 (MQ=255)
     tCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCt                                                                                                                                              >  7:61889/1‑148 (MQ=255)
           cTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAACCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAACTGCTTCaagc                                                                                                                                        <  1:238224/148‑3 (MQ=255)
              gTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGa                                                                                                                                     <  2:244043/148‑1 (MQ=255)
              gTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGAt                                                                                                                                    >  8:115068/1‑149 (MQ=255)
              gTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAACTGCTTCaagcaga                                                                                                                                     <  4:251702/148‑6 (MQ=255)
                   gctggAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGAttt                                                                                                                                  <  2:249526/146‑1 (MQ=255)
                   gctggAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTtcaa                                                                                                                               <  5:334131/149‑1 (MQ=255)
                       gAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTtcaatcaa                                                                                                                           >  3:250902/1‑149 (MQ=255)
                       gAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTtcaatc                                                                                                                             <  2:242083/147‑1 (MQ=255)
                          gATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGcc                                                                                                                        <  8:61889/149‑1 (MQ=255)
                          gATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGcc                                                                                                                        >  2:59340/1‑149 (MQ=255)
                           ctAGGGCTTCTCGTCGCTGATGAAGTTTGATTTTATTATTTTGAAGCGCGGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTAGGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAAGTTGATTTCAATCAAGcc                                                                                                                        <  8:224454/147‑1 (MQ=255)
                                    ctcGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGGTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTg                                                                                                               <  5:10899/148‑1 (MQ=255)
                                       gtcgCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTgg                                                                                                              <  2:228017/146‑1 (MQ=255)
                                         cgCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGAGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGgccgc                                                                                                          >  6:179116/1‑148 (MQ=255)
                                          gCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGgccgcc                                                                                                         >  2:295791/1‑148 (MQ=255)
                                                   gTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGtttt                                                                                               <  5:81548/149‑1 (MQ=255)
                                                   gTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGtttt                                                                                               <  6:9988/149‑1 (MQ=255)
                                                    tttGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTGCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGttt                                                                                                >  1:249777/1‑148 (MQ=255)
                                                    tttGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAAt                                                                                                     <  5:232310/142‑1 (MQ=255)
                                                    tttGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAAt                                                                                                     >  6:232310/1‑142 (MQ=255)
                                                      tGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTtat                                                                                             <  4:45285/148‑1 (MQ=255)
                                                      tGATTTTAATATTTTGAAGCGATGAGGCCTAGTCAAGGATGCTTTCGTCGTTGGGGCTCTGTACGAGGGAATCCACAGCGGCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTtat                                                                                             <  7:229034/148‑1 (MQ=255)
                                                       gATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATTGTTTtata                                                                                            <  5:179116/148‑1 (MQ=255)
                                                          tttATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATAAAGCCGAATTCGTGGCCGCCAAATGGTTTTATAtttt                                                                                        >  7:207732/1‑149 (MQ=255)
                                                             tttattTTCAGCGAAGAGGCGTATTCAATGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCACCAAATGGTTTTATATTTTTAtt                                                                                    <  4:232263/148‑1 (MQ=255)
                                                             attattTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTAt                                                                                     >  7:136271/1‑149 (MQ=255)
                                                                attTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTTGTGCTCAGTACGATGGAATCAACAGCTTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTc                                                                                   <  7:115068/148‑1 (MQ=255)
                                                                 ttttGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCg                                                                                  >  6:208496/1‑148 (MQ=255)
                                                                       aGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGc                                                                            <  4:250902/148‑1 (MQ=255)
                                                                                gCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTtata                                                                                            <  7:70341/123‑1 (MQ=255)
                                                                                gCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTtata                                                                                            >  8:70341/1‑123 (MQ=255)
                                                                                gCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGATGcc                                                                  >  4:79560/1‑149 (MQ=255)
                                                                                  gTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTAGCGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGgccgcc                                                                                                         >  1:225806/1‑109 (MQ=255)
                                                                                  gTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGgccgcc                                                                                                         <  2:225806/108‑1 (MQ=255)
                                                                                        aaGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATTGTTTTATATTTTTATTCGGAATGCATTCGCTGCCTTTTCgg                                                           >  4:307024/1‑148 (MQ=255)
                                                                                        aaGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggc                                                          >  6:256204/1‑149 (MQ=255)
                                                                                        aaGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCgg                                                           >  4:69592/1‑148 (MQ=255)
                                                                                          ggATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCtttt                                                              <  6:163874/143‑1 (MQ=255)
                                                                                          ggATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCtttt                                                              >  5:163874/1‑143 (MQ=255)
                                                                                           gATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTc                                                             <  7:283660/143‑1 (MQ=255)
                                                                                           gATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTc                                                             >  8:283660/1‑143 (MQ=255)
                                                                                                tttCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGCATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggaagagatcg                                                  >  7:218658/1‑143 (MQ=255)
                                                                                                tttCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggca                                                         <  8:218658/142‑1 (MQ=255)
                                                                                                     gtcgtTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGaaaa                                               <  1:59340/148‑1 (MQ=255)
                                                                                                     gtcgtTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAac                                              <  1:295791/149‑1 (MQ=255)
                                                                                                     gtcgtTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAAATTCATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGaaaa                                               >  5:241607/1‑148 (MQ=255)
                                                                                                           ggTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACAcgtc                                         <  1:39553/148‑1 (MQ=255)
                                                                                                                  aGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGtttttt                                  <  5:208496/148‑1 (MQ=255)
                                                                                                                  aGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGtttttt                                  <  8:136271/148‑1 (MQ=255)
                                                                                                                   gTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCa                                <  5:256204/149‑1 (MQ=255)
                                                                                                                           ggAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCata                              <  7:20058/144‑1 (MQ=255)
                                                                                                                           ggAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCata                              >  8:20058/1‑144 (MQ=255)
                                                                                                                                         cTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCa            >  3:2613/1‑148 (MQ=255)
                                                                                                                                            gttGAAGCAACTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATaaa        >  5:229450/2‑149 (MQ=255)
                                                                                                                                             ttGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAAc       >  3:183644/1‑149 (MQ=255)
                                                                                                                                             ttGAAGCAACTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGACGTTTTTTCATATTTCATTTCAGTGGCCCATAAAc       >  3:244849/1‑149 (MQ=255)
                                                                                                                                                  gCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCttt  >  5:257779/1‑149 (MQ=255)
                                                                                                                                                      ttGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggcagggca                                                    <  3:222575/95‑1 (MQ=255)
                                                                                                                                                      ttGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggcagggca                                                    >  4:222575/1‑95 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
ACTGGTCAGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCTTT  >  NC_002947/2007091‑2007383

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: