Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,367,361 G→C 5.6% E223Q (GAG→CAG)  PP_2964 → transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,367,3610GC 5.6% 218.9 / 4.0 71E223Q (GAG→CAG) PP_2964transposase
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base C (2/2);  total (37/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 = 3.50e-01

GAAGTTACCTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCTCGCTGCC  >  NC_002947/3367214‑3367506
                                                                                                                                                   |                                                                                                                                                 
gcaGTTACCTTCAACCTACTCGGGACTCGAGCGCCAATCGAGGTGATCTACCAGCTGAAGGAGCTGCTCGACACCGACTTCGCCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGa                                                                                                                                                  <  4:284263/147‑1 (MQ=255)
   gTTACCTTCAACCTAGTCGGGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACCCCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGt                                                                                                                                               >  2:242292/1‑149 (MQ=255)
       ccTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTCCAGGTGCtt                                                                                                                                           <  8:172082/149‑2 (MQ=255)
        cTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGa                                                                                                                                          >  8:139195/1‑149 (MQ=255)
        cTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTCCAGGTGCttg                                                                                                                                          <  5:142424/149‑3 (MQ=255)
          tCAACCTAGTCGTGACTCGAGCGACAATCGAGGTGATCTACAAGCTGAACGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACATGGAGCTCCTGa                                                                                                                                          >  8:48488/1‑147 (MQ=255)
           cAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATc                                                                                                                                        <  4:218548/148‑1 (MQ=255)
           cAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATc                                                                                                                                        <  3:290277/148‑1 (MQ=255)
                        aCTCGAGCTCCAATCGAGTTTATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCTCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCt                                                                                                                           <  6:134325/148‑1 (MQ=255)
                            gAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGc                                                                                                                        >  8:301388/1‑147 (MQ=255)
                            gAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCt                                                                                                                       <  1:249575/148‑1 (MQ=255)
                                  cAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACCCCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGc                                                                                                                >  4:331500/1‑149 (MQ=255)
                                       gatgTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGTCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCgg                                                                                                           <  7:48488/146‑1 (MQ=255)
                                           tGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGAc                                                                                                       <  8:238577/149‑1 (MQ=255)
                                                tACAAGCTGAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAg                                                                                                 >  6:274609/1‑149 (MQ=255)
                                                  cAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACCGGCACCCAAGAAGCCATCTCGCTAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTg                                                                                                    <  8:272344/145‑1 (MQ=255)
                                                  cAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTTGCCTGGCTAGAGAGCTTCGGGGACCTGGAg                                                                                                 <  6:214207/148‑1 (MQ=255)
                                                  cAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAg                                                                                                 <  1:242292/148‑1 (MQ=255)
                                                        gAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCt                                                                                          >  1:49557/1‑149 (MQ=255)
                                                        gAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACACAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACATGGCCTGGCTAGAGAGATTCGGGGACGAGTACGCCATCt                                                                                          >  8:315209/1‑149 (MQ=255)
                                                           ggAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTg                                                                                       <  2:740/149‑1 (MQ=255)
                                                             aGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATGGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAg                                                                                     <  7:315209/149‑1 (MQ=255)
                                                             aGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAg                                                                                     <  7:341742/149‑1 (MQ=255)
                                                             aGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCAACTCCCCAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAg                                                                                     <  1:320812/149‑1 (MQ=255)
                                                              gctgctCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAgg                                                                                    <  2:49557/149‑1 (MQ=255)
                                                                tgctCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGc                                                                                                                        >  3:88848/1‑111 (MQ=255)
                                                                tgctCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGc                                                                                                                        <  4:88848/111‑1 (MQ=255)
                                                                   tCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGAtc                                                                               >  1:244177/1‑149 (MQ=255)
                                                                       cacCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGAtcgtc                                                                            >  8:149673/1‑148 (MQ=255)
                                                                        acCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGa                                                                          >  4:84395/1‑149 (MQ=255)
                                                                         ccGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGAt                                                                         <  2:244177/149‑1 (MQ=255)
                                                                          cGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATg                                                                        >  6:338776/1‑149 (MQ=255)
                                                                          cGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATg                                                                        >  5:297444/1‑149 (MQ=255)
                                                                          cGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATg                                                                        >  3:128293/1‑149 (MQ=255)
                                                                           gACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGc                                                                       <  5:274609/149‑1 (MQ=255)
                                                                              ttCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCCCCCAAGACGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGcc                                                                      >  8:367445/1‑147 (MQ=255)
                                                                              ttCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAg                                                                                                 <  5:116897/120‑1 (MQ=255)
                                                                              ttCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAg                                                                                                 >  6:116897/1‑120 (MQ=255)
                                                                              ttCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGcc                                                                      >  1:141303/1‑147 (MQ=255)
                                                                                cggccggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGCTGCTGATCCACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGcc                                                                      >  1:173942/1‑145 (MQ=255)
                                                                                 ggccggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGa                                                                  <  7:301388/148‑1 (MQ=255)
                                                                                  gccggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGTGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGAtg                                                                >  6:284002/1‑148 (MQ=255)
                                                                                   ccggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAg                                                                                     <  7:339208/127‑1 (MQ=255)
                                                                                   ccggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAg                                                                                     >  8:339208/1‑127 (MQ=255)
                                                                                   ccggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTc                                                               <  7:219606/149‑1 (MQ=255)
                                                                                     ggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGc                                                             <  4:128293/149‑1 (MQ=255)
                                                                                     ggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGc                                                             <  6:297444/149‑1 (MQ=255)
                                                                                        gcCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCAc                                                          >  7:135757/1‑149 (MQ=255)
                                                                                           agagCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACttt                                                       >  3:349980/1‑149 (MQ=255)
                                                                                              gCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTtgct                                                    >  1:67812/1‑149 (MQ=255)
                                                                                                gATTTCAATGGGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGTGACCTGGAGGCCATCTTTGAGGGGTTCGTCGATGCCAAGATTCGCCACTTtgctgc                                                  <  5:332021/149‑1 (MQ=255)
                                                                                                      aaTGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGc                                            >  4:245995/1‑149 (MQ=255)
                                                                                                       aTGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGTCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGTCAAGATTCGCCACTTTGCTGCCCAAGcc                                           >  5:342683/1‑149 (MQ=255)
                                                                                                        tGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCcg                                          >  4:142436/1‑149 (MQ=255)
                                                                                                           gcTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCgcgg                                       <  2:67812/149‑1 (MQ=255)
                                                                                                              gAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGTGATCGTCGATGCCAAGATTCGCAACTTTGCTGCCAAAGCCgaggcgt                                    >  4:254136/1‑148 (MQ=255)
                                                                                                                   agtcACCCAATAAGCCATCTCTCATGCACCTGGACGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGTCCATCTTTGAGGGGATCTTGGATGCCAAGATTCGCCACTATGCTGCCCAAGCAGCGGCGGCGGAc                               <  2:56099/146‑1 (MQ=255)
                                                                                                                   aGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGAc                               >  4:193922/1‑149 (MQ=255)
                                                                                                                                  cATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAg                <  3:245995/149‑1 (MQ=255)
                                                                                                                                  cATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAg                <  7:367445/149‑1 (MQ=255)
                                                                                                                                  cATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAg                <  7:139195/149‑1 (MQ=255)
                                                                                                                                  cATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAg                <  5:55439/149‑1 (MQ=255)
                                                                                                                                   aTCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAg                >  3:63228/1‑148 (MQ=255)
                                                                                                                                   aTCTCGCAAGCACCTGGAGGAGCTGATCGACCAGTAGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAgg               >  6:347524/1‑149 (MQ=255)
                                                                                                                                     ctcGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGcagtcctt              >  2:189775/1‑141 (MQ=255)
                                                                                                                                      tcGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACt            >  5:372739/1‑149 (MQ=255)
                                                                                                                                       gatcaGCACCTCCAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTg           >  5:154669/5‑149 (MQ=255)
                                                                                                                                         tcaGCACCTCCAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCt         >  3:227395/3‑149 (MQ=255)
                                                                                                                                           aGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCACAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTCCTCg       >  8:177473/1‑149 (MQ=255)
                                                                                                                                             cACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCTCGCt     <  3:193922/149‑1 (MQ=255)
                                                                                                                                                cTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCTCGCTGcc  >  1:157439/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAAGTTACCTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCTCGCTGCC  >  NC_002947/3367214‑3367506

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: