Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,460 A→T 5.7% intergenic (+114/+28) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,4600AT 5.7% 494.4 / 15.2 157intergenic (+114/+28)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (71/73);  new base T (5/4);  total (79/78)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 160 total aligned reads displayed.

CAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAACT  >  NC_002947/4340311‑4340609
                                                                                                                                                     |                                                                                                                                                     
cAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTaa                                                                                                                                                       <  2:18514/150‑1 (MQ=255)
 aGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTaaa                                                                                                                                                      <  2:975000/150‑1 (MQ=255)
     gggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcg                                                                                                                                                  >  2:4307917/1‑150 (MQ=255)
      ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                                 <  1:4488513/150‑1 (MQ=255)
        cAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTATGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAACCGCGa                                                                                                                                               >  1:2613146/1‑150 (MQ=255)
         aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaa                                                                                                                                              >  2:683226/1‑150 (MQ=255)
          aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATCGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                                                                                                             <  2:3021159/150‑1 (MQ=255)
          aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                                                                                                             <  1:1201692/150‑1 (MQ=255)
          aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                             >  2:2898277/1‑150 (MQ=255)
           cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAg                                                                                                                                            <  2:2451729/150‑1 (MQ=255)
           cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGATAg                                                                                                                                            <  2:1164820/150‑1 (MQ=255)
           cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                                                                                                            <  2:2962114/150‑1 (MQ=255)
           cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                                                                                                            <  2:3536663/150‑1 (MQ=255)
           cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                                                                                                            <  1:394283/150‑1 (MQ=255)
            gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAACgg                                                                                                                                           >  2:2376124/1‑150 (MQ=255)
            gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAACCGCGAAAgg                                                                                                                                           >  1:3264725/1‑150 (MQ=255)
            gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGAAACCCCCGACAgg                                                                                                                                           >  1:2322859/1‑150 (MQ=255)
              aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGc                                                                                                                                         <  1:4084701/150‑1 (MQ=255)
               cTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGcc                                                                                                                                        <  1:1940959/150‑1 (MQ=255)
                ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTAGCGGGTAAACCCGGGAAAGGGCCg                                                                                                                                       >  1:2013747/1‑150 (MQ=255)
                 tAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGTAAAACCCCGGAAGGGCCgg                                                                                                                                      >  1:4354614/1‑150 (MQ=255)
                  aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGGTACCGGGGAAAACCGCGACAGGGCCGGc                                                                                                                                     >  1:3111107/1‑150 (MQ=255)
                    aGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCaa                                                                                                                                   >  1:1062167/1‑150 (MQ=255)
                        cATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGTAAACCGCGCAAGGGCCCGcccaaac                                                                                                                               >  1:1499619/1‑144 (MQ=255)
                          tCTGCAGTGAGGCGATGAGCAATAGACCCGCGTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCTGGATTCAATGAATCGGCCAGGGGTGCGGCTGGGCCGGCCCTTTGGCGGGTAAACCCGCGAATGGGCCGGCAAATCCTg                                                                                                                             <  2:3889022/150‑1 (MQ=255)
                               aGTGAGGCGATGAGCAATAGACCCGCCTGGGCGTTTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                        <  2:3456735/150‑1 (MQ=255)
                               aGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                        <  2:1011657/150‑1 (MQ=255)
                               aGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                        <  1:3361750/150‑1 (MQ=255)
                               aGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTATATTTGCGTAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCCGGTCGGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                        <  1:4459642/150‑1 (MQ=255)
                                gTGAGGCGATTAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGACCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAgg                                                                                                                       <  2:1287634/150‑1 (MQ=255)
                                 tgcggCGATGAGCAATAGCCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCTGGCTTCAATGAATCGGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGGGAAAGGGCCGGCAAATCCTGCTCAggg                                                                                                                      <  1:2959261/147‑1 (MQ=255)
                                    ggCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTc                                                                                                                   >  1:2629778/1‑150 (MQ=255)
                                       gATGAGCAATCGACCCCCCTGGGCGGGACTTTTGCTTATATTTGCGGAAACAGCCCGCAGGCTCCAACGAATCTGCCCGCGGTGCGGCCGGCCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCaaa                                                                                                                <  1:2246294/150‑1 (MQ=255)
                                       gATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGGAACCCGCGAAAGGGCCCGCCAATCCAGCGGAGGGCGCaaa                                                                                                                >  1:3603175/1‑150 (MQ=255)
                                         tGAGCAATAGACCCGCCGGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAAGGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAcc                                                                                                              <  2:3578670/150‑1 (MQ=255)
                                            gCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt                                                                                                           >  2:2362182/1‑150 (MQ=255)
                                             cAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAACCCCGACAGGGCCGGCCAATCCTGCTCAGGGCTCAAAACCGTc                                                                                                          >  2:126057/1‑150 (MQ=255)
                                              aaaaGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAACCCGCGAAGGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                         >  1:1819178/4‑150 (MQ=255)
                                              aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTACACCCGCGAAAGGGCCGGCACATCCTGCTCAGGGCTCAAACCCGTcc                                                                                                         >  2:2509234/1‑150 (MQ=255)
                                                     ccGCCTGGGCGGTTATTTTGCTTCTATTTGCGGAAACAGCCCGCTGGCTTCAATGAATCTGCCAGCGGTGCGGCAGTTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGa                                                                                                  <  1:3164028/150‑1 (MQ=255)
                                                     ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCACATCCTGCTCAGGGCTCAAACCCTTCCACACTGa                                                                                                  >  2:2831050/1‑150 (MQ=255)
                                                     ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGa                                                                                                  >  2:1419038/1‑150 (MQ=255)
                                                     ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAAACCAGCTCAGGGCTCAAACCCGTCCACAATGa                                                                                                  >  1:3698885/1‑150 (MQ=255)
                                                     ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGGACACCCGCGAAAGGGCCGGCCAATACTGCTCCGGGCTCAAAACCGTCCCCCATGa                                                                                                  >  1:942761/1‑150 (MQ=255)
                                                       gcgtGGGCGGTTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGATTCAATTAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGTTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATc                                                                                                <  1:4378029/147‑1 (MQ=255)
                                                       gcgtGGGCGGTTATTTTGCTTCTTTTTGCGTAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTGCGCGGTTAACCCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATc                                                                                                <  1:2109073/147‑1 (MQ=255)
                                                             ggggtttttttGCTTCTTTTTGCGGAACCAGCCCGCTGGCTTCATTTATTGTGCCTGGGGTGGGGCTGGTCGGGCCCTTGGGGGGTTAACCCGGCTAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCa                                                                                          <  1:27116/143‑1 (MQ=255)
                                                             gcgtttCTTTTGCTTTTATTTGCGTAAACAGCCGGCAGGATTCATTGAATCTGCCAGGGGTGCGGCAGTTCCGGCCCTTTCGGGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCa                                                                                          <  1:3604042/145‑1 (MQ=255)
                                                             gCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCa                                                                                          <  2:3762488/150‑1 (MQ=255)
                                                               gggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAAGATCACCTCAGc                                                                                        >  2:2726580/1‑150 (MQ=255)
                                                               gggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGTAAACCCCCGAAAGGGCCCGCAAATACAGCTCAGGGCTCAAACCCATCAACAATGATAACCTCAGc                                                                                        >  2:620500/1‑150 (MQ=255)
                                                                   cTTTTGCTTCTATTTGCGTAAACAGCCCGCTGGCTTCATTGAATCTGCCTGGGGTGCGGCTGGTCCGGCCCTTGGGCGGGTAAACCCGCTAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTc                                                                                    <  1:2396021/150‑1 (MQ=255)
                                                                    ttttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCg                                                                                   <  2:461985/150‑1 (MQ=255)
                                                                          ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAccc                                                                             >  2:316093/1‑150 (MQ=255)
                                                                          ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAACCTTCCACAATGATCACCTCAGCCATCGCCAccc                                                                             >  1:4490841/1‑150 (MQ=255)
                                                                              aTTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTg                                                                         <  2:3308998/150‑1 (MQ=255)
                                                                               tttGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTgg                                                                        >  2:4074164/1‑150 (MQ=255)
                                                                                   cGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGAt                                                                    <  2:3412937/150‑1 (MQ=255)
                                                                                    ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATg                                                                   <  2:892729/150‑1 (MQ=255)
                                                                                    ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGAAAACCCCCAAAGGGCCGGCCAAAACAGCGCAGGGCGCAAAACCATCCACAAATATCACCACAGCCATCGCCCCCCCCAGGCGCTg                                                                   >  1:2164241/1‑150 (MQ=255)
                                                                                      aaaCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGc                                                                 <  2:2052043/150‑1 (MQ=255)
                                                                                            cccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATg                                                           <  2:4246565/150‑1 (MQ=255)
                                                                                            cccGCAGGCTTCAATGAATCGGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATg                                                           <  2:4246602/150‑1 (MQ=255)
                                                                                              cGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAACCCGCGAAAGGGCCGGCAAAACCTGCTCAGGGCGCAAACCCGTCCCCAATGCTCACCGCAGCCGTCGCCACCCCGTGGCGATGGCGGCCTGcc                                                         >  2:590862/1‑150 (MQ=255)
                                                                                               gCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAGCGGGTAAACCCCCGAAAGGGCCGGCAAATACTGCTCAGGGGGCAAAACCGTCCACAATGATCACCACAGCCTTAGCCACCCCGTGGCGATGGCTGCATGccg                                                        >  1:1484860/1‑149 (MQ=255)
                                                                                                  ggCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGGAAACCCCGAAAGGGCCCGCAAAACCGGCGCAGGGGGCAAACCCATCCAAAAAGATCACCTCAGCCGTCGCCCCCCCCTGGCGCTGGCGGCATGCctcat                                                     >  1:2527444/1‑148 (MQ=255)
                                                                                                   gCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                    <  2:687162/150‑1 (MQ=255)
                                                                                                        aaTGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACa                                               <  2:3079895/150‑1 (MQ=255)
                                                                                                         aTGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTTAACCCGCGTAAGGGGCGGCAAATACTGCGCAGGGCTCAAAACCGTCCACAAAGATGACCACAGCCATCGCCACCCCGTGGCGGTGGCTGCATGGCACATGGTGCAg                                              >  2:1941289/1‑150 (MQ=255)
                                                                                                                tGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAg                                       <  2:1466739/150‑1 (MQ=255)
                                                                                                                tGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAACCATCCACAATGATCACCTCAGCCTTCGCCACCCCATGGCGATGGCTGCATGCCGCGTGGTACAGCGCggcg                                       >  1:3359444/1‑148 (MQ=255)
                                                                                                                 gCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTACACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGc                                      <  2:1235239/150‑1 (MQ=255)
                                                                                                                  cccgCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCg                                     <  1:1031781/147‑1 (MQ=255)
                                                                                                                  ccAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCg                                     <  2:4435456/150‑1 (MQ=255)
                                                                                                                  ccAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCg                                     <  2:1648074/150‑1 (MQ=255)
                                                                                                                   cAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGAAACCCCCGCAACGGCCCGCCAATACAGCGCAGGGCTCAAACCCATCCACAAAGGTAACCGCAGCCATCGCCACCCCGTGGCGATGGCTGCATGCCGCGTGGTACCGCGCGGGGCgg                                    >  2:1262515/1‑150 (MQ=255)
                                                                                                                       ggTGCGGCAGGTCCGGGCCATTAGCGGGTAAAACCGCGAATGGGGCGGGAAATGCTTGTGATGGGTGAATACCTTGGACAATTATGACCTCAGGCTTCGCCACCCCTTGGCGATGGCTGGATTGCTGTTTGTTCAGCGGGGAGGGGTAAc                                >  2:1577131/1‑150 (MQ=255)
                                                                                                                       ggTGCGGCAGGTCCGGCCCATTAGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCAGCACAGGGCGCAAACCCTTCCACAATGATCACCTCAGCCGTCGCCCCCCCGTGGCGATGGGGGGATGCCTCGTGGTACAGCGCGGGGCGGGAAc                                >  2:1312659/1‑150 (MQ=255)
                                                                                                                        gTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACa                               <  2:3346984/150‑1 (MQ=255)
                                                                                                                         tGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTACCaa                              <  2:2322859/150‑2 (MQ=255)
                                                                                                                                 gTCCGGCCCTGGCGCGGGTACACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCACCAACCCGTCCACACTGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGc                      <  2:2562834/150‑1 (MQ=255)
                                                                                                                                    cGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTg                   >  2:1721659/1‑150 (MQ=255)
                                                                                                                                    cGGCCCTTTCGCGGGTAAACCCGAGAACGGGCCGCCAAATCCAGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCATGGCGATGGATGCATGCCTCTTGGTACAGCGCGGCGCGGTAACAGGCGACGGccag                   >  2:1721687/1‑148 (MQ=255)
                                                                                                                                     ggCCCTTTCGCGGGTAACCCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgc                  <  1:3558249/150‑1 (MQ=255)
                                                                                                                                      gCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTACCAGGCGACGGCCTgcg                 >  2:2723196/1‑150 (MQ=255)
                                                                                                                                       cccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc                <  1:3261959/150‑1 (MQ=255)
                                                                                                                                       cccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc                >  1:3703489/1‑150 (MQ=255)
                                                                                                                                        ccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg               >  1:1639561/1‑150 (MQ=255)
                                                                                                                                        ccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg               >  1:3169794/1‑150 (MQ=255)
                                                                                                                                        ccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg               >  1:1322305/1‑150 (MQ=255)
                                                                                                                                         cTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGt              >  1:3032231/1‑150 (MQ=255)
                                                                                                                                          tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACGGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa             >  2:931979/1‑150 (MQ=255)
                                                                                                                                          tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTGGGCGATGGCTGCCTGCCGCTTGGTACAGCGCGGAGCGGTACCAGGCGACGGCCGGCGCGtc             >  2:1271491/1‑149 (MQ=255)
                                                                                                                                          tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa             >  2:4039272/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGc                      >  1:2267896/1‑140 (MQ=255)
                                                                                                                                           ttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGc                      <  2:2267896/140‑1 (MQ=255)
                                                                                                                                           ttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTAc            >  2:3438494/1‑150 (MQ=255)
                                                                                                                                                  gTAAACCCGCGACAGGGCCGGCACATCCTGCTCAGGGCTCAAACCCTCCACACTGAGCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGaa    <  2:4099796/150‑1 (MQ=255)
                                                                                                                                                  gTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGa     <  2:4472621/150‑1 (MQ=255)
                                                                                                                                                  gTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGa     <  1:4307917/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAACT  >  NC_002947/4340311‑4340609

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