Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,554,556 A→C 7.7% S138R (AGT→AGG)  PP_4041 ← glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,554,5560AC 7.7% 192.3 / 4.7 72S138R (AGT→AGG) PP_4041glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base A (29/37);  new base C (2/4);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGGCATGAAGT  >  NC_002947/4554409‑4554704
                                                                                                                                                   |                                                                                                                                                    
cTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGGGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGAGGGCACt                                                                                                                                                    <  2:422683/150‑1 (MQ=255)
cTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGTGCCGTCGCCGGCCCt                                                                                                                                                    <  1:175747/150‑1 (MQ=255)
cTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACt                                                                                                                                                    >  2:1797387/1‑150 (MQ=255)
cTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTCTGGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACt                                                                                                                                                    <  1:497495/150‑1 (MQ=255)
  gctccGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGTGCCGTCGCCGGCCCt                                                                                                                                                    <  1:1019643/145‑1 (MQ=255)
  gCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGTGCCGTCGCCGGCCCt                                                                                                                                                    >  2:1019643/1‑148 (MQ=255)
  gCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGCGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGa                                                                                                                                                  <  2:155947/150‑1 (MQ=255)
   cACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAg                                                                                                                                                 <  1:741634/150‑1 (MQ=255)
       ggCAAGTGCGAGGGGTGCTAACACCGGCAGTGGAGTTGGGGCGTGCCCGTGCTGACGCGAAAGCTGCACTACATAGGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGTGCCGTCGCCGGCCCTGaccggt                                                                                                                                             <  1:329431/150‑6 (MQ=255)
        gCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCa                                                                                                                                            >  1:271388/1‑150 (MQ=255)
           aGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTg                                                                                                                                         <  1:1046015/150‑1 (MQ=255)
           aGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTg                                                                                                                                         <  1:1797387/150‑1 (MQ=255)
               cGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTc                                                                                                                                     <  2:271388/150‑1 (MQ=255)
                 aGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGa                                                                                                                                   <  1:862543/150‑1 (MQ=255)
                 aGGGCTGGAGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGa                                                                                                                                   >  2:1613086/1‑150 (MQ=255)
                  gggCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAg                                                                                                                                  >  2:994697/1‑150 (MQ=255)
                             gCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCAccc                                                                                                                       >  2:1006628/1‑150 (MQ=255)
                              ccggttGAAGACCTGGCGCTGGCCCGTGGTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCa                                                                                                                      <  2:180892/145‑1 (MQ=255)
                              cAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTCAGCTGGCACGGTGCCGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCa                                                                                                                      <  1:797720/150‑1 (MQ=255)
                                ggCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGTTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAgg                                                                                                                    <  1:1018999/150‑1 (MQ=255)
                                       gACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGt                                                                                                                              <  2:722739/133‑1 (MQ=255)
                                       gACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGt                                                                                                                              >  1:722739/1‑133 (MQ=255)
                                               gcTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGt                                                                                                     >  2:1218443/1‑150 (MQ=255)
                                                cTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTa                                                                                                    <  2:1624036/150‑1 (MQ=255)
                                                  ggCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACGGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCCGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGt                                                                                                  >  2:506554/1‑150 (MQ=255)
                                                   gCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGcc                                                                                                       <  2:733601/144‑1 (MQ=255)
                                                   gCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGcc                                                                                                       >  1:733601/1‑144 (MQ=255)
                                                     ccGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGa                                                                                               >  1:136149/1‑150 (MQ=255)
                                                              cgcgAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGa                                                                                      <  1:749233/150‑1 (MQ=255)
                                                               gcgAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCTGGGGCACTGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGAc                                                                                     >  2:1133380/1‑150 (MQ=255)
                                                                cgAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGTCGTAGTCGAAGCGCATGAcc                                                                                    <  1:1613086/150‑1 (MQ=255)
                                                                  aaaGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAg                                                                                  <  1:994697/150‑1 (MQ=255)
                                                                  aaaGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAg                                                                                  >  1:265249/1‑150 (MQ=255)
                                                                    aGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGt                                                                                <  2:739754/150‑1 (MQ=255)
                                                                       ggCAGCCGGTATGGTGGTCGCGGGCTGGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGCGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTccc                                                                             >  2:817584/1‑149 (MQ=255)
                                                                               gTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGaaa                                                                     <  2:1412387/150‑1 (MQ=255)
                                                                               gTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGaaa                                                                     <  2:6362/150‑1 (MQ=255)
                                                                               gTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGAACTGGGTCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGaaa                                                                     <  1:1453500/150‑1 (MQ=255)
                                                                                  tggtggTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATgg                                                                  <  1:361020/150‑1 (MQ=255)
                                                                                     tggtCGCGGGCTTGCGGGTGCACGGTGCGGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGtt                                                               >  1:449604/1‑150 (MQ=255)
                                                                                       gtCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTc                                                             <  2:136149/150‑1 (MQ=255)
                                                                                        tCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCg                                                            <  1:865408/150‑1 (MQ=255)
                                                                                          gcgGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCCACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGcc                                                          <  1:1870340/150‑1 (MQ=255)
                                                                                           ggggCTTGCGGGGGGACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGATGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCa                                                         <  2:140060/149‑1 (MQ=255)
                                                                                           cgGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCa                                                         <  2:265249/150‑1 (MQ=255)
                                                                                             ggCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGa                                                       <  1:252030/150‑1 (MQ=255)
                                                                                               cTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACa                                                     >  1:682807/1‑150 (MQ=255)
                                                                                                  gCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGc                                                  <  1:1218443/150‑1 (MQ=255)
                                                                                                    gggggCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCt                                                >  1:278898/1‑150 (MQ=255)
                                                                                                    gggggCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCt                                                >  1:273578/1‑150 (MQ=255)
                                                                                                            ggTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGAt                                        >  1:1299888/1‑150 (MQ=255)
                                                                                                             gTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATg                                       <  2:273578/150‑1 (MQ=255)
                                                                                                              tgctgcACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATgc                                      <  1:652383/150‑1 (MQ=255)
                                                                                                                   cacaGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGGCAGACAAGCCTTCGACGATGCGcacc                                 >  2:747728/1‑150 (MQ=255)
                                                                                                                   cacaGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGcacc                                 >  2:749062/1‑150 (MQ=255)
                                                                                                                   cacaGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGcacc                                 >  2:1477631/1‑150 (MQ=255)
                                                                                                                    acaGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGcacca                                >  1:1088069/1‑150 (MQ=255)
                                                                                                                     caGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGcaccac                               >  1:444964/1‑150 (MQ=255)
                                                                                                                      aGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAACCGGTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGGCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCaccacc                              >  1:347345/1‑150 (MQ=255)
                                                                                                                          tCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAAt                          >  1:869783/1‑150 (MQ=255)
                                                                                                                              ccGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGAc                      >  1:358407/1‑150 (MQ=255)
                                                                                                                                 gTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTc                   <  1:1133380/150‑1 (MQ=255)
                                                                                                                                 ctCAGTGCCGTCGCCGGCCCTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCa                                                         <  1:1004502/111‑1 (MQ=255)
                                                                                                                                 ctCAGTGCCGTCGCCGGCCCTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCa                                                         >  2:1004502/2‑112 (MQ=255)
                                                                                                                                   cAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGc                 <  1:1189353/150‑1 (MQ=255)
                                                                                                                                   cAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGc                 <  2:1577192/150‑1 (MQ=255)
                                                                                                                                   cAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGc                 <  1:1189365/150‑1 (MQ=255)
                                                                                                                                      ggCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGGAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCa              >  2:167865/1‑150 (MQ=255)
                                                                                                                                          ggCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGAAGATGCGCACCACCGAATTGAATTCGCCCAgggg          >  1:650004/1‑150 (MQ=255)
                                                                                                                                                ggCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGGCATGaa    <  1:1477631/150‑1 (MQ=255)
                                                                                                                                                 gCACAGACGGTCAGGGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGGCATGAAg   <  1:506554/150‑1 (MQ=255)
                                                                                                                                                  cACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGGCATGAAGt  >  1:53007/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGGCATGAAGT  >  NC_002947/4554409‑4554704

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