Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,007,239 T→C 7.3% 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 7.3% 198.1 / 4.8 68Q41Q (CAA→CAGPP_1790acylneuraminate cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (3/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

AGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCTTTGGCA  >  NC_002947/2007098‑2007387
                                                                                                                                             |                                                                                                                                                    
aGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTc                                                                                                                                                >  3:250652/1‑148 (MQ=255)
aGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTc                                                                                                                                                <  7:77282/148‑1 (MQ=255)
   gCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAACTGCTTCaag                                                                                                                                             <  2:298366/148‑2 (MQ=255)
    cTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACtt                                                                                                                                           >  6:146827/1‑149 (MQ=255)
     tGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGAACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTg                                                                                                                                          >  4:268923/1‑149 (MQ=255)
      ggTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGa                                                                                                                                         <  8:164271/149‑1 (MQ=255)
        tCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGAtt                                                                                                                                       <  8:189692/149‑1 (MQ=255)
        gctggctggAAGATAGGTCTTGACGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGGGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGAtt                                                                                                                                       <  4:250652/148‑1 (MQ=255)
            gctggAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGACGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTtcaa                                                                                                                                   <  5:30359/149‑1 (MQ=255)
             ctggAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTtcaat                                                                                                                                  <  3:11045/149‑1 (MQ=255)
                   gATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGcc                                                                                                                            <  4:172002/149‑1 (MQ=255)
                         tcttctCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTc                                                                                                                      >  3:127690/1‑149 (MQ=255)
                         tcttctCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTc                                                                                                                      >  8:221629/1‑149 (MQ=255)
                             ctcGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGt                                                                                                                    <  7:221629/147‑1 (MQ=255)
                                    cTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCa                                                                                                            >  4:135024/1‑148 (MQ=255)
                                       aTGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATg                                                                                                        >  4:311670/1‑149 (MQ=255)
                                          aaGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGtt                                                                                                     <  2:75300/149‑1 (MQ=255)
                                            gTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGtttt                                                                                                   >  4:209300/1‑149 (MQ=255)
                                                gATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTatat                                                                                               >  3:377504/1‑149 (MQ=255)
                                                  ttttATTATTTTGAAGCGATTAGGCGTTGTCATGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATAttt                                                                                             <  8:138241/149‑1 (MQ=255)
                                                       ttattTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTAt                                                                                         <  4:127690/148‑1 (MQ=255)
                                                         attTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTc                                                                                       <  4:377504/148‑1 (MQ=255)
                                                         attTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCg                                                                                      >  3:154270/1‑149 (MQ=255)
                                                            ttGAAGCGATGAGGCGTAGTCAAGGATGGTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGaa                                                                                   >  1:212639/1‑149 (MQ=255)
                                                                  cGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGtt                                                                             <  3:311670/149‑1 (MQ=255)
                                                                   gATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGCACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTc                                                                            <  3:268923/149‑1 (MQ=255)
                                                                    aTGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTc                                                                                       <  2:322648/137‑1 (MQ=255)
                                                                    aTGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTc                                                                                       >  1:322648/1‑137 (MQ=255)
                                                                     tGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGGAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGc                                                                                <  1:52651/143‑1 (MQ=255)
                                                                     tGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGGAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGc                                                                                >  2:52651/1‑143 (MQ=255)
                                                                           gTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCtt                                                                    >  4:367951/1‑149 (MQ=255)
                                                                           gTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTAAGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGac                                                                                                                                        >  4:171268/1‑80 (MQ=255)
                                                                               tCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCTGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGtt                                                                                                     <  3:330631/112‑1 (MQ=255)
                                                                               tCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCGGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGtt                                                                                                     >  4:330631/1‑112 (MQ=255)
                                                                                 aaGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggc                                                              <  2:212639/149‑1 (MQ=255)
                                                                                 aaGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCgg                                                               <  3:135024/148‑1 (MQ=255)
                                                                                       gCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggcagggca                                                        >  7:6891/1‑149 (MQ=255)
                                                                                         tttCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCgg                                                                                     <  3:231225/118‑1 (MQ=255)
                                                                                         tttCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCgg                                                                                     >  4:231225/1‑118 (MQ=255)
                                                                                         tttCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggcagggcag                                                       >  2:261986/1‑148 (MQ=255)
                                                                                          ttCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGa                                                      >  2:267694/1‑148 (MQ=255)
                                                                                                gtTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAAAAcg                                               >  1:161201/1‑149 (MQ=255)
                                                                                                            tACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCat                                   <  5:146827/149‑1 (MQ=255)
                                                                                                              cGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCCTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCata                                  <  1:267694/148‑1 (MQ=255)
                                                                                                               gATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggca                                                             >  7:219077/1‑120 (MQ=255)
                                                                                                               gATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggca                                                             <  8:219077/120‑1 (MQ=255)
                                                                                                                    aaTCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCg                                                                                      >  5:8106/1‑90 (MQ=255)
                                                                                                                    aaTCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCg                                                                                      <  6:8106/90‑1 (MQ=255)
                                                                                                                        aaCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCgg                                                                                     <  3:140983/87‑1 (MQ=255)
                                                                                                                        aaCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATAtttcatttca                        <  2:161201/148‑1 (MQ=255)
                                                                                                                         aCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGa                                                      <  4:35379/117‑1 (MQ=255)
                                                                                                                         aCAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGa                                                      >  3:35379/1‑117 (MQ=255)
                                                                                                                               tcaagttGAAGCAACTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATatttca                             <  1:39482/131‑1 (MQ=255)
                                                                                                                               tcaagttGAAGCAACTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATatttca                             >  2:39482/6‑136 (MQ=255)
                                                                                                                                   gttGAAGCAACTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATaaa            >  4:313785/2‑149 (MQ=255)
                                                                                                                                   gttGAAGCAACTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATaa             >  4:172038/2‑148 (MQ=255)
                                                                                                                                   cTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATaa             >  6:151654/1‑148 (MQ=255)
                                                                                                                                   cTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATaa             <  5:151654/148‑1 (MQ=255)
                                                                                                                                     tGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGACAAATGGTTTTATATTTTTATTCGGAATCCGTTCGCTGCATTTTCGGCAGGGAAGAAAAAACGTCGTTTTTTCATATTTCATTTCAGTGGCCAATAAAca          >  8:281732/1‑148 (MQ=255)
                                                                                                                                      gAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggcagggc                                                         <  7:22640/101‑1 (MQ=255)
                                                                                                                                      gAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCggcagggc                                                         >  8:22640/1‑101 (MQ=255)
                                                                                                                                       aaGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTAAGTGGCACATAATCGCt        >  8:352733/1‑149 (MQ=255)
                                                                                                                                        aGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCtt       <  4:154270/149‑1 (MQ=255)
                                                                                                                                          cAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCTTTg     <  2:196912/149‑1 (MQ=255)
                                                                                                                                            gTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATatt                                >  4:93948/1‑120 (MQ=255)
                                                                                                                                            gTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATatt                                <  3:93948/120‑1 (MQ=255)
                                                                                                                                            gTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCTTTgg    <  6:132520/148‑1 (MQ=255)
                                                                                                                                            gTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCTTTgg    >  5:132520/1‑148 (MQ=255)
                                                                                                                                             ttGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCTTTGGCa  <  1:261986/149‑1 (MQ=255)
                                                                                                                                             ttGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCAATAAACGCTTTgg    >  8:290235/1‑147 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
AGTGCTGGTCTGGCTGGAAGATAGGTCTTCTCGTCGCTTATGAAGTTTGATTTTATTATTTTGAAGCGATGAGGCGTAGTCAAGGATGCTTTCGTCGTTGGTGCTCAGTACGATGGAATCAACAGCGTCTGCTTGAAGCAGTTGCTTCAACTTGATTTCAATCAAGCCGAATTCGTGGCCGCCAAATGGTTTTATATTTTTATTCGGAATGCGTTCGCTGCCTTTTCGGCAGGGCAGAAAACACGTCGTTTTTTCATATTTCATTTCAGTGGCCCATAAACGCTTTGGCA  >  NC_002947/2007098‑2007387

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: