Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,034,963 A→T 5.3% T334S (ACA→TCA)  PP_0892 → phospholipase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,034,9630AT 5.3% 292.5 / 4.1 94T334S (ACA→TCA) PP_0892phospholipase family protein
Reads supporting (aligned to +/- strand):  ref base A (41/48);  new base T (3/2);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

GTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACC  >  NC_002947/1034816‑1035110
                                                                                                                                                   |                                                                                                                                                      
gTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGa                                                                                                                                                        >  6:273843/1‑147 (MQ=255)
  aCGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAg                                                                                                                                                     >  3:257053/1‑149 (MQ=255)
     gATTTTCAAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAgcgg                                                                                                                                                  <  7:288270/149‑1 (MQ=255)
            gAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTTGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTAGACCCGCCGTTGACAgcgg                                                                                                                                                  <  4:351475/142‑1 (MQ=255)
            gAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTTGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTAGACCCGCCGTTGACAgcgg                                                                                                                                                  >  3:351475/1‑142 (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                                                                                                                           <  4:167065/148‑1 (MQ=255)
                             ccTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAgcgc                                                                                                                           >  2:278523/1‑148 (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                                                                                              >  2:2168/1‑149 (MQ=255)
                                                       acgacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGAc                                                                                              <  4:257053/149‑1 (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‑‑AAGTCGACCt                                                                                           >  8:79716/1‑149 (MQ=255)
                                                              ggTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGCTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGc                                                                                        <  7:201488/147‑1 (MQ=255)
                                                                tCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCAc                                                                                    >  7:106293/1‑149 (MQ=255)
                                                                  aGCGGCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTg                                                                                  <  3:213896/149‑1 (MQ=255)
                                                                         gTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTg                                                                                  <  6:146700/142‑1 (MQ=255)
                                                                         gTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTg                                                                                  >  5:146700/1‑142 (MQ=255)
                                                                         gTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATgcg                                                                           >  8:4970/1‑149 (MQ=255)
                                                                             gTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGcc                                                                       >  8:336569/1‑149 (MQ=255)
                                                                             gTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGcc                                                                       >  1:58981/1‑149 (MQ=255)
                                                                             gTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGccc                                                                      <  5:273843/149‑1 (MQ=255)
                                                                               gCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCaga                                                                                                              >  6:237170/1‑111 (MQ=255)
                                                                               gCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCaga                                                                                                              <  5:237170/111‑1 (MQ=255)
                                                                                        aCCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGgcg                                                             >  5:49952/1‑148 (MQ=255)
                                                                                          cACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAg                                                          <  6:369056/149‑1 (MQ=255)
                                                                                           aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  8:112534/149‑1 (MQ=255)
                                                                                           aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  4:330490/149‑1 (MQ=255)
                                                                                           aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  4:301365/149‑1 (MQ=255)
                                                                                           aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  6:38585/149‑1 (MQ=255)
                                                                                           aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                         <  7:4970/149‑1 (MQ=255)
                                                                                            cTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAggg                                                        >  3:172688/1‑149 (MQ=255)
                                                                                            cTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAggg                                                        >  4:103906/1‑149 (MQ=255)
                                                                                               gAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTg                                                      >  6:282167/1‑148 (MQ=255)
                                                                                               gAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTCGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTg                                                      >  7:30644/1‑148 (MQ=255)
                                                                                                aaCTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGaa                                                    <  7:79716/149‑1 (MQ=255)
                                                                                                aaCTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGaa                                                    <  8:33314/149‑1 (MQ=255)
                                                                                                 aCTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAg                                                   <  1:278523/149‑1 (MQ=255)
                                                                                                      cgcTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAAc                                               <  6:357105/148‑1 (MQ=255)
                                                                                                         ttCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCCCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCAAGCGCGCCCTTGGTGGAGCAGCGTGAAGAAACGGa                                            >  2:290421/1‑148 (MQ=255)
                                                                                                           cAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATc                                          <  5:282167/148‑1 (MQ=255)
                                                                                                           cAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCt                                         >  2:169676/1‑149 (MQ=255)
                                                                                                            aaCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTg                                        <  5:180816/149‑1 (MQ=255)
                                                                                                            aaCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTg                                        >  6:180816/1‑149 (MQ=255)
                                                                                                            aaCCTTGAAGCCAACATCGAAGCCCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTg                                        <  3:164470/149‑1 (MQ=255)
                                                                                                             aCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATc                                          <  3:225731/146‑1 (MQ=255)
                                                                                                             aCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATc                                          >  4:225731/1‑146 (MQ=255)
                                                                                                              ccTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAGCCTAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTg                                        <  1:290421/149‑1 (MQ=255)
                                                                                                               cTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTgggg                                     >  3:50915/1‑149 (MQ=255)
                                                                                                                ttGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGggc                                    >  4:293468/1‑149 (MQ=255)
                                                                                                                   aaGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGggctgg                                 <  8:30644/149‑1 (MQ=255)
                                                                                                                       cAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTtggtgg                                                               >  5:202214/1‑115 (MQ=255)
                                                                                                                       cAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTtggtgg                                                               <  6:202214/115‑1 (MQ=255)
                                                                                                                       cAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTgg                                       <  6:23407/139‑1 (MQ=255)
                                                                                                                       cAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTgg                                       >  5:23407/1‑139 (MQ=255)
                                                                                                                        aaCATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGACCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGTGCTGGCTCGa                            >  6:211863/1‑149 (MQ=255)
                                                                                                                         aCATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGAt                           <  4:50915/149‑1 (MQ=255)
                                                                                                                          cATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATc                          <  7:336569/149‑1 (MQ=255)
                                                                                                                             cGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGc                       >  7:183250/1‑149 (MQ=255)
                                                                                                                               aaGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGTTCGATCGGCtg                     <  3:62402/149‑1 (MQ=255)
                                                                                                                               aaGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtg                     >  7:300706/1‑149 (MQ=255)
                                                                                                                                  cgcTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCtggtg                  <  8:291727/149‑1 (MQ=255)
                                                                                                                                    cTCGAGCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGTGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTggtggt                <  5:211863/149‑1 (MQ=255)
                                                                                                                                     tCGACCCGCCGTTGAGAGCGGCGGTGGTGTCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTTGGGCTGGCTCGATCGGCTGGTGGTc               <  1:313520/149‑1 (MQ=255)
                                                                                                                                     tCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGAt                           >  3:123226/1‑137 (MQ=255)
                                                                                                                                     tCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGAt                           <  4:123226/137‑1 (MQ=255)
                                                                                                                                     tCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTggtggt                >  6:160787/1‑148 (MQ=255)
                                                                                                                                      cGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCa              <  3:198768/149‑1 (MQ=255)
                                                                                                                                       gACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTc               <  5:70889/147‑1 (MQ=255)
                                                                                                                                        caccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCACCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAAc                                               <  1:261498/113‑1 (MQ=255)
                                                                                                                                        caccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAAc                                               >  2:261498/3‑114 (MQ=255)
                                                                                                                                          ccgccgTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACtt          >  5:92352/1‑149 (MQ=255)
                                                                                                                                          ccgccgCTGTCAACGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCACCGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACtt          >  6:19743/1‑149 (MQ=255)
                                                                                                                                            gccgTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTcc        <  2:58981/149‑1 (MQ=255)
                                                                                                                                              cgTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc      <  4:172688/149‑1 (MQ=255)
                                                                                                                                               gTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTc      >  5:257830/1‑148 (MQ=255)
                                                                                                                                               gTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCg     <  8:183250/149‑1 (MQ=255)
                                                                                                                                                ttGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGa    <  6:92352/149‑1 (MQ=255)
                                                                                                                                                   aCAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGAcc  <  6:257830/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGACC  >  NC_002947/1034816‑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: