Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,034,966 G→A 5.3% A335T (GCG→ACG)  PP_0892 → phospholipase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,034,9660GA 5.3% 282.9 / 2.2 91A335T (GCG→ACG) PP_0892phospholipase family protein
Reads supporting (aligned to +/- strand):  ref base G (39/47);  new base A (3/2);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

ACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACC  >  NC_002947/1034818‑1035110
                                                                                                                                                    |                                                                                                                                                   
aCGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAg                                                                                                                                                     >  3:257053/1‑149 (MQ=255)
   gATTTTCAAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAgcgg                                                                                                                                                  <  7:288270/149‑1 (MQ=255)
          gAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTTGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTAGACCCGCCGTTGACAgcgg                                                                                                                                                  >  3:351475/1‑142 (MQ=255)
          gAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTTGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTAGACCCGCCGTTGACAgcgg                                                                                                                                                  <  4:351475/142‑1 (MQ=255)
             tACCAGGCGTGTTTCCTGCTCGTGTAGTTGGCTGTGGTAGTCGACTGGGGCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGGGCTTCTACCTTGATGCGAACTTCGAAGCGCTCGACCCGCGGTTGACAGCGGCGGtggtgg                                                                                                                                         <  6:8953/148‑1 (MQ=255)
                cAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAg                                                                                                                                     >  8:18038/1‑149 (MQ=255)
                cAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGCTGTCAAcggcgggtcgagc                                                                                                                                        >  1:291053/1‑139 (MQ=255)
                cAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGCTGTCAAcggcgggtcgagc                                                                                                                                        <  2:291053/146‑8 (MQ=255)
                  gCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTTGCCAGc                                                                                                                                    >  2:313520/1‑148 (MQ=255)
                  gCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCt                                                                                                                                   >  5:357105/1‑149 (MQ=255)
                      tGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGa                                                                                                                               <  5:142762/149‑1 (MQ=255)
                        cTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAgc                                                                                                                             >  6:70889/1‑149 (MQ=255)
                           ccTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAgcgc                                                                                                                           >  2:278523/1‑148 (MQ=255)
                           ccTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAgcgc                                                                                                                           <  4:167065/148‑1 (MQ=255)
                                aCCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTc                                                                                                                     >  3:295501/1‑149 (MQ=255)
                                     aaGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGccc                                                                                                                 <  7:18038/148‑1 (MQ=255)
                                       gATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCaga                                                                                                              <  1:187695/149‑1 (MQ=255)
                                       gATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCaga                                                                                                              <  3:335750/149‑1 (MQ=255)
                                        aTGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCagag                                                                                                             >  5:369056/1‑149 (MQ=255)
                                               tgGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCTGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCgaag‑‑aa                                                                                                    <  3:37116/149‑1 (MQ=255)
                                                  tAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTc                                                                                                 >  3:330490/1‑149 (MQ=255)
                                                    gacgacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGa                                                                                               >  4:275250/1‑149 (MQ=255)
                                                    gacgacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACCTCGCCCAGAGCGAAG‑‑AAGTCGa                                                                                               >  5:162706/1‑149 (MQ=255)
                                                     acgacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGAc                                                                                              <  4:257053/149‑1 (MQ=255)
                                                     acgacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGAc                                                                                              >  2:2168/1‑149 (MQ=255)
                                                      cgacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGATCTTGCACTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACTGCGGTGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGAcc                                                                                             <  1:286817/149‑1 (MQ=255)
                                                       gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCgaa                                                                                                         >  2:288535/1‑138 (MQ=255)
                                                       gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCgaa                                                                                                         <  1:288535/138‑1 (MQ=255)
                                                       gacTGGGTCAGCGTCGGTTCGTACAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑t                                                                                           >  8:79716/1‑149 (MQ=255)
                                                            ggTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGCTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGc                                                                                        <  7:201488/147‑1 (MQ=255)
                                                              tCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCAc                                                                                    >  7:106293/1‑149 (MQ=255)
                                                                aGCGGCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTg                                                                                  <  3:213896/149‑1 (MQ=255)
                                                                       gTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTg                                                                                  <  6:146700/142‑1 (MQ=255)
                                                                       gTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTg                                                                                  >  5:146700/1‑142 (MQ=255)
                                                                       gTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATgcg                                                                           >  8:4970/1‑149 (MQ=255)
                                                                           gTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGcc                                                                       >  8:336569/1‑149 (MQ=255)
                                                                           gTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGcc                                                                       >  1:58981/1‑149 (MQ=255)
                                                                           gTGCAACTTCGACCACTGGAACTTGCGCTTCAACC‑TGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGccc                                                                      <  5:273843/149‑1 (MQ=255)
                                                                             gCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCaga                                                                                                              >  6:237170/1‑111 (MQ=255)
                                                                             gCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCaga                                                                                                              <  5:237170/111‑1 (MQ=255)
                                                                                      aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGgcg                                                             >  5:49952/1‑148 (MQ=255)
                                                                                        cACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAg                                                          <  6:369056/149‑1 (MQ=255)
                                                                                         aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  7:4970/149‑1 (MQ=255)
                                                                                         aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  4:330490/149‑1 (MQ=255)
                                                                                         aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  4:301365/149‑1 (MQ=255)
                                                                                         aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  6:38585/149‑1 (MQ=255)
                                                                                         aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  8:112534/149‑1 (MQ=255)
                                                                                          cTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAggg                                                        >  4:103906/1‑149 (MQ=255)
                                                                                          cTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAggg                                                        >  3:172688/1‑149 (MQ=255)
                                                                                             gAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTg                                                      >  6:282167/1‑148 (MQ=255)
                                                                                             gAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTCGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTg                                                      >  7:30644/1‑148 (MQ=255)
                                                                                              aaCTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGaa                                                    <  8:33314/149‑1 (MQ=255)
                                                                                              aaCTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGaa                                                    <  7:79716/149‑1 (MQ=255)
                                                                                               aCTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAg                                                   <  1:278523/149‑1 (MQ=255)
                                                                                                    cgcTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAAc                                               <  6:357105/148‑1 (MQ=255)
                                                                                                       ttCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCCCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCAAGCGCGCCCTTGGTGGAGCAGCGTGAAGAAACGGa                                            >  2:290421/1‑148 (MQ=255)
                                                                                                         cAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATc                                          <  5:282167/148‑1 (MQ=255)
                                                                                                         cAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCt                                         >  2:169676/1‑149 (MQ=255)
                                                                                                          aaCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTg                                        <  5:180816/149‑1 (MQ=255)
                                                                                                          aaCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTg                                        >  6:180816/1‑149 (MQ=255)
                                                                                                          aaCCTTGAAGCCAACATCGAAGCCCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTg                                        <  3:164470/149‑1 (MQ=255)
                                                                                                           aCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATc                                          <  3:225731/146‑1 (MQ=255)
                                                                                                           aCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATc                                          >  4:225731/1‑146 (MQ=255)
                                                                                                            ccTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGCCTAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTg                                        <  1:290421/149‑1 (MQ=255)
                                                                                                             cTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTgggg                                     >  3:50915/1‑149 (MQ=255)
                                                                                                              ttGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGggc                                    >  4:293468/1‑149 (MQ=255)
                                                                                                                 aaGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGggctgg                                 <  8:30644/149‑1 (MQ=255)
                                                                                                                     cAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTtggtgg                                                               <  6:202214/115‑1 (MQ=255)
                                                                                                                     cAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTtggtgg                                                               >  5:202214/1‑115 (MQ=255)
                                                                                                                     cAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTgg                                       <  6:23407/139‑1 (MQ=255)
                                                                                                                     cAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTgg                                       >  5:23407/1‑139 (MQ=255)
                                                                                                                      aaCATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGACCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGTGCTGGCTCGa                            >  6:211863/1‑149 (MQ=255)
                                                                                                                       aCATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGAt                           <  4:50915/149‑1 (MQ=255)
                                                                                                                        cATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATc                          <  7:336569/149‑1 (MQ=255)
                                                                                                                           cGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGc                       >  7:183250/1‑149 (MQ=255)
                                                                                                                             aaGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGTTCGATCGGCtg                     <  3:62402/149‑1 (MQ=255)
                                                                                                                             aaGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtg                     >  7:300706/1‑149 (MQ=255)
                                                                                                                                cgcTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtggtg                  <  8:291727/149‑1 (MQ=255)
                                                                                                                                  cTCGAGCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGTGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTggtggt                <  5:211863/149‑1 (MQ=255)
                                                                                                                                   tCGACCCGCCGTTGAGAGCGGCGGTGGTGTCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTTGGGCTGGCTCGATCGGCTGGTGGTc               <  1:313520/149‑1 (MQ=255)
                                                                                                                                   tCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGAt                           >  3:123226/1‑137 (MQ=255)
                                                                                                                                   tCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGAt                           <  4:123226/137‑1 (MQ=255)
                                                                                                                                   tCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTggtggt                >  6:160787/1‑148 (MQ=255)
                                                                                                                                    cGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCa              <  3:198768/149‑1 (MQ=255)
                                                                                                                                     gACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTc               <  5:70889/147‑1 (MQ=255)
                                                                                                                                      caccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCACCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAAc                                               <  1:261498/113‑1 (MQ=255)
                                                                                                                                      caccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAAc                                               >  2:261498/3‑114 (MQ=255)
                                                                                                                                        ccgccgTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACtt          >  5:92352/1‑149 (MQ=255)
                                                                                                                                        ccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCACCGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACtt          >  6:19743/1‑149 (MQ=255)
                                                                                                                                          gccgTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTcc        <  2:58981/149‑1 (MQ=255)
                                                                                                                                            cgTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc      <  4:172688/149‑1 (MQ=255)
                                                                                                                                             gTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc      >  5:257830/1‑148 (MQ=255)
                                                                                                                                             gTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCg     <  8:183250/149‑1 (MQ=255)
                                                                                                                                              ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGa    <  6:92352/149‑1 (MQ=255)
                                                                                                                                                 aCAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGAcc  <  6:257830/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACC‑TCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACC  >  NC_002947/1034818‑1035110

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