Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,977,125 T→C 7.7% S263S (AGT→AGC)  PP_1770 → 3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,977,1250TC 7.7% 181.2 / 4.9 66S263S (AGT→AGC) PP_17703‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base T (26/34);  new base C (2/3);  total (28/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 = 7.33e-01

CGCAATGAAAACCTCGAGATCTTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAA  >  NC_002947/1976979‑1977265
                                                                                                                                                  |                                                                                                                                            
cGCAATGAAAACCTCGAGATCTTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGa                                                                                                                                            >  6:43966/1‑149 (MQ=255)
   aaTGAAAACCTCGAGATCTTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGAtct                                                                                                                                         >  4:167446/1‑149 (MQ=255)
      gAAAACCTCGAGATCTTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGCGATcgcgt                                                                                                                                      >  3:62291/1‑148 (MQ=255)
             tCGAGATCTTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTg                                                                                                                               >  4:310559/1‑149 (MQ=255)
             tCGAGATCTTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTg                                                                                                                               >  5:130132/1‑149 (MQ=255)
                   tCTTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCAGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGcgcgac                                                                                                                         >  1:284443/1‑149 (MQ=255)
                     ttCCTGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGAAGCGATCCGACCATGTGGCATGACATATTCCTCGCCAACAGCGACGATGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGACATGcgcgacgc                                                                                                                       >  6:274593/1‑149 (MQ=255)
                       ccGGTGCGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGCTCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGcgcgacgcga                                                                                                                     >  7:341059/1‑149 (MQ=255)
                          gTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCg                                                                                                                  >  5:271059/1‑149 (MQ=255)
                          gTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCg                                                                                                                  >  3:232259/1‑149 (MQ=255)
                             cGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGCGATCGCGTCGCGCAGGGCGTCGAGATCACTg                                                                                                               <  4:62291/149‑1 (MQ=38)
                              gCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGCGATCGCGTCGCGCAGGGCGTCGAGATCACTgc                                                                                                              <  2:371460/149‑2 (MQ=38)
                                tGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCTGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAg                                                                                                            <  2:106421/149‑1 (MQ=255)
                                  cGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAg                                                                                                            <  5:43966/147‑1 (MQ=255)
                                  cGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAg                                                                                                            <  2:353414/147‑1 (MQ=255)
                                  cGGGTGGGGTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAg                                                                                                            <  5:223846/147‑1 (MQ=255)
                                   gggtgggtTCCTCGATTTAACAAGAATCGCCGTGAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAgggg                                                                                                         <  2:284443/149‑1 (MQ=255)
                                     gtgggtTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGAt                                                                                                       >  2:172114/1‑149 (MQ=255)
                                        ggtTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATgg                                                                                                     <  1:15276/148‑1 (MQ=255)
                                          ttCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCa                                                                                                  <  8:163629/149‑1 (MQ=255)
                                             cgcgATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCg                                                                                                                                       <  6:336860/109‑1 (MQ=255)
                                                      acaAGAATCGCCGGCAGCGATCCGACCATGTGGCTTGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTgg                                                                                       <  3:85293/148‑1 (MQ=255)
                                                      acaAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTgg                                                                                       <  4:206854/148‑1 (MQ=255)
                                                       caAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGgt                                                                                     <  2:258253/149‑1 (MQ=255)
                                                       caAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGgt                                                                                     >  2:292435/1‑149 (MQ=255)
                                                       caAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGgt                                                                                     <  2:156629/149‑1 (MQ=255)
                                                       caAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGgt                                                                                     >  6:49626/1‑149 (MQ=255)
                                                        taGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCTCCAACCGCGACGCTGTCCTGCGCTCGCTTTATACATTTCGCTGGGATCTCGTCGCGCTGGGAGTCGAGATCACTGa                                                                                                              <  2:54840/122‑1 (MQ=255)
                                                        aaGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGgtg                                                                                    <  3:234977/149‑1 (MQ=255)
                                                         aGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCCACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGgtgt                                                                                   <  6:271059/149‑1 (MQ=255)
                                                          gAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTa                                                                                  >  6:174059/1‑149 (MQ=255)
                                                                ccGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCAcc                                                                            >  1:348723/1‑149 (MQ=255)
                                                                  ggCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGAt                                                                                                       >  4:128620/1‑120 (MQ=255)
                                                                     aGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCa                                                                                               >  7:1411/1‑125 (MQ=255)
                                                                     aGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCa                                                                                               <  8:1411/125‑1 (MQ=255)
                                                                         aTCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACgg                                                                   <  7:94216/149‑1 (MQ=255)
                                                                           ccGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGt                                                                 >  8:236509/1‑149 (MQ=255)
                                                                                  tgtgGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCcgc                                                          <  3:167446/149‑1 (MQ=255)
                                                                                   gtgGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCcgc                                                          >  4:110791/1‑148 (MQ=255)
                                                                                   gtgGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCcgc                                                          >  3:215077/1‑148 (MQ=255)
                                                                                    tgGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGa                                                        <  1:172114/149‑1 (MQ=255)
                                                                                    tgGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGa                                                        <  3:310559/149‑1 (MQ=255)
                                                                                         tGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACAAGCTGCTGGGTGTATTCACCCGCGCACGGCTTGCCCGCGAGCAtt                                                   >  7:169013/1‑149 (MQ=255)
                                                                                           aCATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCAttt                                                  <  4:215077/148‑1 (MQ=255)
                                                                                            cATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCa                                                <  5:174059/149‑1 (MQ=255)
                                                                                            cATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCa                                                <  1:213455/149‑1 (MQ=255)
                                                                                            cATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCa                                                <  8:341059/149‑1 (MQ=255)
                                                                                            cATCTTCCTCGACAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTCATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGTTTGCCCGCGAGCATTTCa                                                <  5:274593/149‑1 (MQ=255)
                                                                                             ctCTTCCTCGCCAACCGCGACGCTGTCCTTCGCACGCTTGATACATTTCGCAGTGATCTCGCCGCCCTGCGCGTCGCGAGCGCTGATGGGGATGGGCACGAGCTGCTGGGTGTATTGACCCGCGCACGGGTTGCCCGCGATCATTTCa                                                <  6:381769/147‑1 (MQ=255)
                                                                                                ttCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTa                                             >  3:50972/1‑148 (MQ=255)
                                                                                                 tCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTaa                                            >  6:123244/1‑148 (MQ=255)
                                                                                                  ccTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTaaa                                           <  3:110791/148‑1 (MQ=255)
                                                                                                  ccTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTa                                             <  4:50810/146‑1 (MQ=255)
                                                                                                    tCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGa                                                        <  6:267156/133‑1 (MQ=255)
                                                                                                    tCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGa                                                        >  5:267156/1‑133 (MQ=255)
                                                                                                       ccAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTg                                     <  4:50972/149‑1 (MQ=255)
                                                                                                       ccAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTg                                     <  1:147724/149‑1 (MQ=255)
                                                                                                                 cGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGGCggg                           <  5:5027/149‑1 (MQ=255)
                                                                                                                  gCTGTCCTGCGCACGCTTGATACATTTCGCCGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGc                          <  4:210601/149‑1 (MQ=255)
                                                                                                                        cTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTAtgt                    <  4:232259/149‑1 (MQ=255)
                                                                                                                             cACGCTTGATAAATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGAGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACg               >  8:348186/1‑149 (MQ=255)
                                                                                                                                 cTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTAt           <  6:130132/149‑1 (MQ=255)
                                                                                                                                 cTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTAt           >  3:144061/1‑149 (MQ=255)
                                                                                                                                          tttCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCaa  >  4:334025/1‑149 (MQ=255)
                                                                                                                                            cttcAGCGATCGCGTCGCGCAGGGCGTCGAGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTAtgtg                   >  1:102165/4‑130 (MQ=37)
                                                                                                                                            cttcAGCGATCGCGTCGCGCAGGGCGTCGAGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTAtgtg                   <  2:102165/127‑1 (MQ=37)
                                                                                                                                                  |                                                                                                                                            
CGCAATGAAAACCTCGAGATCTTCCGGTACGCTGCGGGTGGGTTCCGCGATTTCACAAGAATCGCCGGCAGCGATCCGACCATGTGGCATGACATCTTCCTCGCCAACCGCGACGCTGTCCTGCGCACGCTTGATACATTTCGCAGTGATCTCGACGCCCTGCGCGACGCGATCGCTGAAGGGGATGGGCACCAGCTGCTGGGTGTATTCACCCGCGCACGGGTTGCCCGCGAGCATTTCAGTAAAATCCTGGCCCGCCGGGCCTATGTGGACGCTATGAACGCCAA  >  NC_002947/1976979‑1977265

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