Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,790,620 A→T 7.1% intergenic (‑24/+33) PP_5073 ← / ← hemE hypothetical protein/uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,790,6200AT 7.1% 239.6 / 5.5 81intergenic (‑24/+33)PP_5073/hemEhypothetical protein/uroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base T (4/2);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGC  >  NC_002947/5790479‑5790765
                                                                                                                                             |                                                                                                                                                   
cAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGc                                                                                                                                              >  5:88475/1‑149 (MQ=255)
  aaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCg                                                                                                                                            <  6:88475/149‑1 (MQ=255)
      cgcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGtt                                                                                                                                        <  3:113407/149‑1 (MQ=255)
       gcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGtt                                                                                                                                        >  5:253731/1‑148 (MQ=255)
        ccAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTa                                                                                                                                      <  8:60110/149‑1 (MQ=255)
         cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTAc                                                                                                                                     >  2:336812/1‑149 (MQ=255)
                   cAAGGGGGCGAGGATGAGGGCGAGATAACTGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagtgc                                                                                                                           <  7:331013/149‑1 (MQ=255)
                         ggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAg                                                                                                                     >  3:8224/1‑149 (MQ=255)
                          gCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGt                                                                                                                    >  5:39744/1‑149 (MQ=255)
                                aTGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCa                                                                                                                                                     >  6:167338/1‑110 (MQ=255)
                                 tGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggt                                                                                                             <  1:336812/149‑1 (MQ=255)
                                 tGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGACGGGTTTACTGTGCAGTGCATCAAGTCACtggt                                                                                                             >  8:208845/1‑149 (MQ=255)
                                   aGGGCGAGATAACGGATTTTCATGATAATCGGTTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCTGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtgg                                                                                                           <  4:25817/149‑1 (MQ=255)
                                       cGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACt                                                                                                       >  2:100976/1‑149 (MQ=255)
                                         agaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagtgc                                                                                                                           >  3:44104/1‑127 (MQ=255)
                                         agaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagtgc                                                                                                                           <  4:44104/127‑1 (MQ=255)
                                         agaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgc                                                                                                     <  7:208845/149‑1 (MQ=255)
                                           aTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgcgc                                                                                                   >  3:324127/1‑149 (MQ=255)
                                               cGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGATGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAg                                                                                                                     <  5:232423/127‑1 (MQ=255)
                                               cGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTAGATGGGACCAGATTAGCGAGGATCGCTGAATTAAAACTTAACGTACAAAAAACACGGCAACGAAGCCGGGTTTACTCTGAAGTGCATCAAg                                                                                                                     >  6:232423/1‑127 (MQ=255)
                                                ggATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGa                                                                                              >  8:316567/1‑147 (MQ=255)
                                                     ttCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTc                                                                                        <  6:39744/148‑1 (MQ=255)
                                                      tCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgtg                                                                                                                                 >  4:248244/1‑108 (MQ=255)
                                                      tCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgtg                                                                                                                                 <  3:248244/108‑1 (MQ=255)
                                                      tCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATc                                                                                                                        >  4:183149/1‑117 (MQ=255)
                                                      tCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATc                                                                                                                        <  3:183149/117‑1 (MQ=255)
                                                        aTGATATTCGCGTCCTGTTTCGGTTTGGTTGGGGCCACATTAGCGGCGATCGCTGAAGTCAACGGTAACTTACAAAAAAACCGGCATCGAAGCCGGGTTTAGTGTGCAGTGCCTGAAGGCACTGGTGGTCCTGCGC‑‑CGACAGTTCGTg                                                                                     >  6:140333/1‑148 (MQ=255)
                                                        aTGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAATAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGc                                                                                    <  1:106010/149‑1 (MQ=255)
                                                         tGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCa                                                                                   <  6:253731/149‑1 (MQ=255)
                                                          gATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtgg                                                                                                           >  7:57812/1‑126 (MQ=255)
                                                          gATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtgg                                                                                                           <  8:57812/126‑1 (MQ=255)
                                                              aTCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGc                                                                               <  2:63096/148‑1 (MQ=255)
                                                               tCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGt                                                                             <  4:324127/149‑1 (MQ=255)
                                                                 gggTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGtt                                                                            <  4:8224/148‑1 (MQ=255)
                                                                       tGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAAt                                                                     >  5:121534/1‑149 (MQ=255)
                                                                             ggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCAccc                                                               >  5:88192/1‑149 (MQ=255)
                                                                                   tttGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGaa                                                                      <  8:178617/134‑1 (MQ=255)
                                                                                   gATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTc                                                                                        >  6:144395/1‑118 (MQ=255)
                                                                                   gATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTc                                                                                        <  5:144395/118‑1 (MQ=255)
                                                                                   gATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGaa                                                                      >  7:178617/1‑136 (MQ=255)
                                                                                     tGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtg                                                                                                               >  3:283243/1‑95 (MQ=255)
                                                                                     tGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtg                                                                                                               <  4:283243/95‑1 (MQ=255)
                                                                                     tGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTc                                                       <  6:88192/149‑1 (MQ=255)
                                                                                         aCCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTAc                                                                                                        >  2:208854/1‑98 (MQ=255)
                                                                                         aCCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCggg                                                   <  6:121534/149‑1 (MQ=255)
                                                                                           cAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGa                                                                                              <  8:272892/104‑1 (MQ=255)
                                                                                           cAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGa                                                                                              >  7:272892/1‑104 (MQ=255)
                                                                                                  gCGACGATCGCTGAATTCAAACTTAACTTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTcc                                          >  1:10565/1‑149 (MQ=255)
                                                                                                               aaTTCAAACTTAACGTACATAATATCCGGCTTCGAAGGCGGGTTTAGTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGGTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGt                             <  5:140333/149‑1 (MQ=255)
                                                                                                               aaTTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGcc                               <  1:100976/147‑1 (MQ=255)
                                                                                                                 ttCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCAccc                                                               <  8:295998/113‑1 (MQ=255)
                                                                                                                 ttCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCAccc                                                               >  7:295998/1‑113 (MQ=255)
                                                                                                                 ttCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTgg                           >  1:227334/1‑149 (MQ=255)
                                                                                                                   cAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGc                          >  7:28678/1‑148 (MQ=255)
                                                                                                                    aaaCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGt                                     >  4:56102/1‑136 (MQ=255)
                                                                                                                    aaaCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGt                                     <  3:56102/136‑1 (MQ=255)
                                                                                                                    aaaCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGcc                         >  2:349611/1‑148 (MQ=255)
                                                                                                                    aaaCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGcc                         >  1:187362/1‑148 (MQ=255)
                                                                                                                      aCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAg                      >  5:316717/1‑149 (MQ=255)
                                                                                                                          aaCGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCaa                       <  5:22953/144‑1 (MQ=255)
                                                                                                                          aaCGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCaa                       >  6:22953/1‑144 (MQ=255)
                                                                                                                           aCGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGt                                     >  6:61779/1‑129 (MQ=255)
                                                                                                                           aCGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGt                                     <  5:61779/129‑1 (MQ=255)
                                                                                                                              tACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAaca              >  4:291103/1‑149 (MQ=255)
                                                                                                                                cAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGCCCAAGGTTGAAcaca            >  8:24862/1‑149 (MQ=255)
                                                                                                                                 aaaaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACAt           <  6:316717/149‑1 (MQ=255)
                                                                                                                                 aaaaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACAt           <  2:10565/149‑1 (MQ=255)
                                                                                                                                 aaaaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACAt           >  7:59916/1‑149 (MQ=255)
                                                                                                                                   taaaCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCTGGGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACAt           >  1:274542/2‑149 (MQ=255)
                                                                                                                                   aaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGaa                <  1:234723/142‑1 (MQ=255)
                                                                                                                                   aaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGaa                >  2:234723/1‑142 (MQ=255)
                                                                                                                                   aaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAac               >  1:268248/1‑143 (MQ=255)
                                                                                                                                   aaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAac               <  2:268248/143‑1 (MQ=255)
                                                                                                                                   aaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATgg         >  8:257870/1‑149 (MQ=255)
                                                                                                                                   aaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTCGAACACATgg         >  3:151989/1‑149 (MQ=255)
                                                                                                                                    aaaCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTtgat                                                                         <  3:261400/84‑1 (MQ=255)
                                                                                                                                    aaaCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTtgat                                                                         >  4:261400/1‑84 (MQ=255)
                                                                                                                                      aCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGc        >  8:194449/1‑147 (MQ=255)
                                                                                                                                      aCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGc        >  3:312211/1‑147 (MQ=255)
                                                                                                                                      aCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGc        >  4:303277/1‑147 (MQ=255)
                                                                                                                                        ccGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGt    <  8:28678/149‑1 (MQ=255)
                                                                                                                                          ggCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGc  <  8:59916/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGC‑‑CGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGC  >  NC_002947/5790479‑5790765

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