Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,789,594 G→C 6.1% *63S (TGA→TCA)  PP_2441 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,789,5940GC 6.1% 207.0 / 3.0 65*63S (TGA→TCA) PP_2441hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base C (2/2);  total (33/32)
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.86e-01

CAAAGCCCGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGCCTCGATTTCC  >  NC_002947/2789447‑2789741
                                                                                                                                                   |                                                                                                                                                    
cAAAGCCCGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGa                                                                                                                                                     <  4:92254/149‑1 (MQ=255)
       cGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCa                                                                                                                                              <  1:33274/149‑1 (MQ=255)
             gcaAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTc                                                                                                                                        >  3:150191/1‑149 (MQ=255)
               aaGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCa                                                                                                                                      <  4:188737/149‑1 (MQ=255)
                 gagCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGc                                                                                                                                    >  8:1721/1‑149 (MQ=255)
                          cTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCTGAGCCGTCACGGCTcggggtc                                                                                                                                         <  6:70443/135‑7 (MQ=255)
                          cTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCTGAGCCGTCACGGCTcggggtc                                                                                                                                         >  5:70443/1‑129 (MQ=255)
                              aCTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGgcgcg                                                                                                                      >  3:77973/1‑149 (MQ=255)
                                  gggCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTc                                                                                                                  >  4:49649/1‑149 (MQ=255)
                                  gggCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTc                                                                                                                  >  5:331988/1‑149 (MQ=255)
                                   ggCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTc                                                                                                                  >  2:97074/1‑148 (MQ=255)
                                    gCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCa                                                                                                                <  8:131614/149‑1 (MQ=255)
                                        ccTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAg                                                                                                                                     <  4:84709/125‑1 (MQ=255)
                                        ccTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAg                                                                                                                                     >  3:84709/1‑125 (MQ=255)
                                           ggAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGcgccg                                                                                                          >  2:298085/1‑148 (MQ=255)
                                             aaaTCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCaa                                                                                                       >  5:315251/1‑149 (MQ=255)
                                                    cGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCgggc                                                                                                <  3:11326/149‑1 (MQ=255)
                                                     gCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCGGCAATTTCCCCGACACCGAGCCGTGACGGCTCAGACTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGAAATGCgggct                                                                                               >  6:326855/1‑149 (MQ=255)
                                                          cGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCgggcttcgag                                                                                          <  4:281005/149‑1 (MQ=255)
                                                           ggTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCgggcttcgagg                                                                                         <  4:150191/149‑1 (MQ=255)
                                                                cgcTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCgggcttcgagggcttc                                                                                    >  5:178908/1‑149 (MQ=255)
                                                                 gcTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCgggcttcgagggcttcg                                                                                   <  1:320571/149‑1 (MQ=255)
                                                                     tgtTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCTCCGCAATGCGGGCTTCGAGGGCTTCGATg                                                                                <  4:58638/148‑1 (MQ=255)
                                                                     tgtTCCTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGAt                                                                                 <  4:210477/147‑1 (MQ=255)
                                                                             gACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGtcttct                                                                       <  2:217492/149‑1 (MQ=255)
                                                                             gACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGtcttct                                                                       <  7:89004/149‑1 (MQ=255)
                                                                                   cgcCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGacaca                                                                 >  2:107748/1‑149 (MQ=255)
                                                                                       cTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTgc                                                             >  6:311412/1‑149 (MQ=255)
                                                                                       cTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTgc                                                             <  6:178908/149‑1 (MQ=255)
                                                                                             aCCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCAc                                                       <  3:49649/149‑1 (MQ=255)
                                                                                             aCCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCAc                                                       >  3:101254/1‑149 (MQ=255)
                                                                                              ccGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTCACgg                                                                                                                                                  <  1:288935/58‑1 (MQ=255)
                                                                                              ccGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTCACgg                                                                                                                                                  >  2:288935/1‑58 (MQ=255)
                                                                                               cGGGAGGTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGAGCCCGAGCCGTGTCGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTg                                                     <  8:64481/149‑1 (MQ=255)
                                                                                                 ggAGCTGGCCGAGCGCAAGAGCTGCGGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCATGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGt                                                   <  4:77973/149‑1 (MQ=255)
                                                                                                 ggAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGt                                                   <  6:315251/149‑1 (MQ=255)
                                                                                                     cTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGt                                                   <  7:248174/145‑1 (MQ=255)
                                                                                                     cTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGt                                                   >  8:248174/1‑145 (MQ=255)
                                                                                                     cTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTgcgc                                               <  1:107748/149‑1 (MQ=255)
                                                                                                        gCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGATGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTgcgcgcg                                            >  4:259663/1‑149 (MQ=255)
                                                                                                          cGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGtcttct                                                                       <  7:240659/120‑1 (MQ=255)
                                                                                                          cGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGtcttct                                                                       >  8:240659/1‑120 (MQ=255)
                                                                                                          cGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGcgcc                                          <  1:97074/149‑1 (MQ=255)
                                                                                                          cGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGcgcc                                          >  1:139954/1‑149 (MQ=255)
                                                                                                           gAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCCCTGGTGCGCGcgccg                                         >  5:48240/1‑149 (MQ=255)
                                                                                                             gcgcAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCgccgcc                                       >  2:66491/1‑149 (MQ=255)
                                                                                                               gcAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCtt                                     >  8:254806/1‑149 (MQ=255)
                                                                                                               gcAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCAtt                                     >  6:58317/1‑149 (MQ=255)
                                                                                                                cAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTc                                    <  6:331988/149‑1 (MQ=255)
                                                                                                                    agCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCctgc                                >  7:240121/1‑149 (MQ=255)
                                                                                                                         ccGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCg                           >  3:136876/1‑149 (MQ=255)
                                                                                                                         ccGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCg                           >  5:209650/1‑149 (MQ=255)
                                                                                                                             aaTTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCa                       <  6:209650/149‑1 (MQ=255)
                                                                                                                              aTTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAg                      >  5:125850/1‑149 (MQ=255)
                                                                                                                                ttCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGc                     >  1:163845/1‑148 (MQ=255)
                                                                                                                                    ccgaccccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATg                                                                                >  6:102479/1‑85 (MQ=255)
                                                                                                                                    ccgaccccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATg                                                                                <  5:102479/85‑1 (MQ=255)
                                                                                                                                     cgaccccgaGCCGTGACGGCTCAGCCTTCCAGCTCATGCCCGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGaa                >  1:220823/1‑149 (MQ=255)
                                                                                                                                       accccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCg             <  6:125850/149‑1 (MQ=255)
                                                                                                                                        ccccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCgag                  <  4:273440/143‑1 (MQ=255)
                                                                                                                                        ccccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCgag                  >  3:273440/1‑143 (MQ=255)
                                                                                                                                         cccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGcc           <  1:298085/149‑1 (MQ=255)
                                                                                                                                         cccgaGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGcc           >  6:270104/1‑149 (MQ=255)
                                                                                                                                               ccGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGCCTCGAtt     <  4:136876/149‑1 (MQ=255)
                                                                                                                                                  tGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGCCTCGATTTcc  <  2:139954/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAAAGCCCGGGCAGCAAGAGCAGGACCTGCACTGGGGCTACCTGGAAATCTACGCCGACGGTCGCGCTGTGTTCGTCGACCAGCGCCCTACTGACCGGGAGCTGGCCGAGCGCAAGAGCTGCCGCAATTTCCCCGACCCCGAGCCGTGACGGCTCAGCCTTCCAGCTCAGCCAGGCGCGCTTCCAGCGCCGCAATGCGGGCTTCGAGGGCTTCGATGCGGTCTTCTGACACACTGCCGCCACTGGTGCGCGCGCCGCCTTCCTGCTGCCGTGCAGCGAGAATCGCCTCGATTTCC  >  NC_002947/2789447‑2789741

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: