Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,912,562 G→C 5.1% C303S (TGC→TCC) ‡ PP_3452 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,912,5620GC 5.1% 261.4 / 3.4 79C303S (TGC→TCC) ‡PP_3452GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base C (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

GTGGGCGACGAAGTGCTCAAGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCAAGCACACTGGCAAGG  >  NC_002947/3912414‑3912702
                                                                                                                                                    |                                                                                                                                            
gTGGGCGACGAAGTGCTCAAGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCtgtg                                                                                                                                              >  1:89218/1‑149 (MQ=255)
 tGGGCGACGAAGTGCTCAAGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGc                                                                                                                                             <  2:79338/149‑1 (MQ=255)
   ggCGACGAAGTGCTCAAGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGc                                                                                                                                             >  1:315461/1‑147 (MQ=255)
     cgacgaAGTGCTCAAGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAgc                                                                                                                                         >  5:9164/1‑149 (MQ=255)
      gacgaAGTGCTCAAGCGCGTGGCCCTACATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAgcg                                                                                                                                        <  2:184457/149‑1 (MQ=255)
        cgaAGTGCTCAAGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGcgcg                                                                                                                                      >  7:28710/1‑149 (MQ=255)
             tGCTCAAGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTgcgc                                                                                                                                 <  8:80373/149‑1 (MQ=255)
                tCAAGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAg                                                                                                                              >  4:95024/1‑149 (MQ=255)
                  aaGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAgg                                                                                                                             <  7:102070/148‑1 (MQ=255)
                   aGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGGCACAGcgcgg                                                                                                                                     <  5:161395/139‑2 (MQ=255)
                   aGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGGCACAGcgcgg                                                                                                                                     >  6:161395/1‑138 (MQ=255)
                    gcgcGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCGGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGATGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACg                                                                                                                          >  7:54602/1‑149 (MQ=255)
                                    atCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGt                                                                                                          >  8:124931/1‑149 (MQ=255)
                                     tCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGt                                                                                                          >  4:52321/1‑148 (MQ=255)
                                     tCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGt                                                                                                          >  7:94700/1‑148 (MQ=255)
                                     tCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTc                                                                                                         >  8:257611/1‑149 (MQ=255)
                                      cAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCa                                                                                                        >  6:183147/1‑149 (MQ=255)
                                       aaGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCAc                                                                                                       <  2:89218/149‑1 (MQ=255)
                                       aaGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCAc                                                                                                       >  3:192558/1‑149 (MQ=255)
                                        aGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCAcc                                                                                                      <  7:257611/149‑1 (MQ=255)
                                        aGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCAc                                                                                                       >  8:19124/1‑148 (MQ=255)
                                         gACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCAcc                                                                                                      >  2:114884/1‑148 (MQ=255)
                                          aCAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCAGTCACCAt                                                                                                    <  7:21019/149‑1 (MQ=255)
                                            aaGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCa                                                                                                  <  5:183147/149‑1 (MQ=255)
                                              gCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAg                                                                                                 >  7:17074/1‑148 (MQ=255)
                                                 tGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCAtt                                                                                             <  6:305588/149‑1 (MQ=255)
                                                 tGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCAtt                                                                                             <  1:114884/149‑1 (MQ=255)
                                                   gCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATtgg                                                                                           <  1:121023/149‑1 (MQ=255)
                                                   gCCGATCAGGCGTTAGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATtgg                                                                                           >  8:263366/1‑149 (MQ=255)
                                                    ccGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATtggc                                                                                          <  4:311231/149‑1 (MQ=255)
                                                       aTCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATtggcctg                                                                                       <  3:95024/149‑1 (MQ=255)
                                                       aTCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATtggcctg                                                                                       <  4:192558/149‑1 (MQ=255)
                                                       aTCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATtggcctg                                                                                       <  7:263366/149‑1 (MQ=255)
                                                       aTCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATtggcctg                                                                                       <  6:21480/149‑1 (MQ=255)
                                                       aTCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATtggcctg                                                                                       <  5:73343/149‑1 (MQ=255)
                                                        tCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCa                                                                                                  <  2:108774/137‑1 (MQ=255)
                                                        tCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCa                                                                                                  >  1:108774/1‑137 (MQ=255)
                                                              gTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAgg                                                                                >  2:222715/1‑149 (MQ=255)
                                                                  gCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGGCGCTGGCGAGCACCCGGTGCCCGTCAGCATTGGCCTGGCCAAGGTCGa                                                                            >  8:237107/1‑149 (MQ=255)
                                                                    ccGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCg                                                                                                                    <  4:69118/107‑1 (MQ=255)
                                                                    ccGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCg                                                                                                                    >  3:69118/1‑107 (MQ=255)
                                                                     cGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATc                                                                                                   >  6:310736/1‑123 (MQ=255)
                                                                     cGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATc                                                                                                   <  5:310736/123‑1 (MQ=255)
                                                                            gAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGGCGACCGGGga                                                                     >  6:110201/1‑146 (MQ=255)
                                                                                 gAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCt                                                             >  4:234817/1‑149 (MQ=255)
                                                                                 gAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCt                                                             >  3:33951/1‑149 (MQ=255)
                                                                                 gAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCt                                                             >  7:121560/1‑149 (MQ=255)
                                                                                   agagTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTa                                                           <  7:124931/149‑1 (MQ=255)
                                                                                      gTTCGCGATCTTCATCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAACATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTacca                                                        <  2:298466/149‑1 (MQ=255)
                                                                                       ttCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTaccac                                                       >  5:296563/1‑149 (MQ=255)
                                                                                          gcgATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCAt                                                                                                    >  6:157383/1‑101 (MQ=255)
                                                                                          gcgATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCAt                                                                                                    <  5:157383/101‑1 (MQ=255)
                                                                                          gcgATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACcgc                                                    >  6:117063/1‑149 (MQ=255)
                                                                                                ttcttcGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCtt                                                            >  4:269215/1‑135 (MQ=255)
                                                                                                ttcttcGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCtt                                                            <  3:269215/135‑1 (MQ=255)
                                                                                                 tcttcGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTc                                             >  8:167684/1‑149 (MQ=255)
                                                                                                  cttcGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCgcgc                                                  <  4:33951/143‑1 (MQ=255)
                                                                                                  cttcGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGAGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTcc                                            >  1:43654/1‑149 (MQ=255)
                                                                                                    tcGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCcgc                                          >  2:196971/1‑149 (MQ=255)
                                                                                                     cGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCAGCGCggtggtaagg                                          <  6:9164/148‑5 (MQ=255)
                                                                                                      gAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCgcgc                                        <  2:315461/149‑1 (MQ=255)
                                                                                                       aaGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGTGAGAGAGCCCTACCACCGCGCTGGTCCGCGcc                                       <  8:17074/149‑1 (MQ=255)
                                                                                                       aaGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGcc                                       <  3:234817/149‑1 (MQ=255)
                                                                                                            cTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCa                                  <  8:121560/149‑1 (MQ=255)
                                                                                                               gATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACcgcg                                                   >  4:175127/1‑129 (MQ=255)
                                                                                                               gATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACcgcg                                                   <  3:175127/129‑1 (MQ=255)
                                                                                                               gATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAgg                                >  2:271112/1‑148 (MQ=255)
                                                                                                                     cAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATtata                         <  5:110201/149‑1 (MQ=255)
                                                                                                                          ggCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAg                    >  2:15741/1‑149 (MQ=255)
                                                                                                                          ggCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAg                    >  8:34895/1‑149 (MQ=255)
                                                                                                                              gAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCaa                <  8:94700/149‑1 (MQ=255)
                                                                                                                               agcagcGAGTGGCAGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCAAg               <  8:54602/149‑1 (MQ=255)
                                                                                                                                gcagcGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCAGGGAGAGAGCCTTACCAGCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCAAGc              >  5:276755/1‑149 (MQ=255)
                                                                                                                                 cagcGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCAAGca             >  6:36957/1‑149 (MQ=255)
                                                                                                                                 cagcGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCAAGca             <  2:43654/149‑1 (MQ=255)
                                                                                                                                       gTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCAAGCACACTgg       >  2:204681/1‑149 (MQ=255)
                                                                                                                                         gcacgcgCTGGCACAGCGCGTGCGCTAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCaa                <  7:131329/135‑1 (MQ=255)
                                                                                                                                         gcacgcgCTGGCACAGCGCGTGCGCTAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCaa                >  8:131329/4‑138 (MQ=255)
                                                                                                                                          gCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCAAGCACACTGGc      >  2:106271/1‑147 (MQ=255)
                                                                                                                                            cgcgCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCAAGCACACTGGCAAgg  <  1:222715/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GTGGGCGACGAAGTGCTCAAGCGCGTGGCCCTGCATATCAAGACAAGCATGGCCGATCAGGCGTTTGCCCGGCTGGGAGGTGAAGAGTTCGCGATCTTCTTCGAAGGGCTGGATCAGCAAGGGGCCGAGCAGCGAGTGGCCGCGCTGTGCCAGCGCGTGCGCGAGGACGCTGGCGAGCACCCGGTCACCATCAGCATTGGCCTGGCCAAGGTCGACCGGGGAGAGAGCCTTACCACCGCGCTGGTCCGCGCCGATCAGGCATTATATGAAGCCAAGCACACTGGCAAGG  >  NC_002947/3912414‑3912702

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