Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,334,996:1 +T 6.0% coding (609/1704 nt) PP_2932 ← amidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,334,9961.T 6.0% 330.4 / 14.1 84coding (609/1704 nt)PP_2932amidase family protein
Reads supporting (aligned to +/- strand):  ref base . (37/42);  new base T (3/2);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

TGCAGGCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGT  >  NC_002947/3334847‑3335142
                                                                                                                                                      |                                                                                                                                                  
tGCAGGCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCcg                                                                                                                                                     >  1:1929986/1‑150 (MQ=255)
tGCAGGCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGACGGCCATGGTGCGGGCATAGGGCACCACCACGACCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCcg                                                                                                                                                     >  1:613321/1‑150 (MQ=255)
   aGGCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCcg‑cgc                                                                                                                                                 >  2:1825316/1‑150 (MQ=255)
     gcgcCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGGCAGTTGCCACGCACCGAAATCACCCcg‑cgcgg                                                                                                                                               >  2:495849/1‑149 (MQ=255)
       gcCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACg                                                                                                                                             <  2:1870005/150‑1 (MQ=255)
        ccAAAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGGTTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgg                                                                                                                                            <  1:379733/150‑1 (MQ=255)
         cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CTCGACggg                                                                                                                                           <  1:1523935/150‑1 (MQ=255)
         cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGCCggg                                                                                                                                           >  2:804331/1‑150 (MQ=255)
         cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                           <  1:2095008/150‑1 (MQ=255)
         cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                           <  1:1523852/150‑1 (MQ=255)
         cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTAACACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                           <  1:1946536/150‑1 (MQ=255)
             ggTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGAACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGgtgt                                                                                                                                       >  1:1071502/1‑150 (MQ=255)
                cgccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAgg                                                                                                                                    <  2:14779/150‑1 (MQ=255)
                   cgcgCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgcgc                                                                                                                                 >  1:750538/1‑150 (MQ=255)
                    gcgcTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGGGTAGGCGca                                                                                                                                >  1:327041/1‑150 (MQ=255)
                       cTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGAGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAg                                                                                                                             <  2:1998215/150‑1 (MQ=255)
                                gTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGATATCACCCCG‑CGCGTCGGGGTGTAGGCGCACAGGCCgttgtt                                                                                                                    <  2:1465791/150‑1 (MQ=255)
                                 tCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCTGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGttgttg                                                                                                                   <  2:396628/150‑1 (MQ=255)
                                 tCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGttgttg                                                                                                                   <  1:600846/150‑1 (MQ=255)
                                  cATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAgg                                                                                                                            <  2:2323358/140‑1 (MQ=255)
                                  cATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAgg                                                                                                                            >  1:2323358/1‑140 (MQ=255)
                                   aTCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGa                                                                                                                 <  1:1619811/150‑1 (MQ=255)
                                    tCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGcc                                                                                                                          >  1:1941887/1‑140 (MQ=255)
                                    tCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGcc                                                                                                                          <  2:1941887/140‑1 (MQ=255)
                                    tCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGAc                                                                                                                >  1:2131898/1‑150 (MQ=255)
                                    tCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGAc                                                                                                                >  2:2147222/1‑150 (MQ=255)
                                       gCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgggg                                                                                                                                          <  2:470154/121‑1 (MQ=255)
                                       gCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgggg                                                                                                                                          >  1:470154/1‑121 (MQ=255)
                                          accaccACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACAACCACGTCCATGGTCGACAAGGACGGCCAGTTGCCAAGCAAAGAAATAATTGCG‑TGCGAGGGGGTGTAGGCGCAGAGGCCGTTGTTGGACGCCggg                                                                                                          >  1:1117171/1‑150 (MQ=255)
                                           ccaccaCGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGc                                                                                                         <  1:433928/150‑1 (MQ=255)
                                            caccacGTCGAGGATTTCAAGGAGGTCGGCTATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCCGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGcc                                                                                                        >  2:644723/1‑150 (MQ=255)
                                            caccacGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGcc                                                                                                        >  2:1420951/1‑150 (MQ=255)
                                                    cGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGc                                                                                                <  2:766268/150‑1 (MQ=255)
                                                     gAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACg                                                                                                     >  2:774518/1‑144 (MQ=255)
                                                     gAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACg                                                                                                     <  1:774518/144‑1 (MQ=255)
                                                       ggATTGCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTgg                                                                                             >  1:1365258/1‑150 (MQ=255)
                                                        gATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAg                                                                                                                                     >  1:1523491/1‑109 (MQ=255)
                                                        gATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAg                                                                                                                                     <  2:1523491/109‑1 (MQ=255)
                                                        gATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGAGAAACGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGa                                                                                            <  1:883508/150‑1 (MQ=255)
                                                          tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgatttc                                                                                                                                   >  2:679675/1‑104 (MQ=255)
                                                          tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgatttc                                                                                                                                   <  1:679675/109‑6 (MQ=255)
                                                          tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAcc                                                                                          <  1:947405/150‑1 (MQ=255)
                                                            tCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAg                                                                                        <  2:1538518/150‑1 (MQ=255)
                                                             cAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTAGGCGTTAGCGGCCCGTTGCCACGCAGCGCAATCACACCG‑CGCGACGGGGTGGGGGCTCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAgg                                                                                       >  2:2101074/1‑150 (MQ=255)
                                                              aaGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCA‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTg                                                                                              >  2:108083/1‑142 (MQ=255)
                                                              aaGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCA‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTg                                                                                              <  1:108083/142‑1 (MQ=255)
                                                                     tCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCg                                                                               <  1:476552/150‑1 (MQ=255)
                                                                      cGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCgg                                                                              <  2:1433169/150‑1 (MQ=255)
                                                                       ggCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                      >  1:1153840/1‑141 (MQ=255)
                                                                       ggCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                      <  2:1153840/141‑1 (MQ=255)
                                                                          cATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggccag                                                                          <  1:1420951/150‑1 (MQ=255)
                                                                           aTGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggccagg                                                                         <  1:419442/150‑1 (MQ=255)
                                                                                    gCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAg                                                                                                                             >  2:513123/1‑89 (MQ=255)
                                                                                    gCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAg                                                                                                                             <  1:513123/89‑1 (MQ=255)
                                                                                      aTAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCt                                                              >  1:2186430/1‑150 (MQ=255)
                                                                                         gggCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGaa                                                           <  1:2147222/150‑1 (MQ=255)
                                                                                           gaaccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGc                                                         <  2:1071502/148‑1 (MQ=255)
                                                                                           gCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGc                                                         >  2:545688/1‑150 (MQ=255)
                                                                                              ccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTgg                                                      >  1:2390159/1‑150 (MQ=255)
                                                                                               caccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAACAGGTCTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggc                                                     >  2:302533/1‑150 (MQ=255)
                                                                                                 ccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggcgg                                                   <  1:1825316/150‑1 (MQ=255)
                                                                                                        cATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggc                                            >  1:1748400/1‑150 (MQ=255)
                                                                                                         aTGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggcg                                           <  2:613321/150‑1 (MQ=255)
                                                                                                           ggTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggcggt                                         <  1:230605/150‑1 (MQ=255)
                                                                                                             tCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTAc                                       <  1:2361676/150‑1 (MQ=255)
                                                                                                               ggCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTCTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCg                                     >  2:171646/1‑150 (MQ=255)
                                                                                                                        ggCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTg                            <  2:1748400/150‑1 (MQ=255)
                                                                                                                         gCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTgg                           >  2:922472/1‑150 (MQ=255)
                                                                                                                         gCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTgg                           >  1:1358625/1‑150 (MQ=255)
                                                                                                                         gCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTgg                           <  2:2186430/150‑1 (MQ=255)
                                                                                                                           cAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGAc                         <  1:20051/150‑1 (MQ=255)
                                                                                                                             gTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGa                       >  2:1743698/1‑150 (MQ=255)
                                                                                                                                  cACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCAGGAACCGGGTTCTTCGGCCAGGCCGGCGGCGAAAAAGCCGGCATTGCTTTCTGTACCTGCACCGTTGGACGAGCCCg                  >  1:532393/1‑150 (MQ=255)
                                                                                                                                   aCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGa                 >  2:486390/1‑150 (MQ=255)
                                                                                                                                   aCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGCCCGTTGTTGGACGCCGGGCCCCGGCCGCTGGACCCGGTCTCTTCGGCCCGGATCTCGGCCCTCGACCAGTCTGCGGGGGCCGTCCCGGCACCGTTGGACGAGCCCGa                 >  2:2367419/1‑150 (MQ=255)
                                                                                                                                         cGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGaa           >  1:915425/1‑150 (MQ=255)
                                                                                                                                         cGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGaa           >  2:423849/1‑150 (MQ=255)
                                                                                                                                          gAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGANGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAg          <  1:495849/150‑1 (MQ=255)
                                                                                                                                          gAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAg          >  1:1089192/1‑150 (MQ=255)
                                                                                                                                           aaaTCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAgg         <  2:1358625/150‑1 (MQ=255)
                                                                                                                                            ctacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAg                                                                                        <  2:204091/67‑1 (MQ=37)
                                                                                                                                            ctacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAg                                                                                        >  1:204091/4‑70 (MQ=37)
                                                                                                                                             tacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCTCGGCCACAGCCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgc       >  1:1123668/3‑150 (MQ=255)
                                                                                                                                               cACCCCG‑CGCGACGGGGTGGAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgcgc     <  2:1119646/150‑1 (MQ=255)
                                                                                                                                                  cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGa                                                            <  2:1229788/92‑1 (MQ=255)
                                                                                                                                                  cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGa                                                            >  1:1229788/1‑92 (MQ=255)
                                                                                                                                                  cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt  >  1:1168993/1‑150 (MQ=255)
                                                                                                                                                  cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt  >  2:1637018/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TGCAGGCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGT  >  NC_002947/3334847‑3335142

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