Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,034,960 T→C 6.3% L333L (TTG→CTG)  PP_0892 → phospholipase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,034,9600TC 6.3% 270.4 / 5.4 94L333L (TTG→CTG) PP_0892phospholipase family protein
Reads supporting (aligned to +/- strand):  ref base T (42/46);  new base C (3/3);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

GGCGTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGA  >  NC_002947/1034813‑1035108
                                                                                                                                                   |                                                                                                                                                       
gtcGTACGGATTTTCGAGTACCAGCCGTGCTTCTTGCAACTTAATATCGCTTTGGTAGACGACTGGGTCATCGTCGGTTCGTGCAACTTCGACCACTTGAACTTGCGCTTCAACCTTTAAGCCAACATCGAAGCGCTCGACCCGCAGtt                                                                                                                                                        <  2:357752/147‑1 (MQ=255)
   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                                                                                                                         >  2:278523/1‑148 (MQ=255)
                                ccTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAgcgc                                                                                                                         <  4:167065/148‑1 (MQ=255)
                                  tGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGtt                                                                                                                                                        <  8:159588/115‑1 (MQ=255)
                                  tGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGtt                                                                                                                                                        >  7:159588/1‑115 (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                                                                                                       <  1:288535/138‑1 (MQ=255)
                                                            gacTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCgaa                                                                                                       >  2:288535/1‑138 (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                                                       <  6:38585/149‑1 (MQ=255)
                                                                                              aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                       <  7:4970/149‑1 (MQ=255)
                                                                                              aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                       <  4:330490/149‑1 (MQ=255)
                                                                                              aCTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAgg                                                       <  4:301365/149‑1 (MQ=255)
                                                                                               cTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAggg                                                      >  4:103906/1‑149 (MQ=255)
                                                                                               cTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAggg                                                      >  3:172688/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                                        >  4:225731/1‑146 (MQ=255)
                                                                                                                aCCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATc                                        <  3:225731/146‑1 (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)
                                                                                                                                                   |                                                                                                                                                       
GGCGTACGGATTTTCGAGTACCAGCCGTGCTTCCTGCACCTGAAGATGGCTGTGGTAGACGACTGGGTCAGCGTCGGTTCGTGCAACTTCGACCACTGGAACTTGCGCTTCAACCTTGAAGCCAACATCGAAGCGCTCGACCCGCCGTTGACAGCGGCGGTGGTGGCCAGCTTCGAGCGCGACTTCGCCCAGAGCGAAG‑‑AAGTCGACCTCGCCCACTGGCATGCGCGCCCTTGGTGGCGCAGGGTGAAGCAACGGATCTGGGGCTGGCTCGATCGGCTGGTGGTCAACTTCCTCGA  >  NC_002947/1034813‑1035108

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: