Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,066,370 T→A 5.4% L154Q (CTG→CAG)  PP_2678 → hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,066,3700TA 5.4% 294.9 / 5.0 93L154Q (CTG→CAG) PP_2678hydrolase
Reads supporting (aligned to +/- strand):  ref base T (45/43);  new base A (2/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGACCACGC  >  NC_002947/3066224‑3066519
                                                                                                                                                  |                                                                                                                                                     
ccAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCa                                                                                                                                                    >  2:716962/1‑150 (MQ=255)
 cAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAg                                                                                                                                                   >  1:2353401/1‑150 (MQ=255)
  aGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGc                                                                                                                                                  >  1:756745/1‑150 (MQ=255)
  aGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGc                                                                                                                                                  >  2:594370/1‑150 (MQ=255)
  aGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGc                                                                                                                                                  >  2:2073264/1‑150 (MQ=255)
   ggCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGcc                                                                                                                                                 <  2:1319052/150‑1 (MQ=255)
   ggCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGcc                                                                                                                                                 >  1:607784/1‑150 (MQ=255)
       tttGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGc                                                                                                                                             <  2:787066/150‑1 (MQ=255)
       tttGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGc                                                                                                                                             <  1:925382/150‑1 (MQ=255)
           ccGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGc                                                                                                                                         <  1:232114/150‑1 (MQ=255)
             gATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCACGCCCGGCTGCAGCACCCgggtcggcag                                                                                                                                       <  2:454218/150‑8 (MQ=255)
                gtgCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAg                                                                                                                                    <  1:780560/150‑1 (MQ=255)
                 tgCCGATCGGTGCCCTGGCCGGCACCCGCGACATGCTGCGCGAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAgg                                                                                                                                   <  2:529019/150‑1 (MQ=255)
                     gATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtggg                                                                                                                               >  1:342908/1‑150 (MQ=255)
                       tCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                             <  2:607784/150‑1 (MQ=255)
                        cGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGGTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGgtgggtgg                                                                                                                            <  1:1983073/150‑1 (MQ=255)
                        cGGTGCCCTGGCCGGCACCCGCCACCGGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGgtgggtgg                                                                                                                            <  1:733086/150‑1 (MQ=255)
                          gTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGCTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgg                                                                                                                                               <  2:2203150/129‑1 (MQ=255)
                          gTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgg                                                                                                                                               >  1:2203150/1‑129 (MQ=255)
                              ccTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCg                                                                                                                      >  2:2309873/1‑150 (MQ=255)
                                tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                              <  1:1678501/140‑1 (MQ=255)
                                tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                              >  2:1678501/1‑140 (MQ=255)
                                tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAg                                                                                                                                    >  2:139797/1‑134 (MQ=255)
                                tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAg                                                                                                                                    <  1:139797/134‑1 (MQ=255)
                                  gCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgggtc                                                                                                                                            >  2:593231/1‑122 (MQ=255)
                                  gCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgggtc                                                                                                                                            <  1:593231/124‑3 (MQ=255)
                                    aagCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAg                                                                                                                >  2:1902570/3‑150 (MQ=255)
                                       cACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATTCTCGGCACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGGGGTGGGTGGGCACCTGCTGAGGtt                                                                                                             >  1:1489512/1‑150 (MQ=255)
                                        aCCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTg                                                                                                            >  1:2389573/1‑150 (MQ=255)
                                          ccGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGACGGGGACTGGATGCTCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgct                                                                                                          >  1:1499851/1‑150 (MQ=255)
                                           cgcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtg                                                                                                                                 >  1:2075293/1‑126 (MQ=255)
                                           cgcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtg                                                                                                                                 <  2:2075293/126‑1 (MQ=255)
                                            gcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgg                                                                                                        <  2:1155675/150‑1 (MQ=255)
                                            gcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgg                                                                                                        <  2:952946/150‑1 (MQ=255)
                                              gACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgggg                                                                                                      >  2:1115427/1‑150 (MQ=255)
                                                 ctgctgCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTCGGTGGGCACCGCCTGAGGTTgctggggctg                                                                                                   >  2:663418/1‑150 (MQ=255)
                                                  tgctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGGCTGGTTGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGGGGGCACCGCCTGAGGTTgctggggctgg                                                                                                  >  1:1623594/1‑150 (MQ=255)
                                                  tgctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctggggctgg                                                                                                  >  1:755742/1‑150 (MQ=255)
                                                    ctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggc                                                                                                <  1:2309873/150‑1 (MQ=255)
                                                     tgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggcg                                                                                               <  2:1884751/150‑1 (MQ=255)
                                                      gcgcCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggcgg                                                                                              <  1:716962/150‑1 (MQ=255)
                                                         ccAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCa                                                                                           >  2:789985/1‑150 (MQ=255)
                                                            gCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                                       >  1:2318363/1‑103 (MQ=255)
                                                            gCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                                       <  2:2318363/103‑1 (MQ=255)
                                                            gCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATAc                                                                                        >  2:2185495/1‑150 (MQ=255)
                                                                ggCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGt                                                                                    <  2:2389573/150‑1 (MQ=255)
                                                                   gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                             >  2:1760891/1‑106 (MQ=255)
                                                                   gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                             <  1:1760891/106‑1 (MQ=255)
                                                                   gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGcc                                                                                 >  2:2228099/1‑150 (MQ=255)
                                                                   gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGAAGCTGCCGACCCGGGAGCTGCAGCCGGGCGTGAAGGAGGTGGGTGGGGACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATCGGGGTAcc                                                                                 >  2:1166293/1‑150 (MQ=255)
                                                                     tGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGa                                                                               <  2:1443213/150‑1 (MQ=255)
                                                                       ccATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCa                                                                                                                                                    >  2:53758/1‑79 (MQ=255)
                                                                       ccATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCa                                                                                                                                                    <  1:53758/79‑1 (MQ=255)
                                                                         aTGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTg                                                                           >  1:315153/1‑150 (MQ=255)
                                                                          tGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTgg                                                                          >  1:339742/1‑150 (MQ=255)
                                                                           ggCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGc                                                                         <  2:315153/150‑1 (MQ=255)
                                                                              cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCa                                                                                                                                        <  1:1666049/84‑1 (MQ=255)
                                                                              cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCa                                                                                                                                        >  2:1666049/1‑84 (MQ=255)
                                                                              cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGAt                                                                      >  2:1412448/1‑150 (MQ=255)
                                                                                  aaCATGTATCGCCTGGTCGGCGACTGGATGTGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTc                                                                  <  2:2198294/150‑1 (MQ=255)
                                                                                     aTGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGAt                                                               >  2:847321/1‑150 (MQ=255)
                                                                                     aTGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGAt                                                               >  1:230515/1‑150 (MQ=255)
                                                                                     aTGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGAt                                                               <  1:1412448/150‑1 (MQ=255)
                                                                                       gTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGa                                                             >  1:1751072/1‑150 (MQ=255)
                                                                                         aTCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGaaa                                                           <  2:1903475/150‑1 (MQ=255)
                                                                                          tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                              >  2:1436933/1‑82 (MQ=255)
                                                                                          tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                              <  1:1436933/82‑1 (MQ=255)
                                                                                          tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTg                                                                                                                  <  2:1014651/94‑1 (MQ=255)
                                                                                          tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTg                                                                                                                  >  1:1014651/1‑94 (MQ=255)
                                                                                            gCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGAGCTGGCGATAATCGAGTAAAAGa                                                        >  2:1268737/1‑150 (MQ=255)
                                                                                               tGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTTCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTCAGACGTGTCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCg                                                     <  1:1268737/150‑1 (MQ=255)
                                                                                                      cGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGAGCGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCt                                              <  2:342908/150‑1 (MQ=255)
                                                                                                       gACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGa                                                                                                                 >  2:1726434/1‑82 (MQ=255)
                                                                                                       gACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGa                                                                                                                 <  1:1726434/82‑1 (MQ=255)
                                                                                                        aCTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGa                                                                <  1:650200/130‑1 (MQ=255)
                                                                                                        aCTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGa                                                                >  2:650200/1‑130 (MQ=255)
                                                                                                          tGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTc                                          <  2:1751072/150‑1 (MQ=255)
                                                                                                             aTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGcc                                       <  1:594370/150‑1 (MQ=255)
                                                                                                              tGCGCGGTACCGAGGTGGTGCTGGCGACCCGGGTGCTGCAGGCGGGCGTGCATGAGGTGGGTGGGCACCGCCGGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCg                                      >  2:1853667/1‑150 (MQ=255)
                                                                                                               gcgcGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                                       >  2:1732908/1‑52 (MQ=255)
                                                                                                               gcgcGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                                       <  1:1732908/52‑1 (MQ=255)
                                                                                                                      aCCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTg                              >  2:586242/1‑150 (MQ=255)
                                                                                                                              ggtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGACGGGGCTGGGCGGGCTGACGGGGTCCGTCCTGGCGATTTTCTATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTa                      <  2:2078868/150‑1 (MQ=255)
                                                                                                                              ggtcCTGCCGACCCGGGTGCTGCAGGCGGGCGTGCATGTGTGGGGTGGGCACCGCCTGAGGTTGCTGGGGCTAGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTa                      <  1:1737510/150‑1 (MQ=255)
                                                                                                                               gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCTCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc                     <  1:1341134/150‑1 (MQ=255)
                                                                                                                               gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc                     <  1:1115427/150‑1 (MQ=255)
                                                                                                                               gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc                     <  1:586242/150‑1 (MQ=255)
                                                                                                                                  ctgcCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGCCGGTTTTCGATGAAAAGACCGGGGTGCTTTGCGCCGGCGACCTCGTGCTCTACCaa                  >  2:542769/1‑150 (MQ=255)
                                                                                                                                   tgcCGACCCGGGTGCTGCAGGCGGGCGTGCATGACGTGGGTGCGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGAGCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAAc                 <  2:1082686/150‑1 (MQ=255)
                                                                                                                                   tgcCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAAc                 >  2:2030359/1‑150 (MQ=255)
                                                                                                                                          ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGCCTGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt          >  1:1039658/1‑150 (MQ=255)
                                                                                                                                          ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt          >  2:2274310/1‑150 (MQ=255)
                                                                                                                                          ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt          <  1:2030359/150‑1 (MQ=255)
                                                                                                                                                gctgcAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGACCAc    <  1:847321/150‑1 (MQ=255)
                                                                                                                                                 ctgcAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGc                                               >  1:1369348/1‑106 (MQ=255)
                                                                                                                                                 ctgcAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGc                                               <  2:1369348/106‑1 (MQ=255)
                                                                                                                                                  tgcAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGACCACGc  >  2:335000/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGACCACGC  >  NC_002947/3066224‑3066519

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: